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Abstract 



Although cosmological observations suggest that the fluctuations of seed fields are almost Gaus- 
sian, the possibility of a small deviation of their fields from Gaussianity is widely discussed. The- 
oretically, there exist numerous inflationary scenarios which predict large and characteristic non- 
Gaussianities in the primordial perturbations. These model-dependent non-Gaussianities act as 
sources of the Cosmic Microwave Background (CMB) bispectrum; therefore, the analysis of the 
CMB bispectrum is very important and attractive in order to clarify the nature of the early Uni- 
verse. Currently, the impacts of the primordial non-Gaussianities in the scalar perturbations, where 
the rotational and parity invariances are kept, on the CMB bispectrum have been well-studied. 
However, for a complex treatment, the CMB bispectra generated from the non-Gaussianities, which 
originate from the vector- and tensor-mode perturbations and include the violation of the rotational 
or parity invariance, have never been considered in spite of the importance of this information. 

On the basis of our current studies [l}]7], this thesis provides the general formalism for the 
CMB bispectrum sourced by the non-Gaussianities not only in the scalar-mode perturbations but 
also in the vector- and tensor-mode perturbations. Applying this formalism, we calculate the 
CMB bispectrum from two scalars and a graviton correlation and that from primordial magnetic 
fields, and we then outline new constraints on these amplitudes. Furthermore, this formalism 
can be easily extended to the cases where the rotational or parity invariance is broken. We also 
compute the CMB bispectra from the non-Gaussianities of the curvature perturbations with a 
preferred direction and the graviton non-Gaussianities induced by the parity-violating Weyl cubic 
terms. We also present some unique impacts to the violation of these invariances on the CMB 
bispectrum. 
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List of Figures 

1.1 Evolution of the Universe. The events in the Universe are presented as the function 
of the time and energy scales. In addition, the several cosmological probes, which 
provide us with information about the structure and evolution of the Universe, are 
illustrated. Acronyms: Big bang nucleosynthesis (BBN), large scale structure (LSS), 
baryon acoustic oscillations (BAO), quasi-stellar objects = quasars (QSO), Lyman-alpha 
(Lya), cosmic microwave background (CMB), Type la supernovae (la), hydrogen 21cm- 
transition (21cm). This figure is quoted from Ref. [13] 

1.2 Balance between the scale of the density fluctuation and the co moving horizon. 
These are functions of the comoving scales |r| and the scale factors a. The comoving 
horizon evolves as r = Jdt'/ait') oc —a ^ (at the inflationary epoch), a (at the 
radiation dominated era), and a^^^ (at the matter dominated era). The scales of 
the density contrasts were larger than the Hubble radius until a ~ 10~^. However, 
before inflation operated, all interest scales of the density contrasts were smaller than 
the Hubble radius and therefore the physics is described in terms of the quantum 
mechanics. Similarly, at a very late time, the scales of the cosmological interests 
re-entered in the Hubble radius. This figure is quoted from Ref. [i] [15] 

2.1 Example of the potential of the inflaton 0. An accelerated expansion of the Universe 
occurs when V{(j)) ^ |(9f0)^. Inflation ends at = 0end when ^{dtcpY ~ V. The 
CMB fluctuations are generated from the quantum fluctuations about 60 e-folds 
before the end of inflation. At reheating, the energy density of the inflaton is 
converted into radiations. This figure is quoted in Ref. [l] [TH] 

2.2 Large-field inflation. This shape of the potential can be realized when V{(j)) oc (ff. In 
cases like this, the inflaton field moves over a super-Planckian range during inflation, 
namely, A0 > Mpi, and a large amplitude of the gravitational waves is produced. 
This figure is quoted in Ref. [l] M 

2.3 Small natural inflation. If the periodicity 27r/ is super-Planckian, the model can 
naturally induce a large amplitude of the gravitational waves. This figure is quoted 

in Ref. (l) M 

3.1 CMB anisotropy on the last scattering surface. The red (blue) parts correspond to 

the hot (cold) spots (Copyright 2011 by Daichi Kashino) [31] 

3.2 Interaction between several components in the Universe [3S] 

3.3 Geometry of Thomson scattering. Blue (Red) sohd and two dashed arrows de- 
note the incident (scattered) wave number vector and its orthogonal unit vectors, 
respectively. We set that x' = x [lO] 

3.4 Geometry for the line-of-sight direction \5U\ 

3.5 The CMB spectra of the // (left top panel), IE (right top one), EE (left bottom one) 
and BB (right bottom one) modes. Here we consider a power-law flat ACDM model 
and fix the tensor-to-scalar ratio as r = 0.1. The other cosmological parameters are 



fixed to the mean values reported in Ref. 27 

3.6 CMB // power spectra sourced from the scalar-mode perturbations in a power-law 
flat ACDM model. The red sohd line is plotted with = 0.112, = 0.728, Wf, = 



0.02249, = 0.967, k = 0.088 27 . The green dashed, blue dotted and magenta 
dot-dashed lines are calculated if Uc decreases to 0.1, fl^ increases to 0.8, and Ub 
increases to 0.028, respectively [60] 
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4.1 Representations of triangles forming the bispectrum, F^^^^^^ {ki, k2, k^), with various 
combinations of the wave numbers satisfying k^ < k2 < ki. This figure is adopted 



from Ref. 19 



4.2 Two dimensional color contour for the shapes of the primordial bispectra. Each 
panel describes the normalized amplitude of S{ki, ^2, kz){k2/kiY{k^/ki)'^ as a func- 
tion of /c2/fci and /cs/fci for k^ < k2 < ki. Since the primordial bispectra shown 
here are nearly scale invariant, the shapes may look similar regardless of the val- 
ues of ki. The amplitude is normalized by a rule that it is unity at the point 
where S{ki, k2, k^){k2/kif'{k^/kif' takes on the maximum value. The top-left panel 



is the local form given by Eq. (4.7), which peaks at the squeezed configuration. 



The bottom-left panel is the equilateral form given by Eq. (4.10), which peaks at 



the equilateral configuration. The top-right panel is the orthogonal form given by 



Eq. (4.11 ), which has a positive peak at the equilateral configuration, and a negative 
valley along the elongated configurations. Note that all of these shapes are almost 
orthogonal to each other. This is quoted in Ref. |TT] |69] 

5.1 The CMB /// (left top panel), IIE (right top one), I EE (left bottom one), and 
EEE (right bottom one) bispectra induced from the local-type (red solid line), 
equilateral-type (green dashed one), and orthogonal- type (blue dotted one) non- 
Gaussianities of curvature perturbations. The three multipoles are fixed as ^l = 
£2 = is = i. Here we consider a power-law flat ACDM model and fix the nonlinear 
parameters as f^l^^ = f^^^ = /nl*^"^ — 100- The other cosmological parameters 

are fixed to the mean values reported in Ref. |9] [76] 

5.2 The CMB /// (left top panel), IIE (right top one), I EE (left bottom one), and 
EEE (right bottom one) bispectra induced from the local-type (red solid line), 
equilateral-type (green dashed one), and orthogonal- type (blue dotted one) non- 
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6.2 Absolute values of the CMB reduced bispectra of temperature fluctuation for di = 
£2 = is = £■ The lines correspond to the spectra generated from tensor-scalar-scalar 



correlation given by Eq. (6.22) with gt^g = 5 (red solid line) and the primordial 
non-Gaussianity in the scalar curvature perturbations with /nl*^' = 5 (green dashed 
line). The other cosmological parameters are fixed to the mean values limited from 



WMAP-7yr data reported in Ref. 10 



6.3 Absolute values of the CMB reduced bispectra of temperature fluctuation generated 



from tensor-scalar-scalar correlation given by Eq. (6.22) {TSS + STS + SST) and 
the primordial non-Gaussianity in the scalar curvature perturbations {SSS) as a 
function of is with ii and £2 fixed to some values as indicated. The parameters are 
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1 INTRODUCTION 



1 Introduction 

1.1 Big Bang Universe 

The discovery of Bubble's law that light from distant galaxies is systematically shifted toward the 
red end of the spectrum supported the Friedmann Universe and guided us to the picture of an 
expanding Universe [l]. Next, the observed abundance of the light elements (H, He, and Li) and 
the CMB radiation match the predictions of the Big Bang scenario that the hot Universe cools 
gradually via the spatial expansion |2) |3). Table [lT| and Fig. [L1| show the major events in the 
thermal history of the Universe as a function of the time and energy scales. At this point, we will 
give a summary of these milestones in the evolution of the Universe on the basis of Ref. jlj. 

From 10~^°sec to today, the evolution of the Universe is largely-understood by the theories of 
the particle physics and gravity. Here, we start by discussing the Universe at lOOGeV, namely, 
the epoch of the electroweak phase transition (10~^°sec). Above lOOGeV, the electroweak sym- 
metry is restored and the Z and bosons behave as massless particles. Then, the interactions 
are strong enough to keep quarks and leptons in the thermal equilibrium. Below this energy 
scale, since the electroweak symmetry is broken, the Z and bosons acquire a finite mass, and 
the cross section of the weak interactions decreases as the Universe cools down. As a result, at 
around IMeV, neutrinos decouple from matters such as electrons. Promptly after, at Isec, the 
temperature of the Universe drops below the electron mass and there is efficient annihilation of 
electrons and positrons. Then, an initial matter-antimatter asymmetry of about 10~^ survives and 
the resulting photon-baryon fluid maintains a thermal equilibrium. Around O.lMeV, the nuclear 
binding processes become important, and the light elements (H, He, Li, and Be) are formed from 
protons and neutrons during the Big Bang nucleosynthesis (~ 200sec). The consistency with the 
observational abundances of H, He, and Li backs up the Big Bang theory. The energy density of 
matter is comparable to that of radiation around leV (lO^^sec). The charged particles, namely 
protons and electrons, and photons are strongly coupled in the plasma via the Coulomb and Thom- 
son scatterings, and the density contrasts propagate as the cosmic sound waves. Around O.leV 
(380000years), protons and electrons combine into neutral hydrogen (and helium) atoms. Then, 
photons decouple and become the free-streaming CMB. After 13.7 billion years, these photons 
are cooled to around 3K (300/ieV) because of the red shift and reveal the snapshot of the early 
Universe. The anisotropics in the CMB intensity and polarizations provide the evidence for the 
fluctuations in the density of the primordial seed flelds. 

These tiny density contrasts, namely 5 = [p(x, t) — p]/p ~ 10~^ with p and p being the energy 
density at an any point and its mean value respectively, grow via gravitational instability to form 
large-scale structures observed in the late-time Universe. A competition between the pressure 
of radiation and the gravitational attractive force affects the growth of the structure. During 
the radiation dominated era, the growth slows down, namely 5 oc Ina with a(t) being the scale 
factor that quantifles the space-time expansion. The clustering becomes more efficient after matter 
dominates over the background density (and the pressure vanishes), namely 5 oc a. At first, small- 
scale contrasts become nonlinear as (5 > 1, and form gravitationally bound objects that decouple 
from the overall expansion. This explains the hierarchical structure formation; fast, small-scale 
structures (like stars and galaxies) are formed and then they merge into larger structures (clusters 
and superclusters of galaxies). Around redshift z~25 (l + 2; = a~^), high energy photons arising 
from the first stars begin to ionize hydrogen in the intergalactic medium. This process, so called 
reionization, lasts until 2 6. Meanwhile, the most massive stars burn up the nuclear fuels 
and explode as supernovae. In these explosions, the heavy nucleons (C, O, etc) necessary for the 
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Figure 1.1: Evolution of the Universe. The events in the Universe are presented as the function 
of the time and energy scales. In addition, the several cosmological probes, which provide us 
with information about the structure and evolution of the Universe, are illustrated. Acronyms: Big 
bang nucleosynthesis (BBN), large scale structure (LSS), baryon acoustic oscillations (BAO), quasi-stellar 
objects = quasars (QSO), Lyman-alpha (Lya), cosmic microwave background (CMB), Type la supernovae 
(la), hydrogen 21cm-transition (21cm). This figure is quoted from Ref. j4j 

creation of life are generated. At z ~ 2, the second accelerated expansion starts. This may be 
because an unknown energy with negative pressure, the so called "dark energy", that starts to 
dominate the Universe. Further, the background space-time is accelerating and the growth of the 
structure ends as 5 ~ const. 

1.2 Paradigm in the early Universe 

The evolution of the Universe from 10~^°sec (ITeV) after its birth to today as discussed above, is 
supported by the observational facts and theoretical physics as the standard model of the particle 
physics and general relativity. Prior to this epoch, the temperature of the Universe exceeded ITeV 
and we have not had any direct experimental discoveries so far. Hence, the physics is as speculative 
as it is fascinating. 

To explain the CMB anisotropics, we require an input of the primordial seed fluctuations. We 
believe that these primordial fluctuations were generated in the very early Universe (~ 10~^'^sec) 
during the inflationary era. We will explain how the microscopic quantum fluctuations of fields 
get stretched by the inflationary expansion to the macroscopic fluctuations larger than the horizon 
at that time. After a metric perturbation exits the horizon, no causal physics can affect it, and 
it remains frozen with constant amplitude until it re-enters the horizon at a later time during the 
conventional (non-accelerating) Big Bang expansion^ The fluctuations associated with the cos- 



^This is valid only when there are no anisotropic stress fluctuations. 
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JNeutrmo decoupling 
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10^ - 10^ yrs 
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8 X 10^ yrs 
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We born 


14 X 10^ yrs 


100 fieV 






Table 1.1: Major events in the Universe. This is quoted from Ref. H 



mological structures re-enter the horizon when the Universe is about 100000 years old, right before 
the recombination (see Fig. 1.2). Inside the horizon, the perturbation amplitudes are changed by 
the causal physics, which means the CMB anisotropies, the galaxies, and clusters of galaxies are 
generated. Since we can calculate the evolution of perturbations after they re-enter the horizon, 
we can use the late time observations of the CMB and the LSS to infer the input spectrum of the 
primordial fluctuations. If we assume that this spectrum was produced by inflation, this leads to 
an observational probe of the physical conditions for the 10~'^^sec epoch. 



1.3 Concept of this thesis 

Conventionally, the information of the primordial density fluctuations has been extracted from the 
two-point functions (power spectra) of the CMB fluctuations. There is a statistical property that 
although a non-Gaussian variable generates both even and odd-point correlations; a Gaussian vari- 
able generates only even-point correlations. Hence, it is hard to discriminate between the Gaussian 
and non-Gaussian signals in the CMB power spectrum. Theoretically, whether the primordial seed 
fluctuations are Gaussian depends completely on the inflationary models. Therefore, the check of 
the non-Gaussianity of the primordial fluctuations will lead to a more precise comprehension of the 
early Universe. To extract the non-Gaussian signals from the CMB anisotropy, we should focus on 
the higher order correlations of the CMB fluctuations such as the CMB three-point correlations 
(bispectra). Owing to the recent precise observation of the Universe, the CMB bispectra are be- 
coming detectable quantities. As a result, the CMB bispectra are good measures of the primordial 
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Figure 1.2: Balance between the scale of the density fluctuation and the comoving horizon. These 
are functions of the comoving scales |r| and the scale factors a. The comoving horizon evolves as 
r = oc —a ^ (at the inflationary epoch), a (at the radiation dominated era), and a^'"^ 

(at the matter dominated era). The scales of the density contrasts were larger than the Hubble 
radius until a ~ 10~^. However, before inflation operated, all interest scales of the density contrasts 
were smaller than the Hubble radius and therefore the physics is described in terms of the quantum 
mechanics. Similarly, at a very late time, the scales of the cosmological interests re-entered in the 
Hubble radius. This flgure is quoted from Ref. (3] 



non- Gaussianity. 

The primordial non-Gaussianities originating from the scalar components and their effects on 
the CMB bispectrum have been well-studied (Refs. |5}|6]). However, for some situations the vector 
components (vorticities) and tensor ones (gravitational waves) also act as non-Gaussian sources. 
This indicates that unknown signals, unlike the scalar case, may also appear in the CMB bispectra. 
To study these impacts in detail, we produced the general formulae for the CMB temperature and 
polarization bispectra from the scalar, vector and tensor non-Gaussianities [7|[8]. Next, utilizing 
these formulae and computing the practical CMB bispectra, we obtained new constraints on some 
primordial non-Gaussian sources and learned more about the nature of the early Universe [9 -13 



This thesis aims to discuss the CMB bispectra induced by the primordial scalar, vector, and 
tensor non- Gaussianity on the basis of our recent studies (7 13 . More concrete organization of 



this thesis is as follows. In Sees. [2] and |3| we demonstrate how to generate the seed fluctuations 
in the inflationary era on the basis of some review papers and present formulae for the scalar, 
vector, and tensor modes of the CMB anisotropics as mentioned in Ref. [7j. We also review some 
observational flndings obtained by the analysis of the CMB power spectra. In Sees. |4] and [5} we 
describe the general formulae of the CMB bispectra generated from the primordial scalar, vector, 
and tensor non-Gaussianities [sj. We then discuss the applications to the non-Gaussianities in two 
scalars and a graviton correlator [s] (Sec. [6]), involving the violation of the rotational or parity 
12 13] (Sees. [7]and|8]), and sourced by the primordial magnetic flelds (9}jll] (Sec. |9]). 



mvariance 



Finally, summarize this thesis and discuss some future issues (Sec. 10). In the appendices, we 
describe some mathematical tools and the detailed calculations required for the conduct of our 
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2 Fluctuations in inflation 

Inflation expresses an exponential growth of tlie scale factor of the Universe in the early time, 
namely, a ~ e^*. In Einstein gravity, this requires p ~ —p with p and p being the pressure and 
energy density, and is often realized by the existence of a enigmatic scalar field, inflaton. We believe 
that the small fluctuations of this field have created the curvature perturbations and the density 
contrasts. Moreover, some vorticities and gravitational waves may also have evolved together. In 
this section, we describe the physical treatment of these fluctuations in the inflationary era in 
accordance with Ref. (ll. 



2.1 Dynamics of the scalar field 



S 



S-ER + Saj 



(2.1) 



In the simplest models of inflation, there exists a single scalar field 0, the so-called inflaton. The 
dynamics of a scalar field (minimally) coupled to gravity is derived from the action 

where R denotes the Ricci scalar and Mpi = (SttG)^^/^ is the reduced Planck mass. This action 
consists of the gravitational Einstein-Hilbert action, ^eh, and the action of a scalar field with 
canonical kinetic term, S^. The potential V{(j)) describes the self-interactions of the scalar field. 
The energy momentum tensor for a scalar field is 



2 5S. 



Then, the field equation of (f) is given by 

5Ss 1 



<p-g,A^d''<pd^<p + V{<p) 



-9 



where = dV/ dcf). Assuming the FLRW metric as 



ds^ 



-dt^ + a?dx^ 



-dr'^ + dx'^ 



(2.2) 



(2.3) 



(2.4) 



with T being the conformal time and the homogeneity of the inflaton, 0(t, x) = 0(t), its energy 
momentum tensor takes the form of a perfect fluid with 



^(9,0)^+^(0), p^ = liMr 



Thus, the parameter for the equation of state of 0: 



V 



(2.5) 



(2.6) 



shows that a scalar field can realize the negative pressure (w^ < 0) and accelerated expansion of 
the Universe [w^ < —1/3) when the potential energy V dominates over the kinetic one {dt4>)'^/2. 
The Friedmann and field equations are written as 



3M2 



-(9*0)2 + y(0) 



d 



. 



where H = dta/a is the Hubble parameter. If the potential has large values, 
significant Hubble friction from the term Hdt4>. 



(2.7) 



experiences the 
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S(j) 









■>CMB 



-'end 
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Figure 2.1: Example of the potential of the inflaton 0. An accelerated expansion of the Universe 
occurs when V"(0) S> |(c?t0)^. Inflation ends at = 0end when ~ V. The CMB fluctuations 

are generated from the quantum fluctuations 50 about 60 e-folds before the end of inflation. At 
reheating, the energy density of the inflaton is converted into radiations. This figure is quoted in 
Ref. Q 



2.2 Slow- roll inflation 

The acceleration equation of the Universe dominated by a homogeneous scalar field is given by 



6M2 



where 



(P0 + Sp^) = H\l - en) 



(2.8) 



The so-called Hubble slow-roll parameter en is defined as 

dtH d\nH 



(2.9) 



dN 



(2.10) 



where dN = Hdt with N being the e-folding number. When ejy < 1, an accelerated expansion 
occurs. For en — ^ 0, the de Sitter limit, namely, — )■ — p<^, realizes. Then, the potential energy 
dominates over the kinetic energy as 



{dtct^f « V(0) . 

Moreover, only if the second time derivative of is small enough 

dh 



df^ 



< \^Hdt(t)\AVs 



(2.11) 



(2.12) 



2.2 Slow-roll inflation 
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an accelerated expansion will be sustained for a sufliciently long period of time. This corresponds 
to the smallness of a second Hubble slow-roll parameter as 

where owing to \t]h\ <^ 1, the fractional change of en per e-fold is small. The slow-roll conditions, 
^H, \Vh\ ^ 1; can be interpreted as the conditions on the shape of the potential as 



<<^)^^[y) ' vi<P)^My-f. (2.14) 

In the slow-roll regime 

e, (2.15) 



Eq. (2.7) is rewritten as 



r.9 V 

- ^ - const. , 

■^^^^Pi (2.16) 



d 

3H ' 

and the former equation leads to the de Sitter space-time as 

a{t) ~ e^* . (2.17) 

The parameters e and rj are called the potential slow-roll parameters, and in the slow-roll approx- 
imation, the Hubble and potential slow-roll parameters are related as 

e/f ~ e , r]H^r]-e. (2.18) 

Inflation ends when the slow-roll conditions are broken as 

e//(0end) = 1 , e(0cnd) ~ 1 • (2.19) 



From Eq. (2.16), the e-folding number prior to the end of inflation is evaluated as 

'^end 



N{(f)) = In 



a 

tend r<f>cnd TJ V 

Hdt = / -7r-d<P ^ / —d<P . (2.20) 



Using the slow-roll parameters, this is expressed as 



To solve the horizon and flatness problems, it is necessary that the total e-folding number exceeds 
about 60, namely, 

iVt,t = ln^> 60. (2.22) 

'^start 
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The precise value depends on the energy scale of inflation and on the nature of reheating after 
inflation. The initial condition of the observed CMB anisotropies are created at N^mb ~ 40 — 60. 
Considering an integral constraint 

V = ^^-b ~ 40 - 60 , (2.23) 
•/-end ^piv2e 

one can calculate the corresponding fleld value 0cmb- 

Finally, for an example, let us present the slow-roll analysis of the simplest model of inflation: 
a single fleld inflation driven by a mass term as 

= ^mV' . (2.24) 



The slow-roll parameters are 



e(0)=r^(0) = 2(^)' . (2.25) 



Mpi 



For e, |?7| ^ 1, we need to consider super-Planckian values for the inflaton as 

> y2Mpi = 0end . (2.26) 

The fleld value of the inflaton is related with the e-folding number before inflation ends is given by 

iVW = ^-^. (2.27) 

The origins of the observed CMB anisotropies are created at 

0cmb = 2^/N~kMpi ~ 15Mpi . (2.28) 



2.3 Reheating 



After inflation ends, begins to oscillate around the minimum of the potential as shown in Fig. 2A_ 
During this phase of the coherent oscillations, behaves as a pressureless matter: 

dtp^ + 3Hp^ = . (2.29) 

If the coupling of the inflaton fleld to other particles exists, the inflaton energy decays in accordance 
with 

dtp^ + {3H + T^)p^ = , (2.30) 

where the decay rate depends on model-dependent physical processes. Eventually, the infla- 
tionary energy density is converted into the relativistic degrees of freedom and the hot Big Bang 
starts. We can flnd the detailed review in Ref. [21. 



2.4 Inflation models 
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v{<t>) 




0end 0CMB 



Figure 2.2: Large-field inflation. This shape of the potential can be realized when oc 0^. 

In cases like this, the inflaton field moves over a super-Planckian range during inflation, namely, 
> Mpi, and a large amplitude of the gravitational waves is produced. This figure is quoted in 
Ref. (l). 



2.4 Inflation models 

The nature of inflation is a mystery and there still remain some questions to be answered. How- 
ever, we believe that inflation has occurred at an enormous energy scale (maybe as high as ~ 10^^ 
GeV), far out of reach of the particle accelerators. Hence, any description of the inflationary era 
requires extrapolations of the known laws of physics, and until recently, a simple phenomenological 
parametrization of the inflationary dynamics was used. In this approach, a variety of the inflation- 
ary potential functions as described in Figs. 2A_ and have been discussed and the experimental 
predictions have been computed. As we will see in the next subsection, the details of the spectra 
of the primordial perturbations will depend on the precise shape of the inflaton potential. 



2.4.1 Single-field slow-roll inflation 

The inflationary model selection amounts to a specification of the action of the inflaton and the 
coupling to the gravity. In the above discussion, we have considered the simplest model, namely, 
the single-field slow-roll inflation, characterized by the action as 



S 



1 , .2 ^ 1 



(2.31) 



From the time when the origins of the CMB fluctuations were generated at 0cmb to the end of 
inflation at 0cnd; the dynamics of the inflaton is determined by the shape of V{(j)). The different 
properties for V{(j)) are classified by determining whether the change of the inflaton field, which is 
expressed as A0 = 0cmb — 0end, is large or small as measured in Planck units. 



Small-field inflation 
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In the small-field models, the change of the field is small (sub-Planckian), namely, < Mp\. 
This leads to a fact that the amplitude of the gravitational waves produced during inflation 
is too small to be detected. The potentials realizing such small-field evolution often arise in 
the mechanisms of the spontaneous symmetry breaking, where the field rolls off an unstable 



equilibrium toward a displaced vacuum as shown in Fig. 2A_ A simple example is the Higgs- 
like potential as 

2 



v{4>) = Vo 



1 - 



More generally, the small-field models can be locally approximated by 



= Vo 



1 



(2.32) 



(2.33) 



where the dots express the higher-order terms becoming important near the end of inflation 
and during reheating. A famous inflationary potential is the so-called Coleman- Weinberg 
potential as [3| |4] 

4 



V^(0) = Vo 



(2.34) 



which appears from the radiatively-induced symmetry breaking in the electroweak and grand 
unified theories. Although the original values of the parameters Vq and /i based on the 
SU{5) theory are incompatible with the small amplitude of the fluctuations in inflation, this 
potential remains a popular phenomenological model (e.g. Refs. Isl). 



Large-field inflation 

In the large-field models, the inflaton shifts from a large field values to a minimum of the 
potential = 0. If the field evolution is super-Planckian, A(f) > Mpi, the gravitational waves 
induced by inflation should be observed in the near future. 

The prototypical large-field model is the so-called chaotic inflation where a single monomial 



term dominates the potential as described in Fig. 2.2 



V{<j)) = Ap0^ . (2.35) 

For such a potential, the slow-roll parameters become small for super-Planckian field values, 
(p ^ Mpi. Note that the slow-roll conditions are independent of the coupling constant Xp. 
However, to arrange for a small amplitude of the density fluctuations, the inflaton self- 
coupling should be very small as Ap <^ L This condition automatically guarantees that the 
values of the potential are sub-Planckian as V <^ Mpj, and the quantum gravity effects are 
not necessarily important. 

The natural inflation, whose potential takes 



V(0) = Vo 



cos 



/ 



+ 1 



(2.36) 



is one of the most elegant inflationary models. This potential shown in Fig. 2.3 often appears 
if the inflaton field is taken to be an axion. Depending on the parameter /, 6(j) can become 
small or large. However, it is particularly attractive to consider the natural inflation for 
the large-field variations, 27r/ > Mpi, since for axions a shift symmetry can be employed to 
protect the potential from the correction terms even over the large field ranges. 



2.4 Inflation models 
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v{4>) 




Figure 2.3: Small natural inflation. If the periodicity 27r/ is super-Planckian, the model can 
naturally induce a large amplitude of the gravitational waves. This flgure is quoted in Ref. pi. 



2.4.2 Beyond the single-field slow-roll inflation 

There are a variety of possibilities for realizing the inflationary expansion. Inflation is often stud- 
ied as a framework for a theory of the early Universe, but it is not a unique theory. A lot of 
the phenomenological models have been proposed with the different theoretical motivations and 
observational predictions. In these models, the simplest inflationary action (2.31) is extended in 
several ways: 

• Non-minimal coupling to gravity 



The action (2.31) is so-called minimally coupled model in the sense that there is no direct 
coupling between the inflaton fleld and the metric. In contrast, we also can imagine a non- 
minimal coupling between the inflaton and graviton. However, in practice, the non-minimally 
coupled theories can be transformed to the minimally coupled form by a fleld redeflnition. 

Modified gravity 

Similarly, the possibilities that the Einstein-Hilbert action is modifled at the high energies 
can be entertained. However, the simplest examples for this UV modiflcation of gravity, so- 
called f{R) theories, can also be expressed by a minimally coupled scalar fleld with potential 
V{<P). 

Non-canonical kinetic term 



The action (2.31) has a canonical kinetic term 

V(0) 



£,A, — X 



X 



1 



(2.37) 



Then, inflation can only occur if the potential V{(j)) is very flat. However, we can imagine 
the high-energy theory with the non-canonical kinetic terms as 



C^ = F{<P,X)-V{<t>) 



(2.38) 
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where X) is the function composed of the inflaton field and its derivatives. For actions 
like this, inflation can be driven by the kinetic term and occur even in a steep potential. In 
Sec. [6} we treat the inflationary model with the non-canonical kinetic terms. 

• Multi-field inflation 

If there are two or more fields during inflation, then the possibilities for the inflationary dy- 
namics and evolutions the cosmological perturbations expand dramatically. Several physical 
properties of the multi-field inflationary models are reviewed in e.g., Ref. [6]. The inflationary 
model mentioned in Sec. [7] is one of the models of the multi-field inflation. 



2.5 Scalar perturbations 



Here, we focus on the full computation of the quantum-mechanical scalar fluctuations generated 
during inflation and their relation to cosmological perturbations. Our calculation follows closely 
the treatment by Maldacena |7j. 

We consider the single-field slow-roll inflation defined by 



(2.39) 



To fix time and spatial reparametrizations, we choose the comoving gauge for the metric gij and 
the scalar field 6 as 



'J J ' ^i"' J 







(2.40) 



In this gauge, the inflaton field is unperturbed and all scalar degrees of freedom are parametrized 
by the gauge- invariant comoving curvature perturbation TZ{t,x.). There is an important property 
that TZ remains constant outside the horizon if there exist no extra anisotropic stresses [^[^ Hence, 
we can restrict our computation to the correction functions of TZ at horizon crossing. 



2.5.1 Quadratic action 



One expands the action (2.39) to second order in TZ: 

2 



^(2) = 1 / dV 

Defining the Mukhanov variable 



H 



zTZ 



{dtTZf 



dtcf) 

I — 
H 



(2.41) 



(2.42) 



and expressing with the conformal time r lead to the action for a canonically normalized scalar as 



drdx 



(2.43) 



^On superhorizon scales, this 7?. is consistent with TZ in Refs. fsl'O^, C in Refs. [7j[To|, —7?. in Refs. [T]|Tl], and 
— C in Ref. ^12j. In a numerical code CAMB y^[l4j, the primordial scalar-mode power spectrum is given by this 
Tl. 

^In Sec. [oj we will show that due to the finite anisotropic stresses of the primordial magnetic field, the curvature 
perturbations (and gravitational waves) do not remain constant even on the superhorizon 
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where ' = 8^- Performing the Fourier expansion of v as 



(2vr)' 



rfkirie 



ikx 



we can derive the field equation 



Wk + - - h^k = 



(2.44) 



(2.45) 



The Mukhanov equation (2.45) is hard to solve in full generality because the function z depends 
on the background dynamics. Then, one may find the numerical solution. However, in the pure 
de Sitter space-time, we can obtain the analytic solution. 



2.5.2 Quantization and initial condition 



One can perform the quantization of the field v in completely analogy with the treatment of the 
quantum harmonic oscillator. 

At first, we express the field v and its conjugate momentum v with the quantum operator as 



V ^ V 



rf^k 

(27r)3 
(2^ 



Vkirjak + vJrja 



„ikx 



This means that the Fourier component fk is rewritten as the operator: 



fk t^k = VkMcLk + vJrja 



-k ' 



(2.46) 



(2.47) 



where the creation (a|^) and annihilation (cik) operators satisfy the canonical commutation relation 
as 



Ok, a|^, 

if the mode functions are normalized as 

{Vk,Vk) = 

Note that these relations arise from 



h 



{27Tf6{k - k') 



ih . 



1 . 



(2.48) 



(2.49) 



(2.50) 



Equation (2.49 ) is used as one of the boundary conditions on the solutions of Eq. (2.45 ). The second 
boundary condition that fixes the mode functions completely comes from the vacuum selection. 
We set a vacuum state for the fiuctuations: 



ak |0) = . 



(2.51) 



As the vacuum, one often choose the so-called Bunch-Davies Vacuum denoting the Minkowski 
vacuum of a comoving observer in the far past, when all comoving scales were far inside the 
Hubble horizon, i.e., r — t- — oo or |A;r| ^ 1. In this limit, the mode equation (2.45) reaches 

Vk + k\k = . (2.52) 



26 



2 FLUCTUATIONS IN INFLATION 



This is identical to the equation of a simple harmonic oscillator with the time-independent fre- 
quency. Therefore, we can impose the initial condition as 



-ikr 



lim Vk 



2k 



(2.53) 



Owing to the two boundary conditions (2.49) and (2.53), we completely fix the mode functions on 
all scales. 



2.5.3 Solution in de Sitter space-time 

Considering the de Sitter limit: e — > {H = const.) and 



z a 



z a T 



2 ' 



the mode equation reduces to 



Vk+ [k'^ r Ufc = . 



(2.54) 



(2.55) 



Taking into account the two boundary conditions (2.49) and (2.53), one can gain this solution as 

-ikr 



e I 
V2k V 

Using this, we can express the time evolution of the primordial curvature perturbation as 



7^(k, 



r 



Hvk H 



adfd) adfi 



[vk{T)di, + vl{T)al^^] 



(2.56) 



(2.57) 



2.5.4 Power spectrum of curvature perturbations 

We then compute the power spectrum of curvature perturbations: 



\n=l 



\n=l 

= (2vr)35(^k„ 

\n=l 

On superhorizon scales, |/cr| <^ 1, this becomes 

J]7^(k„)\ = {2'nfPn{ki)5 (j^k^ 



^l^fci(^) 
adfd) 



2kf \d, 



[1 + ihrf 



(2.58) 



\n=l 



vn=l 



( H, 



2k^ \{d, 



HV * 



(2.59) 



Here we evaluate all quantities at horizon crossing, namely = —1/k. Note that since TZ = const, 
on superhorizon scales, this power spectrum becomes the initial condition for the CMB power 
spectrum of the scalar mode when a given fluctuation mode re-enters. 



2.6 Tensor perturbations 



27 



2.6 Tensor perturbations 

As the quantization of the scalar perturbation have been discussed in detail, the corresponding 
calculation for tensor perturbations will appear almost trivial. 

Expanding the Einstein-Hilbert action, one can obtain the second-order action for tensor fluc- 
tuations as 



5(2) 



pi 



drdPxa^ 



hijhij — dihijdih. 



(2.60) 



Except a normalization factor of Mpi/2, this is the same as the action for a massless scalar field in 
the FLRW Universe. We obey the Fourier expansion as 



(27r)3 



5^/.W(k,r)eS;Hk)e> 



kx 



(2.61) 



A=±2 



where e-^"* is the transverse-traceless polarization tensor which has two circular states A = ±2. The 



convention and useful properties of this tensor are described in Appendix |Dj Taking into account 



the normalization: e-^''(k)ej-^ ■*(— k) = 2(5a,a', the tensor quadratic action becomes 



5(^)= |rfr(^)Y(0^[l^^'^^'^'^)l'-^>^'^(>^'^)l 



Defining the canonically normalized field: 



pi 



V2 



we have 



A=±2 



d^k 1 

(27r)3 2 



,W|2 



dr 



a, (A) 
-\vl 



(2.62) 



(2.63) 



(2.64) 



Here, the third term in the square bracket is the topological term, hence this never affects the 
equation of motion. Then, through 5S/5v^'^^ = 0, we obtain the mode equation of the gravitational 
wave: 



. 



(2.65) 



This is consistent with the equation of the scalar mode (2.55), therefore, the pure de Sitter solution 
is given by Eq. (2.56 ). Hence, in the de Sitter limit, the quantized mode function of the gravitational 
wave becomes 



h'^'\k,T) = h{r)dl'^ + hl{T)dii^ 



where d^-^'' and d^'^^'^ satisfy 



hkir) 



4'^ |0) = 



Mpi 



(1 + ikr) 



(2.66) 



"k ' "k' 



{2nf6{k - k')5x, 



(2.67) 



A' • 



28 



2 FLUCTUATIONS IN INFLATION 



Then, the power spectrum of primordial gravitational waves on the superhorizon limit, namely, 
\kT\ <^ 1, is calculated as 



n /^^'"Hkn) ) = ( 



n=l 



(2.68) 



,n=l 



Note that like the curvature perturbation, the gravitational wave also becomes constant outside 
the horizon if there are no anisotropic stress fluctuations. 



2.7 Consistency relations in the slow roll limit 



As a measure of the amplitude of the primordial gravitational wave, one often use the tensor-to- 
scalar ratio as 



r = 



2Ph{k) 
Pnik) 



Comparing Eq. (2.59) with Eq. (2.68), we find a consistency relation 



16e 



H* 



16e* 



(2.69) 



(2.70) 



Note that using Hdt = dN, the tensor-to-scalar ratio r relates directly to the evolution of the 
infiaton function of the e-folding number N 



(2.71) 



M2 \dN 



The total field evolution between A'cmb and A^^cnd can therefore be written as 



A(f) 

w 



pi 



dN 



(2.72) 



During the slow- roll inflation, r{N) does not evolve much and one may gain the following approx- 



imate relation, so called, the Lyth bound 15 



pi 



0.01 



1/2 



(2.73) 



Therefore, the large values of the tensor-to-scalar ratio, r > 0.01, imply A0 > Mpi, namely, a 
consequence of the large-field inflation. 

The scale dependence of the spectra is quantified by the spectral indices as 



dlnPn 
dink 



Tit - 3 = 



dlnPh 
dink 



(2.74) 
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From Eqs. (2.59) and (2.68), the right-hand sides are expanded as 

dlnH^: din en 



d\n.k 
dlnPh 
dink 



dN 
dlnH^ dN 
dN dink 



dN 
-3 . 



dN 
dink 



-3 



(2.75) 



From the definition of the Hubble slow-roll parameters (2.10) and (2.13), we obtain 

dlnH^ 



d In en* 



dN '^^ ' dN 

The factor dN/{dlnk) is evaluated by recalling the horizon crossing condition k = a^^H^, or 



(2.76) 



dink = dN + dlnH^ , 



hence we have 



dN 
dink 



dink 

dN 



-\ -1 



1 + 



dlnH^: 
dN 



-\ -1 



1 + gh* 



(2.77) 



(2.78) 



As the result, to first order in the Hubble slow-roll parameter, we find the consistency relations: 



Us - I 

nt 



2r] 



H* 



-2e 



H* 



^ -2e, . 



Then, from (2.70), we also find the consistency relation between r and nt as 

r = —Srit . 



(2.79) 



(2.80) 



In the slow-roll limit, the measurements of the scalar and tensor spectra relate directly to the 
shape of the potential V^(0). This means that if is a measure of of the scale of the potential, 
e of its first derivative V^, r] of its second derivative V^,^, etc. Therefore, the measurements of 
the amplitude and the scale-dependence of the cosmological perturbations encode the information 
about the potential driving the infiationary expansion. This allows to reconstruct a power series 
expansion of the potential around </)cmb corresponding to the time when the origins of the CMB 
fiuctuations exited the horizon. 
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Figure 3.1: CMB anisotropy on the last scattering surface. The red (blue) parts correspond to the 
hot (cold) spots (Copyright 2011 by Daichi Kashino). 



3 Fluctuations in cosmic microwave background radiation 

Cosmic microwave background (CMB) radiation is composed of photons which have decoupled from 
electrons in the epoch of the hydrogen and helium recombination at 2; = 1089 and it is observed 
as the perfectly black body radiation whose averaged temperature is 2.725K. Historically, in 1949, 
Alpher and Herman predicted its existence as relics of the big bang Universe and its first detection 
came in 1964. More precisely, however, the CMB involves the spatial fluctuations of O(10~^)K 



(see Fig. 3.1). We had to wait the detection of the CMB anisotropy until the data of the COBE 
experiment were released in the 1990s. 

Theoretically, the density contrast of the CMB is computed in the system where photons, 
neutrinos, baryons, dark matters and dark energy exist in the gravitational potential. Compared 
with the observational data, the values of several key parameters have been well-determined. The 
WMAP experiment established the facts that the Universe is close to spatially flat and the present 
structure grew from the nearly scale-invariant primordial fluctuations. These consequences are 
almost consistent with the prediction of the standard slow-roll inflation. Furthermore, we have a 
compelling evidence that the Universe is dominated by dark energy and dark matter, which implies 
that 96% of the total energy of the Universe remains unknown. Nowadays, some bare anomalies 
such as the preferred direction and the parity violation are furthermore being discussed [l], and 
we expect to extract more detailed information from the new precise measurements [2]. 

In addition to the intensity of the CMB, the polarizations also lead to better understandings. 
The curl-free component of the polarizations, E mode, reflects the recombination history, in par- 
ticular, the reionization of the Universe. The curl component of the polarizations, B mode, is 
generated from the primordial vector and tensor perturbations. Hence, the detection of the B- 
mode polarization will provide clues as to inflation and the physics beyond the standard model of 
the particle physics. 
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3 FLUCTUATIONS IN COSMIC MICROWAVE BACKGROUND RADIATION 



In this section, we describe the original formalism of CMB fluctuations including intensity and 
polarization anisotropies from the scalar, vector, and tensor modes partially on the basis of our 
paper [s], some publications [4]-[7] and some academic dissertations [8 10 , and summarize current 
outputs from the analysis of the CMB power spectra. 

3.1 Einstein equations 

Here, we derive the zeroth and first-order Einstein equations. Let us consider the flat (K = 
0) FLRW metric and small perturbations in the synchronous gauge (for open and closed cases, 
see [njg): 

ds^ = a^[-dT^ + {5ij + hij)dx'dx^] . (3.1) 
We have the inverse metric to first order in perturbations as 

9'' = --,. ^7°' = 0, (f^ = ^{S^^-h^r)_ (3.2) 
The Einstein equation with the cosmological constant A can be written as 

= - ^(5^/? = SttGT^ - A5^ , (3.3) 

where the left-hand side denotes the curvature of space-time and the right-hand one is the energy 
momentum tensor. The Ricci tensor R^y and Ricci scalar, namely a contracted form of the Ricci 
tensor, R, are expressed with the Christoffel symbols as 

p pa pa pa , pa p/3 pa p/3 /o A\ 

where = da- The Christoffel symbols in a metric space without torsion are given by 

r = 2^ (diMK,!^ + 9i^K,iJ. ^ 9ij,v,k) ■ (3-5) 
Up to first order, we can express as 

r°oo = ! 



r*nn — r'^in — 



r,k = l:{dkh',+djh\-d'h,^ 



2 ' (3.6) 

I;. 



k) 



therefore each component of the Ricci tensor is calculated as 
a'R% = 3 - + ^ {h\ + nh\) , 



(3.7) 



3.1 Einstein equations 
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Here, l-L = a/a = aH is the Hubble parameter in terms of conformal time with H being the 
observable Hubble parameter. Then the Ricci scalar is also given by 

i? = = 4j f 6- + h\ + mh\ - v^h\ + d^d^K'^ . (3.8) 



a^ \ a 



Contracting the Einstein equation (3.3) allows one to eliminate the Ricci scalar and reduce the 
Einstein equation to 

i?^ = SttG fr^ - J5^T%) + . (3.9) 



2 



Hence, in vacuum, we have = 0. 



3.1.1 Homogeneous contribution 



At zeroth order, the 00 and ii components of Eq. (3.3) lead to the Friedmann constraint equation 
and the Raychaudhuri evolution equation, respectively. Substituting Eq. (3.7) into Eq. (3.9), these 
are obtained as 



o''^ ° ^ (3.10) 

2--n' = -—a'f\ + a'A. 
a 3 



The physical meaning of these equations can be illustrated with the perfect fluid form as follows. 
The energy momentum tensor of the perfect fluid is given by 

= (p + p)m'^m.+p5^ , (3.11) 

hence the above equations change to 

, ,2 SttG 2 / A 



' \ (3.12) 

2^-n' = -8nGa^p 



a \ 8nG 

Note that we may identify the cosmological constant as a component of the perfect fluid as 

= ^diag(l, -1, -1, -1) . (3.13) 

To use a different phrase, an unperturbed perfect fluid of density and pressure are given by pa = 
A/{S7iG),pa = — Pa- If we use w = p/p, then w\ = —1. 
For convenience, we change the Friedmann equation to 

SvrG 2- 87rG'_ ^ 8nG _ 

i 

where in third equality, we decompose the total energy density in the Universe into individual 
species pj. Introducing a quantity which means the ratio between the energy density of each 
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w 


p(a) 


a(t) 


a(r) 




rad dom 


1/3 




^1/2 


r 





mat dom 







^2/3 







curv dom 






t 




— oo 


A dom 


-1 


a" 






— oo 



Table 3.1: FLRW solutions dominated by radiation, matter, curvature, or a cosmological constant. 



species and the critical density in the Universe at the present time, ilj, and which is expressed as 
f2j = 8nGpio/{3HQ), and using the scaling relation as pi = pjo/a"S this equation is rewritten as 



H 



(3.15) 



For radiations, matters and cosmological constant, we have ra, = 4, 3, 0, respectively. In this 
notation, we can also include a curvature term as a component of rii = 2. In Table |3.1[ we 



summarize the solutions of Eq. (3.15) if the cosmological fluid consists of a single component. 



3.1.2 Perturbed contribution 



At first order, 00 and ij components of Eq. (3.9) generate the evolution equations as 



1 



and 00 and iO components of Eq. (3.3) generate the constraint equations as 



2nh\ + d^dih'j - V^h\ 



■leTrGa^^T^ 



(3.16) 



(3.17) 



From here, let us express these equations with the variables in the helicity states. To do it, we 
decompose all kind of vectors and tensors, such as metric, velocities and energy momentum tensors, 
into each helicity part in accordance with the formulae: 




A=±l / 



(3.18) 



A=±l 



A=±2 



3.1 Einstein equations 
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where we define the projection vectors and tensors as 



O 



(±2). fx _ (±2), 



(3.19) 



or(k)-er(k) 



The polarization vector and tensor, e['^^\e\f'^\ satisfy the divergenceless and transverse-traceless 
conditions as 



4er)(k) = Mr(k) = er^(k)=0. 



J±2) 



(±2), 



(3.20) 



The prescription for the scalar-vector-tensor decomposition and explicit forms of the polarization 
vector and tensor are presented in Appendix D Then, from Eqs. (3.16) and (3.17), we can rewrite 
the evolution equations as 

kso + H/iiso = -SnGa'' {5T% + 5Tr) , 

and the constraint equations as 

UK, + h\h,, - /i(°)) = ^7,GaHT\ , 



(3.21) 



kill 



247rGa25rj°) , 



Here, we have obeyed the convention as 



(3.22) 



levrGa^^Ti^i) 




A=±l 



(3.23) 



A=±l 



A=±2 



Our perturbation quantities are related to the variables for the scalar mode in the synchronous 
gauge of Ref. p3], namely h and r], as 



(3.24) 



Hence, we understand the correspondence to the gauge-invariant variables by Bardeen ($^, ^h) 
[Is] and Kodama-Sasaki (^E', $) [Ig]: 



$j7 = $ 



2A;2 



2k 



(3.25) 
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3.2 Boltzmann equations 

The distribution function of several species evolves in accordance with the Boltzmann equation as 

dr dr dx'^ dr dpf^ dr \dr J ^ ^ 

where r is the conformal time, is the proper momentum of species, and a subscript C denotes 
the collision term. In this Boltzmann equation, there exist two contributions: the gravitational 
redshift and the effect of scattering, which correspond to the third term of the first equality and the 
term of the second equality, respectively. For convenience, we introduce the comoving momentum 
and energy as = ap\e = a^pp' + m^. Setting a unit vector parallel to the fluid momentum as 
q = gn and expanding the distribution function up to first order: 

/(x,g,n,r) = /(°)(g)[l + /(^)(x,g,n,T)] , (3.27) 

the above Boltzmann equation is rewritten as 

^/ .(0) (Qf^^ . df'^ dx^ 



dr V dr dx^ dr 



^?Z5!!^£^/(0)?£!!^i+/a.^^^+/(.,^^.(|/) . (3.28) 
oq dr oq dr oq dr on^ dr \ot J ^ 



To estimate dq/dr^dn" /dr, we consider the geodesic equation as 



/7 PM 

P°— + V^pP'^P^ = , (3.29) 
or 

where P'^ is the canonical momentum as 

The contraction is given by P^P^j, = — (e/a)^ = —w?. From // = component, we obtain 

dT~ 2^'''' ar • ^^-"^^^ 
Similarly, from the spatial components, we have 

2-^ = h'hmhn^ %^ ^-hmhnd^'h''' + ^-hmhnd'h^'' ^ 0{h) . (3.32) 

dr or or e e 

Furthermore, since we have the zeroth order expression as dx'^/dr — {q/e)n\ the Boltzmann 
equations up to first order are expressed as 

df^ ^ q.^df^ _ 1 .,.,.5^ain/(°) _ J_ fdn 
dr ^ e"" dx^ 2'''' dr dlnq ' f(^) [dr J ^ ' ^'^''^'^^ 

The general expression for the energy momentum tensor is given by 

/pfj.p 
V^d'^^f ■ (3-34) 



3.3 Stokes parameters 
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Substituting the relations: 



1 



-g = a 



1 - -h' 



(3.35) 



and noting that 

J fl^hjdVtn = ' J ^id^n = J flihjflkdVtn = , (3.36) 

the homogeneous and hnearized components of the energy momentum tensor are obtained as 

ef(^\l + fW)q^dqdQ^ 

ef^'k'dq - ^ / / ef'^''>f^'k'dqdQr^ = T% + 5T% , 







-P = 

4tt 




T\ = - 



a'* 
1 




1 

An 
3^ 





q 



(3.37) 



qn^fWfWq^dqdn^ = 5T\ , 
2 / 1. ..A /_ 1 



n 



n, + -n"hj,] + f^^)qHqdn^ 

q' 



S'jJ \f'h'dq+^^ J J '-n^n,f^'^f'^'\'dqdn^^T,+6T 



These components correspond to the density contrast 6, velocity f * and anisotropic stress 11*^ of 
fluid asQ 

6T% . 6T'n 5T^i 



-P 



-6 



P 



5T\ ^ 

p+p P+P 



-v' , Tj = p6'' 



3 ' 



p 



(3.38) 



Therefore, equating the integral of Eq. (3.33) over qi,q2,q3 with Eqs. (3.37) and (3.38), we can 
see that the Boltzmann equation (3.33) becomes the differential equations with respect to 6, f*, 



for each species. These equations correspond to the Euler and continuity equations. Generally, as 
the species of the cosmological fluid which mainly generate the inhomogeneity of the cosmological 
structure, there exist baryon, photon, neutrino and cold dark matter (CDM), hence we can trace the 
evolution of their fluctuations due to solving these Boltzmann equations coupled with the Einstein 



equations (3.21) and (3.22). Between baryons and photons, Thomson scattering is effective, so 



that their Boltzmann equations have the collision term. On the other hand, for neutrinos and 
CDMs, since there are no short-length interactions, the right-hand side of the Boltzmann equations 



vanishes. All these species couple with the metric via gravity. This relation is illustrated in Fig. 3.2 



3.3 Stokes parameters 

Here, we introduce the Stokes parameters to characterize the polarization states of the radiation 
field. For simplicity, at first, we consider a plane electromagnetic wave propagating along the z 
axis. The Fourier decomposition of the radiation field is expressed as 



/oo 
dk {±E^e'^^ +yEye''^y) 
-oo 



^i{kz—u)t) 



(3.39) 



^The anisotropic stress of the magnetic field is often normalized by photon's energy density as Eq. (9.2 1 
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% 



proton 



dark matter 



Figure 3.2: Interaction between several components in the Universe. 



where E^, Ey and (px, (py are the real quantities describing the amplitudes and phases in the x — y 
plane, respectively, and u = kc denotes the frequency of the wave. 
The Stokes parameters are given by 



I = \Ex\'^ + l-E'yP 

Q = \E!x\'^ — \Ey\'^ 
U ^ -2Re[E:Ey] 
V ^ -2lm[E:Ey] 



(3.40) 



where these are all real quantities. For the monochromatic wave, P = Q"^ + U"^ + V"^ is satisfied. 
/ measures the intensity of radiation and is always positive. The other parameters represent the 
polarization states and can take ether positive or negative values. Q and U quantify the magnitude 
of the linear polarization, and V parametrizes the circular polarization. While I and Q are parity- 
even quantities, U and V are parity-odd ones. 

In order to see the transformation rule of Q and U under rotation of axes, let us introduce the 
new coordinate {x',y'), which is related to the original coordinate {x,y) by the rotation around 
the z axis as 



x' \ f cosip sinip \ f X 
y' J \ — sin ip cos ip J \ y 



(3.41) 



3.4 Boltzmann equations for photons 
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Then, the radiation field is converted into 

E(z, t)= f dkf idE'e'^- + y'E'e'^A ^i{kz-ut) 

7-00 V " y J (3 42) 

{E',±tE'y) =e^'^{E^±iEy) ■ 

hence we have 

Q' ±iU' = e^^'^{Q±iU) . (3.43) 

This imphes that the hnear combination, Q ± iU , are the spin-±2 quantities. Therefore, the 
anisotropy of the hnear polarization should be expanded with the spin-2 spherical harmonics. 



3.4 Boltzmann equations for photons 

Here, for quantifying the CMB anisotropy, let us focus on the linearized Boltzmann equation for 
photons. The distribution function of photons is given by 



/ 



exp 



P 



r[i + e(x,n,t)] 



- 1 



(3.44) 



hence we have 14 



/ 



(0) 



x,g,n,rj 



_q_df^ 
■/(o) dq 



e 



ain/(o) 
c^lng 



(3.45) 



where 6 = AT/T. Substituting Eq. (3.45) into Eq. (3.33), the Boltzmann equation in terms of 
is expressed as 



Q- 



df{0) /QQ ^^QQ 



dq \ dr 



n 



: dhjj 

~d7 



dl 
dr 



(3.46) 



c 



The first two terms in the bracket denote the free-streaming of the photon, whereas the remaining 
third term in the bracket account for the gravitational redshift. As for the CMB polarization, 
there exist no gravitational effects because Q, U and V themselves are first-order quantities and 
the third term in the left-hand side of Eq. (3.33) does not appear. Hence, we summarize the 
Boltzmann equations oiQ,Q,U and V: 



dr dx^ 



dU ^,dU 

— h n 1— 

OT ax* 

dV ^,dV 
-T, \- n T— 

OT (7X* 



(3.47) 



Ut 
Vt 



where G-r, Qt, Ut and Vr denote the collision terms of Thomson scattering and ' = d/dr is the 
derivative with respect to the conformal time. Next we consider the contribution of these terms. 
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Figure 3.3: Geometry of Thomson scattering. Blue (Red) solid and two dashed arrows denote the 
incident (scattered) wave number vector and its orthogonal unit vectors, respectively. We set that 
x' = X 



3.5 Thomson scattering 

The process of scattering off a photon by a charged particle without the energy exchange of photons 
is called the Rayleigh scattering. In particular, when the charged particle is an electron, the 
process is known as Thomson scattering. During the epoch of recombination, electrons scattered 
off photons by Thomson scattering. Here, we consider an incoming plane wave of the radiation 
with a wave number vector kj parallel to the z axis and an outgoing radiation scattered off by an 
electron with a wave number vector ks. We take the plane spanned by ki and kg as the scattering 



plane as shown in Fig. 3.3 



The differential cross section of Thomson scattering is given by 

|ki-ksP, (3.48) 



do 3o"T'|,'- 2 



(M Stt 

where d£L = d{cos6)d(j) and ax is the cross section of Thomson scattering. This equation quantifies 
the change of the intensity by the scattering. For simplicity, we consider the case of the x' axis 
parallel to the x axis. Here, we suppose the incident radiation with the polarization states /' = 
{I'y,, I'^,,U' ,V'), where /' = + 1'^, and Q' = 1'^, — I'^,. When there is no dependence on the 



azimuthal angle, namely 0' = 0, from the notation of Stokes parameters (3.40) and (3.48), we 
obtain 



/, = ^ cos^ er, , 4 = , u='^cos eu' , (3.49) 



where we have normalized these equations so that the number of photons is conserved during a 
single scattering. For a general case with a non- vanishing azimuthal angle 0, Q' and U' are replaced 
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as 



Q' ± lU' ^ e^^''^{Q' ± iU') . (3.50) 

Then, the changes of the Stokes parameters between the incident radiation from n' = z [9' = 
0, (f)' = 0), and the scattered radiation with h = {6 



Aefn' = z,n) 



A{Q±iU){n' = z,n) 



1 

An 
1 

4:71 



are given by 17 
3 



1 + cos2 e)e' ^^'""^ (fi' + i'^' 



s=±2 



— sm' 
4 



s=±2 



(3.51) 



where we use = A////4. Using the exphcit formulae of the spin-0 and spin-2 spherical harmonics 
described in Table A. 2, we can extend this expression to the form corresponding to an arbitrary 
direction of n': 



Ae(n',n) = 5^[|l 



Y2Un)Y*Jn') + Yor.ih)YLin') } Q' 



- E ^ Vf ^2„.(n).r4(nO (q' + '-tU') 

s=±2 



(3.52) 



(Ag±iAf/)(n',n) = ^ 



10 



Yo^(h) 



We will use Eq. (3.52) in the frame satisfying k || z, where k is the wave number vector. Integrating 
these equations over all directions n', we express the scattered fields as 



9r(n) = 



Q±tU) (n) 

T 



-K 



Q{h) - J dQ'AQ{n',h) -^ry,■h , 
{Q±iU){n)- I dn\AQ±iAU){n',h) 



(3.53) 



where we define the differential optical depth as k = aaxrieXe with UeXe being the density of ionized 
electrons, and its total value at time r is given by 



K r = 



k{T')dT' 



(3.54) 



with To being the present conformal time. 

From here, we discuss the polarization property in more detail. One of the key points in 
Eq. (3.52) is that the temperature anisotropy generates the polarization of the CMB photons. 
Then, what mode of the temperature anisotropy is related to the generation of the polarization? 
For simplicity, we suppose that the incident radiation field is unpolarized, Q' = U' = V = and 
consider the case for n = z. Integrating Eq. (3.52) over all incident radiation, we gain 

3(Tt / 271 



Q±iU{z) 



47r V 15 



221 



(3.55) 



^The Stokes parameter, V , which means the circular polarization of photon, can be ignored because it cannot 
be generated through Thomson scattering if this is initially absent. 
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When the incident temperature (intensity) anisotropy is expanded with the spherical harmonics 



as G'(^',0') = T.im^'imYim{e\<P'), Eq. (3.55) is replaced with 




g±2f/(z) = -^W— a'22. (3.56) 



Thus, if there exists no quadrupole moment (£ = 2) in the unpolarized radiation field, the total 
scattered radiation along the z direction would be never polarized. Long before recombination, in 
the thermal equilibrium, the polarization states of photons are equally populated and the incident 
radiation should not have any polarization. Therefore, there are only the unpolarized radiations 
before recombination. Allowing the polarization at the last scatters just before the photons begin 
to stream freely, the polarized emission can lead to the multipole anisotropy and one has polarized 
quadrupole and octupole and so on. These effects are automatically involved in the Boltzmann 
equation. 



3.6 Transfer functions 



Here, we derive the CMB anisotropy sourced from scalar-, vector- and tensor-mode perturbations. 
We obey a Fourier transformation as 



X(n,r) 



(2^ 



Ax(r,k,n) , 



(3.57) 



where X = Q,Q±iU and Ax is called the transfer function. Note that we include the factor e**^'" 
in Ax. 

At first, for convenience, we derive the transfer functions of photons when k || z. For k || z, the 
scalar- vector-tensor decomposition of the gravitational redshift and Doppler term in the Boltzmann 
equation of photons (3.47) and (3.53) are given by 

1 „ . . 1 r 1.- 



-n^n^hij(k \\ z, r) 



/iiso(k II z, r) + /i(°)(k II z, r) - cos^ ek,n + 



+ J2 ^sin^^,,„cos^^fc,„e^^'^'=-/iW(k || z,r) 



A=±l 



^sin2^fc,„e^^*'=-/iW(k II z,r) 



A=±2 



= e(°)(k 



1 • 



h\;^ik,T) + h^'\k,r) (cos^^^fc,„- J 



+ sin^fc,„cos^fc,„e^^'^'=-e^^)(k || z)/i(^)(A;, r) 

A=±l 

+ Y sm^4,ne"*'^^-e(")(k II z)h^^\k,T) , 



A=±2 



Vb(k II z,r) ■ n = icosek,^vi'\k \\ z,r) + ^ ^ sin0,,„e^*^'=."^;(')(k || z,r) 

zl V2 



A=±l 



Z COS 9k,ne'\^ II i)vi'\k,T) 

+ J2 -«sin4,„e^*<^^--"e^^Hk || z)vl^\k,T) 



A=±l 



(3.58) 
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where and are the initial stochastic variables of scalar, vector and tensor modes, 

and we use the calculation results from Appendix [Dl as 



or(z) 




(3.59) 



According to Refs. 



14 



18 



19 



if we equate with the comoving curvature perturbation on 
superhorizon scales TZ, the initial conditions of the metric perturbations and the baryon velocity 
are El 



^iso^(^' ''"ini) 



2 

mi J 1 



2(15 + 4i?^ 



-{kn 



2 

mi J 1 



(3.61) 



vl'\k,n^O = ^{knnO' ■ 



In the tensor mode, equating with the primordial gravitational wave on superhorizon scales 
h^'^'^K it is satisfied that 



h^^\k,n. 



2V2 



-(/cTii 



(3.62) 



2(15 + 4/?^ 

Note that the Doppler effect does not affect only the tensor-mode perturbation. We introduce the 
transfer function in the Fourier space as |4| [20] 

z,n)=e(°)(k II z)Af)(r,A;,/ifc,„) , 
e(°)(k II i)A^p\T,k, 



Af(r,k 



(A^^^±.A^^O(r,k 



Ar(r,k 



(A(r)±.A[r^)(r,k 



Af (r,k 
(Ag")±.A[J))(r,k 



z,nj 
z,n) 

z,n) 
z,n) 

z, n) 



5 l^k,n I ) 



E -N^'l^lne^"''-^^"^^ II z)Ar)(r,fc, /.,,„) 



A=±l 



TK^T\^^Kn)^Jl-^4,n^^'^'^"i^^\^ II z)AirHr,fc,/i,,„) (3.63) 

A=±l 

E e^^^'^-e^^Hk II z)Af(r,fc, /.,,„) , 

A=±2 

E ^/i^) e^^<^'=-e^^)(k II z)A(.^Hr,A;,/ifc,„) , 



A=±2 



Then, the parameters in Ref. 14 are given by 



f?(fc,Ti„i) = -Tl 



5 + 4i?, 2' 
12(15 + 4i?,) ^ ™^ , 



(3.60) 



^(fc,Ti„i) = --7^(fcTi„i)^ 
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where /ifc,n = k - fi. Then, from Eqs. (3.52) and (3.53), we can write the collision term of Thomson 
scattering for the scalar mode: 



/ <ifi'Ae<')(n',n) = / H / .ifi'E l^«'°'(k II i) 

m 

X '[YU^)Y*^{u') + lQYUh)Y;^{n')]Kf\T,k,^') 



.s=±2 



(2vr)= 



X 



F4(nOAf (r,fc,/i') 
-Vf |E^^2:n(n')|Aif)(r,A;,/i') 



.s=±2 



(3.64) 



for the vector mode: 



j d(]'Ae(^)(n',n) 



(2vr; 



3 / - ^ ^ 5 V 3 

m A=±l 



^^'EEn/?^^^'Hk II z) 



X ^{r2n^(^i)r;^(n') + 10Fo™(^lWr„(< 

xru(n')Af)(r,A;,/iO 



.s=±2 



/cin'(AQ(^)±.Af/(^))(n',n) = -/^ / dfi' ^^.F^A) E ^^^'^k || z 



zF4(nOriA(nOA;^)(r,A;,/i') 
3 

s=±2 



'(3.65) 
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for the tensor mode: 



771 A — i2 



X 



{Y,^{n)Y*Jn') + 10Yom{n)Y*^{n')} 



.s=±2 



X 



>^2*™(n')>^2A(n')Af(r,A;,/.0 



,s=±2 



Here, we use 



Stt, 



Using the multipole expansion as 
and the Q'-integrals for A = 0: 



IGtt 



±2l^2A (for A = ±1) , 



E(-^) V47r(2Z + l)llo(nOAj%7^)(r, A;) 



/ 



^2m^0 = Si^2Sm,Q , 

J dQ'Y^^Yio = 5ifl5mfl , 
dn'±2Y*^Yio = 



for A = ±1: 



dn'Y*Y^Yio = 



207r 



5m,A 5li 



H,3 



dn'Y*J^2xYio = , 



J dn' 



±2^m±2^A^0 



(3.66) 



2 ^ (3-67) 

1 ^ e'^^ = ^J\±2Y2x (for A = ±2) . 



(3.68) 



(3.69) 



(3.70) 
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and for A = ±2: 



dfl'Y*^Y2xYio = 



2^5 



(3.71) 



Sm,x SiQ =F — A(5i 1 + ^^5z,2 =F -^^S,3 + —5, 



28 
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lA 



we can obtain the anisotropies generated via Thomson scattering for the scalar mode: 



(2^ 

d^k 

(27r)3 



e(°)(k II z) 
e(°)(k II z) 




(3.72) 



2 ' 



for the vector mode: 

dn'AQ^^\n',h 



rf^]'(Ag(^) ±iA[/(^))(n',n 

^(^)(fc,r 
and for the tensor mode: 



dn'Ae^^\n',n 
dQ'iAQ'-^^ ±iAU^^^){n',n 



(2^ 



.A=±l 




E ^\/T^^-(n)^^'Hk II z) 



(2vr) 



.A=±l 



167r 



"±2 



F2A(A)e^'Hk II z) 



10 ^'^ 10 ^'^ 5 7 35 ^'^ ' 



3 T 



3 r 



6 T 



(27r)3 



E \/W>^-(n)e('Hk II z) 



(2vr) 



1 r 
10 



.A=±2 

E 

A=±2 

7 ■ 70 



^(^)(A;,r) , 




±2l2A(n)e^"Hk II z 



^(T) l^(T) J5_^(T) _ 3^(T) 6^(T) 



V;(^)(^,t) ,(3-74) 
A( 



3 T(T) 
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Thus, from the Boltzmann equation (3.47), the gravitational redshift and the Doppler terms (3.58), 
and the collision term of Thomson scattering (3.53) and (3.72) - (3.74), we derive the Boltzmann 
equation for the scalar mode: 



Af + lklik,nA\ 



iS) 



1 . 



his + ( f^ln - ^ ) h''^ 



A^P + tkfik,nA'; 



AiS) 



K 



-K 



2 



iS) 



iS) 



(3.75) 



Af-\{l-^.l^)i.^^^ 
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for the vector mode: 



A^^ + ik^k,n^p 



iV) 



-K 



A 



(V) 



(V) 



(3.76) 



and for the tensor mode: 



(T) 



(3.77) 



The multipole expansion of these equations gives 



A 



(S/V/T) 
I/P,l 



+ 



k 



21 + 1 



(/ + i)A 



(S/V/T) _ jUS/V/T) 
I/P,l '■^I/P,l-1 



.AS/V/T) 
I/P,l 



(3.78) 



where 



y{S) 



V, 



(S) 



PA 



I 



y{V) 
^P,l 

yiT) 
^1,1 



(T) 



P,l 



-K 



1 . 



K 



-K 



A 



(T) 
1,1 



— K 



(3.79) 



Following the line of sight integration 21 , we can give the explicit solution of this Boltzmann 



equation as follows. At first, let us present the derivation for tensor mode. Using de /dr = ke 
and multiplying e^''t^k,nr-K both sides of Eq. (3.77), we have 



d 



A 



dr 



(T) 



(3.80) 



where g{T) = ke~'^ is the visibility function which describes the probability that a given CMB pho- 
ton last scattered at a given time. Through the integral over conformal time and some treatments, 
we obtain each from of the transfer function at r = Tq as 



Ay 



(T) 



TO 



dre 



(T) 



(3.81) 
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Here we use k{to) = 0, k{t = 0) ^ oo and x = k{To — t). In the same manner, we also obtain the 
scalar-mode function: 



Ay 



-L 



+9\ — + AV 



(5) , 33^(^ 1 



■AS) - 



+9 



h(S) 3^{S) ':^g^{S) 



3 



+ 



4 k^ 



iS) 



(3.82) 



and the vector-mode one: 



^0 
^0 

A^P = - dTe-'^'^-"''gip^^^ 
Jo 



-e"" + g 



k 



+ vl ' +g 



(3.83) 



where in the second equality of each equation, we neglect the topological terms. 

Consequently, we can summarize the transfer functions when k || z for the scalar mode: 



Af (To,k|| z,n) 
(A(f)±iA[f))(ro,k|| z,n) 



e(°)(k II z) r dTe-'i"^'-''Sf\k,T) , 
Jo 

±2l2o(A)e^°nk II z) r dre-'^^^"^S^p'\k,r) , 
Jo 



4 /67r 
3 V T 



+ ^7 



P 4 F 4 A; 



(5)' 



3 ^V'^^) 
4~P~ 



(3.84) 



Sl^'\k,T)^^gij^'^{k,T) , 



(5), 
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for the vector mode: 



A;'')(ro,k|| z,n) 



and for the tensor mode: 

(ro,k|| z,n) 

(A[f)±^A[,^))(ro,k|| z,n) 



.A=±l 
.A=±l 




.A./^^F,.(n)e(^)(k II z) 



TO 



dTe-'f''-"''Sf'\k,T) 



16n 



-±2 



Y2x{nK^'\k II z) 



dTe-"''^--''S^p \k,T) , 
(3.85) 



A; 




y2i(n)?'^>(k II z 



E 

.A=±2 



W^,.F..(n)e(^)(k II z) 

LA=±2 



TO 



dTe-'>"''"''SP{k,T 



dTe-'^'>^'-''Sf\k,T) , (3.86) 



Here, we have introduced the source function, s\^i^p^'^\k^T). 

3.7 All-sky formulae for the CMB scalar-, vector- and tensor-mode 
anisotropies 

In this subsection, let us formulate the all-mode CMB coefficients aim in the all-sky analysis on 
the basis of the derivation in Refs. [sj [T]. 

Since the CMB anisotropy is described in the spherical coordinate system, its intensity (/) and 
two polarization [Q and U) fields should be expanded by the spin-0 and spin-2 spherical harmonics, 
respectively, as 



e(^)(n) = 5^afiF,„(n) 



£,m 



(g(^) ±^^(^))(n) = J2a%^ ^,YUn) • 



(3.87) 



i,m 



Here, the index Z denotes the mode of perturbations: Z = S (scalar), = V (vector) or = T (tensor). 
The main difficulty when computing the spectrum of polarization arises from the variance under 
rotations in the plane perpendicular to fi. While Q and U are easily calculated in a coordinate 
system where k || z, the superposition of the different modes is complicated by the behavior of Q 
and U under rotations. However, using the spin raising and lowering operators d , d defined in 
Appendix [A], we can obtain spin-0 quantities. This leads to the rotational invariant fields like the 
intensity one and there are no ambiguities connected with the rotation of coordinate system arise. 
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z ^ n : line-of sight direction 

(photon's traveling direction) 




k : wave number vector 

► y 



A' 



^ y 



Figure 3.4: Geometry for the line-of-sight direction. 



Acting these operators on Q ±iU in Eq. (3.87), we have 



^2(g(^) + ^[/(^))(n) = ^ 

Em 



im 




(3. 



(z) 

Instead of a^2^„, it is convenient to introduce their linear combinations as 



22 



''Elm 
'■Blm 



_1 f _L (^) 



"2,^m "■-21m 



(3.89) 



These fields ii^ and B have parity-even and odd properties, respectively, in analogy with the electric 
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and magnetic fields. Then, from Eqs. (3.88) and (3.89), we can express 



Aif^(ro,k,n) 
Aif)(ro,k,n) 



1 


r(£-2)!i 


1/2 


2 


.(£ + 2)!_ 





i 


r(£-2)!i 


1/2 


2 


.(£ + 2)!_ 





(27r)3 

d"2(A(f)+2A[f))-d2(A 



2.a(^) _,A(f)^ 



(ro,k,n), (3.90) 
(ro,k,n) . 



Here, X discriminates between intensity and two polarization (electric and magnetic) modes. 



respectively, as X = I,E,B. When k || z, using Eqs. (3.84) - (3.86) and the operations derived 



from Eq. (A. 3) as 



^2 [2F2A(n)e-'^'=-"] 
^ ' [-2l^2A(n)e-'^'=-"] 



-d 



A 



1 - t^ln 

A 

1 ~ f^k,r, 



[(l-/x|j2r2A(l^M)e-^'^'='"1 , 
[(l-/.^.^J_2l^2A(f^M)e'^"'=-1 



we obtain more explicit expressions as 



ASfH^o,k|| z,n) 
Ag^H^ckll z,n) 



^^-2)'l''^(0)(k||z 




rfr5jf^(A;,r)^(^)(x)e-^'^'=-' 
^ -Az./^yu(A)e(^Hk II z 



.A=±l 




'o 



Ag'Hro,kl| z,n) 



r(£-2)!i 


1/2 


.(£ + 2)!_ 





E 

.A=±l 




-ZA/^ru(n)e(^)(k II z) 



X 



/"f"o 

/ c/r4^^(fc,r)i3(^)(x)e-'^'=-" , 
Jo 



Ag^)(ro,k|| z,n) 



r(£-2)!i 


1/2 


.(£ + 2)!_ 





.A=±2 



E \/^>^2A(fi)e(^nk II z) 



X 



Ag^)(ro,k|| z,n) 



r(£-2)!i 


1/2 


.(£ + 2)!_ 





.A=±2 



E-^\/^>^2A(fl)e(^)(kl| Z) 



2 V 15 



X 



rfr4''HA;,r)i3(^)(x)e- 



(3.91) 



(3.92) 
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with the operators B defined as 

^(^)(x) = (l + 9|)V, 

£^^\x) =4x+ (12 + x^)d^ + 8xdl + x^dl , 

B^^\x)=x^ + Axd^, + x^dl , 

£^^\x) = -12 + x^il - dl) - 8xd^ , 

B^^\x) = 8x + 2x^d^ . 



From here, we want to show analytical expressions of a^m's- In the above discussion, we have 
analytical formulae of the transfer functions when k || z. This implies that we consider the physics 



in the blue basis of Fig. 3^ and their transfer functions are completely determined by not the angle 
between z and fi, namely but the angle between k and n, namely Qk,n- However, as shown in 



Eq. (3.90), to obtain a^m's, we have to consider the physics in the red basis Fig. 3.4 and perform 



the f2„-integral. Here, instead of the transformation of the transfer functions, it is a better way to 
transform the integration variable in the agrn as f2„ — )■ ^lk,n- This can be done by using the Wigner 
D-matrix -D^|„/, which is the unitary irreducible matrix of rank 2i+l that forms a representation 
of the rotational group. The property of this matrix and the relation with spin-weighted spherical 
harmonics are explained in Appendix [Bj If we consider the rotational matrix 



cos 6k cos (pk — sin (p^ sin 6^ cos (pk 
cos 9k sin (pk cos (pk sin 9k sin (pk 
— sin 9k cos 9k 



(3.94) 



corresponding to the configuration (a = (pk, l3 = 9k,^ = 0) of Eq. (B.3) and satisfying 



the transformation equation (B.l) can be equated with 

= isink))YUnk,. 



(3.95) 



(3.96) 



Using this equation and the relation of the coordinate transformation as dQn = dQk,n, the agrn of 
arbitrary mode is written as 



''Xlm 



(27r)3 

4:71 

2i+l 



E^-!n' iS{Qk)) / dQk,nYU^k,nWx\ro,k\\ z, A) 

m' 

f d^k " 



(27r) 



Y.i-'^r' -m'Yl^i^k) J dQk,nY;^:{^k,nWx\ro,k\\ z,n) 



In the second equality, we have obeyed the relation of Eq. (B.4) in this 

1/2 



An 
2e+l 



i-ir -m'Yl^m ■ 



(3.97) 



(3.98) 
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Using the mathematical results of the fifc,n-integrals 

5; 4. 



2L + 1 



^Lo(^fc,n) 



/ 



^O, Y* p-^l^k,nX 



(3.99) 



(£-2)! a;2 

we can find the general formulae of the a^m for all-mode perturbations: 



''Xlm 



f2^ 



^[sgn(A)]^+-_,r;„(fi,)e(^)(k)ri5(fc) 



(3.100) 



where the helicity of the perturbations is expressed by A: A = for {Z = S), = ±1 for {Z = V) or 
= ±2 for {Z = T), the index x discriminates between the two parity states: x = for X = /, i?, 
X = 1 for X = B, and the time-integrated transfer functions Tx){k) are expressed as|^ 

fTO 



rZ\k) 
rSik) 



drS} '{k,T)je{x 
(£-2)!]^/^ P 



(£ + 2)! 

(£+1)!]^/^ P 
+ 1 



dTS\^^S^^\x)je{x) 



dTS^P{k,T 



^ + 2) 



e-2) 



ii-2)\ 



1/2 



{£ + 2)\ 
(£-2)!^'/' 



TO 



^ Ji{X) 
X 



drSp {k, t)£^' '[x 



X 



(£ + 2)! 

1/2 



"dr4^)(/t,r)i3W(-^^^'^^ 



(3.101) 



X) 



X 



drSPik^T 



JAx) 



x^ 



dTSP{k,T)S^''\ 



x] 



,Jeix) 



x^ 



dTSP{k,T)B^'^\:> 



Note that in the all-sky analysis, due to the dependence of transfer functions on 0fc,n, o^m's depend 
on the helicity state through the spin spherical harmonics. In the above discussion, we take the 
synchronous gauge and derive the aim- However, in the same manner, we can obtain the identical 
form of Eq. (3.100) even in another gauge. In this case, the different points can be confined only 



^In Ref. [s], there are three typos: right-hand sides of Eqs. (B21), (B22) and (B23) must be muhiphed by a 
factor —1, respectively. 
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in the transfer function (3.101). In a numerical code CAMB 13, 19 , these transfer functions of 



polarization modes are expanded as, e.g., 



To 



dr 



g,p(^) + 2#(^) + gi,(^) ^ 4 #(^)+#(^ _ ^^^^^ / _ _6 

k X V 



rZ\k) = -2l dr 

'0 



TO 



A; 



(3.102) 



3.8 Flat-sky formulae for the CMB scalar- 
anisotropies 



vector- and tensor-mode 



The fiat-sky approximation uses the (2D) plane wave expansion of the CMB fluctuation instead 
of the spherical harmonics one, and it is valid if we restrict observed direction (parallel to fi) only 
close to the z axis |4 , 23 . As confirmed in Ref. |4j, the flat-sky power spectra of E- and i?-mode 
polarizations sourced from the primordial scalar and tensor perturbations are in good agreement 
with the all-sky ones for i > 40. On the basis of these studies, we have also compared the all-sky 
power spectra with the flat-sky ones for the J, E, and B modes from scalar, vector, and tensor 
perturbations and found their consistencies at £ > 40. 

As mentioned in the previous subsection, in order to estimate the aim, one must construct 
the transfer functions for arbitrary k. In other words, we want to obtain the transfer functions 
expressed by arbitrary k (whose direction is denoted by flk) and n (denoted by fi„) instead of 
flk,n in Fig- 3.4, In the / modes, only by changing Qk,n to Qk and Qn with the relation (3.95), 



the transfer functions for arbitrary k can be obtained. In the E and B modes, in addition to this 
treatment, one must consider the mixing between Ag and Au under the transformation S{Qk) as 
described in Ref. |j4j. This effect is expressed as 



(A 



[zy 



±^A[f)')(ro, 



n 



(3.103) 



with the mixing angle ip. The angle represents the rotation angle between 6k,n and On, where 
Ois^n and 9n are the unit vectors orthogonal to fi in a particular basis in which k || z and a general 



basis, respectively (see Fig. 3.4) 



In the flat-sky analysis, i.e.. On 0, by using Eqs. (3.84) - (3.86) and by using the limit of ijj 
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as ■?/'—>■ 0ri — 0A; + TT, the transfer functions for arbitrary k are derived as 



Af)(ro,k,n)^e(k) / dTS'/'{k,T)e 



TO 



-iknD 



(A(f)±zA[;0(ro,k,n)^e 



T2i{ 



TO 



sin'^fee(k) / dTS'p\k,T)e 



-ikiiD 



Af)(ro,k,n) ^zsin^, ^'^k)] / dTSi''\k,T)e 
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X 



X 
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(3.104) 



dTS'^p\k,T)e- 



ik-nD 



Af (ro,k,n)^sin2 

(Ag^^ ± iA^P){to, k, n) ^ e^2*('^"-^'=) 
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(l + cos2^^fc)(5^e^^)(k)) 

\A=±2 / 

t2 cose ^e^'^k) 



ik-nD 
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drS^lk, T)e 



-ikiiD 



It is important to note that the (pk dependence which are inherent in the vector and tensor per- 
turbations vanishes in the fiat-sky approximation, besides a trivial (pk dependence due to a spin-2 
nature of the Stokes Q and U parameters. One may explicitly see that 0fc,n dependence vanishes 
in the transfer functions when taking 6n — )■ because the S{flk) matrix rotates the basis with the 
new z axis always being on the x — z plane in a particular basis in which k\\ z (see Fig. 3.4). This 
approximation means that for 0„ ^ 1, it is valid to calculate the CMB fluctuation on the basis of 
vector and tensor perturbations fixed as = 0, namely, (j)k,n = t^- 

In the fiat-sky limit, aim in the all-sky analysis described as Eq. (3.90) is modified by using the 
plane wave as 



d^k 



A(f)(ro,k,n)^aV"^(^) 



Aif)(ro,k,n)^ 



(27r)3 



s=±2 



(3.105) 



A(f)(ro,k,n)^-^ 



s=±2 



(ro,k,n)e" 



-Si{(f)i-<f>n) 



where © is the 2D vector projecting fi to the fiat-sky plane expressed as = (Ocos0„, 9sin0„) Fl 
For example, in order to obtain a}/^, we substitute Eq. (3.104) into Eq. (3.105) and calculate 



^Not confuse with the CMB temperature fluctuation. 
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as follows: 



-! 



TO /*CXD 



dT 



X 



(3.106) 



where D = tq — t is the conformal distance and we have decomposed k into two-dimensional vector 
parallel to the flat sky and that orthogonal to it, k = (k", kz)- In order to obtain the last equation, 
we use following relations which are satisfied under kH = —l/D as 



k = \ lk^ 
sin 9k 



D 



kD ^{kzDf+e ' 



(3.107) 



cos6'fe = sgn(/c^)^/l 



kD 



The other-mode aim's can be calculated in the same manner. Thus, we summarize the all-mode 
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a{£ys: 
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S^P{k = ^kl + {i/D)\T) , 
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(3.108) 
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Sf\k=^M+WW,r) 



S^Pik=^/k[T(£/W,r) 



2i sgn(/c^)Wl - 



£2 



sP{k = PkfTWW,r) 



{k.Df + & 
d source fun 

spectrum and bispectrum by using these formulae |3| |24}]26|. 



(z) 

where we label S)^^ as the modified source function. One can formulate the flat-sky CMB power 



3.9 CMB power spectrum 

To extract several information about the Universe from the observational data, the two-point 
correlation function of the CMB fluctuations (called CMB power spectrum) is often-used. Here, 
we formulate the CMB power spectrum and summarize the constraints on several model parameters 
from the current observational data. 

If we assume the Gaussianity and the symmetry under the parity and rotational transformations 
in the initial stochastic variables, their power spectrum can be expressed as 



/^(A.)(k^)^(A.)*(k2)) = (27r)3p^(A:i)5(ki - k2)5A„ 



1 



A2 



iZ = S) 



(3.109) 



1/2 {Z = V,T) 

This implies that the couplings between the different modes of the perturbation vanish in the 



power spectrum. Then, from the formula of the all-sky aim (3.100), the CMB power spectra of all 
modes are derived as 



\n=l 



(3.110) 
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where we use the relation derived by the reality condition of the metric perturbation: 



(3.111) 



In Fig. 3.5, we plot the CMB intensity and polarization power spectra of the scalar and tensor 



modes. Here we think that the scalar and tensor perturbations are sourced from the primordial 
curvature perturbations (^'^"^ = TZ) and primordial gravitational waves (^^^^■' = /i*-^^-*), respectively. 
The ratio between these power spectra, called the tensor-to-scalar ratio r, is defined as |^ 



r = 



2PT(fc) 
Psik) 



(3.112) 



At first, let us focus on the // spectra in the left top panel. In the scalar spectrum, the dominant 
signal is generated from the acoustic oscillation of the photon and baryon fluid. The first acoustic 
peak is located at £i ~ 220. This scale is corresponding to the angle of the sound horizon at 
the recombination epoch as £i ~ 27rro/rs(2;*). At small scales corresponding to £ 3> £i, due to 
the difference between the photon and baryon speeds, the coupling between photons and baryons 
is ineffective and the acoustic oscillation is highly damped. This effect is well-known as the Silk 
damping 28 . On the other hand, the gravitational blue shift due to the potential decay in the late 
time affects the fluctuations for i <^ ii. This is called the integrate Sachs Wolfe (ISW ) effe ct [29 ] 



(z) 

and arises from the terms proportional to e"*^ of the source function Sf' in Eqs. (|3.84[) and (|3.86p. 



In the tensor spectrum, the ISW effect leads to the d ominant signals and the scattering effec t (the 
second term of the source function Sf^^ in Eq. (3.86)) is subdominant [g]. As shown in Eqs. ( |3.84 ) 
and (3.86), the EE and BB spectra (the bottom two panels) have no gravitational redshift term 
and there is no ISW effect. Instead of it, the scattering term generates the CMB fluctuation. 
The most interesting signature in the polarization power spectra is the enhancement at £ < 10. 
This can be caused by Thomson scattering at reionization of hydrogen which may have occurred 
at z ~ 10 Therefore, these signals can be important to determine the optical depth of the 
Universe, k. The IE spectrum described in the right top panel seems to include the features of 
both the // and EE spectra and has both positive and negative values. Unlike the above three 
cases, the BB spectrum never arises from the scalar mode because the scalar mode has only one 
helicity, namely A = Hence, we believe that the BB spectrum directly tells us the amplitude 
of the primordial gravitational waves depending on the energy scale of inflation. 

Theoretically, the CMB power spectrum depends on the parameters which determine the 
dynamics of the Universe as the energy density of the cosmic fluids, curvature, and the Hub- 



(S) 

ble constant Hq. Figure 3.6 shows the dependence of Cjj ^ on the density parameters of cold 
dark matters, the cosmological constant and baryons as lUc = J]c(i^o/100sec-Mpc/km)2,fiA and 

Ub = r2fe(iJo/100sec ■ Mpc/km)^, respectively. From this figure, one can observe that as Uc de- 

(s) 

creases, the overall amplitude of C)^ ^ enlarges. This behavior is understood as follows. If ujc 
decreases, since the radiation dominated era is lengthened, the gravitational potential for smaller 
k enters the horizon and decays. Thus, Cjj^ at corresponding multipoles as £ ~ fcr is boosted due 
to the gravitational blue shift. This is the so called early ISW effect l29j. Next, focusing on the 
blue dotted line, one can find that if becomes large, C)^ ^ is boosted for £ < 10. This is due to 
the late ISW effect, that is, A dominates the Universe earlier and the potential at larger scales is 



27 



and CAMB 



13 



^This definition is consistent with Eq. (2.69), and the notation of Ref. 

^The vector mode generates the BB spectrum due to its two heiicities. However, due to the decay of the vector 
potential via the Einstein equation, this becomes the subdominant signal. To avoid this, the sources such as cosmic 
strings and magnetic fields need to exist and support. 



3.9 CMB power spectrum 



59 



10000 
1000 
100 
10 

1 

0.1 
0.01 
0.001 











scalar 




tensor 





















\ 

\ 


tensor - 










10 


100 


1000 



Figure 3.5: The CMB spectra of the // (left top panel), IE (right top one), EE (left bottom one) 
and BB (right bottom one) modes. Here we consider a power-law flat ACDM model and fix the 
tensor-to-scalar ratio as r = 0.1. The other cosmological parameters are fixed to the mean values 



reported in Ref. 27 



(S) 

destroyed, hence C)j ^ at corresponding €s is amplified. We also notice that when uib enlarges, the 
ratio of the amplitude between the first and second peaks of the magenta dot-dashed line increases. 
Solving the coupled Boltzmann equations, the acoustic oscillation of the baryon-photon fluid in 
the matter dominated era is roughly given by 

e ~ Q + i?^ $ cos(A;r,) - , (3.113) 

where $ is the potential of the conformal Newtonian gauge, denotes the sound horizon and 
R = 3pfe/(4p^). Then if cob increases and R becomes large, this equation experiences increase in 
amplitude and suppression of the intercept. Hence, the difference of Cj^\ oc 9^ between the odd- 
and even-number peaks increases. These parameters are limited with the others from the current 



observational data set as Table 3.2, Other than these parameters, massive neutrinos and some 



relativistic components also make impacts on the CMB fluctuation (e.g. Refs. 30, 31]). 

The CMB power spectrum also depends on the primordial curvature perturbations and pri- 
mordial gravitational waves. Conventionally, these spectra are parametrized as 

- U.002MPC-1 J > (3.114) 
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Figure 3.6: CMB // power spectra sourced from the scalar-mode perturbations in a power-law 
flat ACDM model. The red solid line is plotted with Uc = 0.112, = 0.728,^6 = 0.02249, = 
0.967, K = 0.088 27 . The green dashed, blue dotted and magenta dot-dashed hues are calculated 



if Uc decreases to 0.1, increases to 0.8, and Ub increases to 0.028, respectively. 



and Eq. (3.112). As shown in Eq. (3.100), the magnitudes of the primordial curvature perturba- 

( s) 

tions As and gravitational waves rAs simply determine the overall amplitude of the Cx/x2 ^^"^ 

(T) 

^XiX2 h ' respectively. The spectral index of the power spectrum of the curvature perturbations 

changes in slope of C*XiX2^i- F'^o'^ the observational data, the constraints on As^Ug^r are given 
as Table |3.2[ Here, we want to note that ng = 1 is excluded at about 3-sigma level. This implies 



the deviation from the exact de Sitter expansion in inflation. As shown in the bound on r, unlike 
the primordial curvature perturbation, the primordial gravitational wave has not been detected 
yet. However, some experimental groups aim to discover the BB spectrum through remove of 
some noisy foreground emission and improvement of the instruments [2| [34}]36] . If these projects 
achieve, it will be possible to judge the existence of the primordial gravitational waves of r < 0.01. 

So far, we discussed under the assumption that the parity and rotational invariances are kept. 
On the other hand, there are a lot of studies which probe the somewhat exotic scenarios where 
these invariances violate. As a theoretical prediction, if parity-violating action such as the Chern- 
Simon term exists in the early Universe, (^(+2)(k)^(+2)*(i^/)^ _^ and the IB 

and EB spectra appear 37-40 . Using the parametrization as 



C,°B:^ = C,£;,,sin(2Aa) , (3.115) 

the parity violation is limited as —5.0° < Aa < 2.8° (95% CL) [27]. The rotational invariance is 
broken if the Universe has the preferred direction and this situation is realized by the anisotropic 
inflation 41-43 . This leads to the direction-depending power spectrum as 

Psiy) = pr{m+9{^-^f] , (3.116) 

and produces the off-diagonal components in the CMB power spectrum as li 7^ ^2- From the CMB 
observational data, the magnitude of the statistical anisotropy has been limited as (7 = 0.15 ±0.039 
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Parameter 



WMAP 7-yr 



WMAP + BAG + Hn 



lOOwfe 

UJc 

Us 
K 

As 

r 



2.249] 



hO.U56 
-0.057 

0.1120 ±0.0056 



n 797+0.030 
'^'-0.029 

0.967 ±0.014 
0.088 ±0.015 
(2.43 ±0.11) X 10-^ 
< 0.36 



2.255 ±0.054 
0.1126 ±0.0036 
0.725 ±0.016 
0.968 ±0.012 
0.088 ±0.014 
(2.430 ±0.091) X 10-^ 
< 0.24 



Ho 


70.4 ± 2.5 km/s/Mpc 


70.2 ± 1.4 km/s/Mpc 




0.0455 ± 0.0028 


0.0458 ±0.0016 




0.228 ±0.027 


0.229 ±0.015 




n 1 Q4c;+0.0056 
U.±0^d_Q 0055 


0.1352 ±0.0036 


2^reion 


10.6 ± 1.2 


10.6 ± 1.2 


to 


13.77 ±0.13 Gyr 


13.76 ±0.11 Gyr 



Table 3.2: Summary of the cosmological parameters of AGDM with finite r model from the WMAP 



7-year data 27 , and the data set combined with the results of the galaxy survey 32 and Hubble 

respectively. Here z^eion denotes the redshift at the reionization epoch, 

I Ub + UJc- 



constant measurement [33 
to is the present time of the Universe, and uj. 



Furthermore, owing to the progress of the observational accuracy, the deviation of the Gaus- 
sianity can be measured. Beyond the power spectrum, this is achieved by using the three-point 
function (bispectrum). In the next section, we discuss about how to extract the information on 
the early Universe from the GMB scalar, vector and tensor bispectrum. 

In addition, we can add other components of fluids in the analysis of the GMB spectrum. From 
Sec. [9| we focus on the effect of the primordial magnetic fields on the GMB. 
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4 Primordial non-Gaussianities 

The study of non-Gaussian impacts in the cosmological fluctuations provides an important infor- 
mation of the early Universe [l]. The primordial non-Gaussianities are measures of the interactions 
in inflation, hence constraining on this will lead to a great deal about the inflationary dynamics. 
It may also puts strong constraints on alternatives to the inflationary paradigm (e.g., Refs. [2]-[8]). 

In the previous part, we expanded the inflationary action to second order in the comoving 
curvature perturbation TZ and the gravitational waves These actions allowed us to compute 

the power spectra P-R,(fc) and Ph{k). If the fluctuations TZ and /i^^^^ obey the exact Gaussian 
statistics, the power spectrum (or two-point correlation function) contains all the information 
However, for the non-Gaussian fluctuations, higher-order correlation functions beyond the two- 
point function contain additional information about inflation. Estimating the leading non-Gaussian 
effects requires the expansion of the action to third order since we must take into account the leading 
non-trivial interaction terms. In this section, we review recent studies about the primordial non- 
Gaussianity based on e.g., Refs. [9 11 . 



4.1 Bispectrum of the primordial fluctuations 

At flrst, we give the deflnition of the bispectrum of the initial perturbations ^^^^ of the scalar 

(A = 0), vector (A = ±1), and tensor (A = ±2) modes. The Fourier transformation of the 
two-point function is the power spectrum 

rie(^"nk„)|) = i2nrPzih)S (j^K^ 5.„..(-l)^^ X = ^1^^ . (4.1) 

Similarly, the Fourier mode of the three-point function is so called the bispectrum 

ne^^")(k„)\ = (2vr)3F^^^^^3(k,,k2,k3)5(X^kJ . (4.2) 



\n=l / \n=l 



Here, the delta function enforcing the momentum conservation is a consequence of the translation 
invariance of the background. If the scale invariance is kept, we have 

F^^^'^-^{bk^,bk2,bks) = 6-^F^^^^^3(ki,k2,k3) . (4.3) 

Moreover, the rotational invariance reduces the number of the independent variables to just two, 
namely k2/ki and k^/ki. 

In order to compute the primordial bispectrum for a speciflc inflation model, we need to treat 
the time evolution of the vacuum in the presence of the interactions carefully. On the other hand, 
for the power spectrum, this is not the leading order effect. As the methodology for calculation 



of the cosmological correlation function, "the in-in formalism" is well-known 12 16 . In practice, 
computing the bispectra can be algebraically very cumbersome. In Sec. [8} using this formalism, 
we actually discuss the computation for the bispectrum of the gravitational waves. 



^AU odd-point correlation functions vanish for the Gaussian fluctuations, while all even-point functions can be 
expressed in terms of the two-point function due to the Wick's theorem. In other words, all connected higher-point 
functions vanish for the Gaussian fluctuations. 
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Figure 4.1: Representations of triangles forming the bispectrum, F^^^'^^^{ki, k2, k^), with various 



combinations of the wave numbers satisfying k^ < k2 < ki. This figure is adopted from Ref. 19 . 



4.2 Shape of the non-Gaussianities 



The delta function in Eq. (4.2) enforces that the three wave number vectors form a closed triangle. 
Depending on the inflation models, the maximal signals arise from the several triangle configura- 
tions. Therefore, specifying the shape of the non-Gaussianity is a powerful probe of the mechanism 



that laid down the primordial fluctuations 17, 18 



We have a useful way to visualize the shape dependence of the bispectrum. We study the 
structure of the bispectrum by plotting the magnitude of F^^^'^^^{ki^k2ik^){k2/kif'{k^/kiY as 
a function of k2/ki and k^/ki for k^ < k2 < ki. To classify various shapes of the triangles, 
let us utilize the following names: squeezed {ki ?a /c2 ^ k^), elongated {ki = /c2 + ^3), folded 
(ki = 2/c2 = 2/C3), isosceles (^2 = ^3), and equilateral (fci = k2 = ks). In Fig. 4.1[ the visual 
representations of these triangles are presented. 

From here, we focus on the three representative shapes of the primordial bispectrum: "local" , 
"equilateral", and "orthogonal". Then, it may be convenient to decompose the non-Gaussianity 
of the curvature perturbations into the magnitude-depending and shape-depending parts: 



F°°°(fci,fc2,A;3) = -hL{'^7r'AsyS{k,,k2,ks) 
5 



(4.4) 



where As is the magnitude of curvature perturbations defined in Eq. (3.114). 



4.2 Shape of the non-Gaussianities 
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4.2.1 Local type 

The simplest way to parametrize the non-Gaussianity was via a non-linear correction to a Gaussian 
curvature perturbation TZg [20j 



local 



^(x) = 7^,(x) + - v 



(^.(x)^>] 



(4.5) 



From this equation, we can see that the non-Gaussianity is localized at a given point in the real 
space. Therefore, we call this the local-type non-Gaussianity, and /nl^^ is called "the local-type 
nonlinear parameter" . Then, the bispectrum of the local-type non-Gaussianity of the curvature 
perturbations is derived as 



^iocai(fci, h) = ^/i^r [Pnih)Pn{k2) + 2 perms.] 



6 



(4.6) 



where we fix the spectral index as ra^ = 1. Hence, equating with Eq. (|4.4|), we can write 

1 



2 perms. 



(4.7) 



Without loss of generality, we shall order the momenta such that < k2 < ki. Then, the 
bispectrum of the local non-Gaussianity dominates in the squeezed limit: fcs ^ fci ~ k2 as shown 



in the top- left panel of Fig. 4.2 In this limit, we obtain 

12 



pOOO 



i_i(fci, k2, ks^O) = -f^oc.^P^(k,)Pnik,) 



(4.8) 



How large is /nl*^^ generated from inflation? If the single-field slow-roll inflation occurs, /nl*^^ 
would be of order the slow-roll parameter, e ~ 10~^ |21 23 . More accurately, in Refs. 12 , 24 , the 



authors found that the coefficient of PnikijPn^k^) from the simplest single- field slow-roll infiation 
with the canonical kinetic term in the squeezed limit is (1 — Us). Therefore, one has a consistency 
relation between the scalar spectral index and the local-type nonlinear parameter in the single-field 
slow-roll infiation: 



local 
NL 



12 



[I 



(4.9) 



which gives /nl*^' = 0.015 for Us = 0.963. In contrast, the large local-type non-Gaussianity will be 



able to be realized in the models with multiple light fields during infiation |25j-|35|, the curvaton 
scenario 



36 38 , inhomogeneous reheating 39, 40 



4.2.2 Equilateral type 



The equilateral bispectrum is parametrized as 41 

1 



S'^^\k,,k2,ks)=6 



ikik2] 



+ 2 perms. > — 



(fci/cs^s) 



+ 



^1^2 ^3 



+ 5 perms. 



(4.10) 



This function approximately expresses the bispectra arising from the infiation models where the 
scalar fields have non-canonical kinetic terms. An example is the so-called Dirac-Born-Infeld (DBI) 



infiation 



42 



43 



equil 



, where /^l 



cx — 1/c^ for c<j <^ 1 with Cg being the effective sound speed at which 
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the scalar field fluctuations propagate. We can also find a lot of other models that can produce 
large /^l"'' 44-48 . The normalized equilateral spectrum S'"^^''^{ki,k2,k^){k2/kiY{k^/kiY , peaks 
at the equilateral limit, namely, ki = k2 = k^ as described in the bottom-left panel of Fig. 4.2 
The local- and equilateral-type bispectra are nearly orthogonal to each other, which denotes that 
both can be measured almost independently. 

In Sec. |8j we confirm that the graviton non-Gaussianity originated from the Weyl cubic action 
is also categorized as the equilateral type. 

4.2.3 Orthogonal type 

In addition to the local- and equilateral-type non-Gaussianities, we have the following parametriza- 
tion: 



S°'''^°^{k,,k2,ks) = Q 



{kik^k.y 



^1^2 ^3 



+ 5 perms. 



• (4.11) 



This is the so called orthogonal-type non-Gaussianity and constructed such that it is nearly orthog- 



onal to both the local-type and equilateral-type non-Gaussianities 49 . This shape approximately 



represents the forms appearing in a linear combination of the higher-derivative scalar-field inter- 
action terms, each of which yields forms similar to the equilateral shape. In Ref. I49j, the authors 



discovered that on the basis of the effective field theory of inflation approach 47 , a certain linear 



combination of similarly equilateral shapes can leads to a distinct shape which is orthogonal to 
both the local and equilateral forms. The orthogonal bispectrum S°'^^^°^{ki, k2, k^){k2/kiY{kz/ki)'^ 
has a positive peak at the equilateral configuration, and negative valley along the elongated con- 



figurations as seen in the top-right panel of Fig. 4.2 



4.3 Observational aspects 

Observational limits on the primordial non-Gaussianities are beginning to reach the interesting lev- 



els. Using the optimal estimators 11, 49 -52 , the latest constraints on the nonlinearity parameters 



from the CMB data (WMAP 7-yr data) are obtained as 53 



10 < f^l""' < 74 , -214 < /^r' < 266 , -410 < < 6 (95% CL) . 



orthog 



(4.12) 



As another approach for extracting the non-Gausssinity from the CMB data, the methods with 
the Minkowski functionals have been developed |54- 



56 



^(_^local 



The PLANCK satelite 57 and the proposed CMBPol mission 58 will give tighter bounds as 

Oil). 



'NL ; 5 and 2. At the level of /^l*^' ~ C'(l), we expect to see a signal from the secondary 
effects not associated with the non-linear effects in inflation. In order not to confuse these effects 
with the primordial signals, some researchers have been estimating in detail how the gravitational 

pi). 
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non-linear evolution of the fluctuations can induce its own non-Gaussianity (see, e.g. 

In addition, the primordial non-Gaussianity also leaves the signatures in the large scale structure 
in the Universe. In general, extracting the information of the primordial non-Gaussianity from the 
observational data of the large scale structure is complicated due to the fact that the gravitational 
nonlinear evolution of the density fluctuations produces the secondary non-Gaussianity completely 
dominating over the signals from the primordial origin. However, recently, the scale-dependent 
bias has been considered as a promising probe of the primordial non-Gaussianity 63 64 . In these 



4.4 Beyond the standard scalar-mode non-Gaussianities 
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Figure 4.2: Two dimensional color contour for the shapes of the primordial bispectra. Each panel 
describes the normalized amplitude of S{ki, k2, fc3)(fc2/^i)^(^3/^i)^ as a function of k2/ki and k^/ki 
for /cs < /c2 < fci- Since the primordial bispectra shown here are nearly scale invariant, the shapes 
may look similar regardless of the values of ki. The amplitude is normalized by a rule that it is 
unity at the point where S{ki,k2^k^){k2/kiY{kz/kiY takes on the maximum value. The top-left 
panel is the local form given by Eq. ( |4.7 ), which peaks at the squeezed configuration. The bottom- 
left panel is the equilateral form given by Eq. (4.10), which peaks at the equilateral configuration. 




The top-right panel is the orthogonal form given by Eq. (4.11), which has a positive peak at the 



equilateral configuration, and a negative valley along the elongated configurations. Note that all 
of these shapes are almost orthogonal to each other. This is quoted in Ref . pTj . 



papers, it has been confirmed that for highly biased tracers of the underlying density field, the 
bias parameter depends on scale and /nl, namely. 



Ps^^,{k) = [b + Ab{k, f^l^')Y Ps{k) . (4.13) 
Application of the method to the luminous red galaxies (LRGs) sample of SDSS leads to a bound 



as 
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- 29 < f^l^^ < 70 (95% CL) . 



(4.14) 



As one see, this limit is comparable to the bounds from the CMB data. The way to extract the 
information on the primordial non-Gaussianity from the matter distribution has continuously been 



studied (see, e.g., Refs. 65 67]). 



4.4 Beyond the standard scalar-mode non-Gaussianities 

Historically, as described above, only the non-Gaussianity of curvature perturbations has been 
well-known studied. However, the non-Gaussianity of vector- or tensor-mode perturbation can be 



generated from the cosmological defects 68 , 69 , the magnetic fields 70 , the nonlinear gravitational 
waves 12, 71-74, and so on. Furthermore, somewhat exotic non-Gaussianities including the 



violation of the rotatoional or parity invariance in the early Universe have recently been discussed 
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(see, e.g., Refs. 72, 75 77]). Hence, for precise comprehension of the early Universe, detailed 
analyses of these signals are necessary. 

This is the main motivation of our studies: construction of the general formulae for the CMB 
bispectrum with not only scalar- but also vector- and tensor-mode contributions, and computation 
and analysis of the CMB bispectrum sources from these novel non-Gaussianities. In the following 
sections, we focus on our studies. 
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5 General formalism for the CMB bispectrum from pri- 
mordial scalar, vector and tensor non-Gaussianities 

In this section, on the basis of the formulation of the CMB anisotropy in Sec. |3j we develop the 
formulae for the CMB bispectrum sourced from scalar-, vector-, and tensor-mode non-Gaussianity. 
These results have been published in our paper 

At first, we should remember an expression of CMB fluctuation discussed in Sec. |3} In the 
all-sky analysis, the CMB fluctuations of intensity, and E and 5-mode polarizations are expanded 
with the spin-0 spherical harmonics, respectively. Then, the coefficients of CMB fluctuations, 
called aim's, are described as 



— $^[sgn(A)]^+^'_,r;„(k)^W(k)ri: 



(5.1) 



where the index Z denotes the mode of perturbations: Z = S (scalar), = V (vector) or = T 
(tensor) and its helicity is expressed by A; A = for Z = S", = ±1 for Z = V or = ±2 for 
Z = T, X discriminates between intensity and two polarization (electric and magnetic) modes, 
respectively, as X = I,E,B and x is determined by it: a; = for X = I,E or = 1 for X = B, 
is the initial perturbation decomposed on each helicity state and T^f) is the time-integrated 



transfer function in each sector given by Eq. (3.101 



Next, we expand ^^^"^ with spin- (—A) spherical harmonics as 

e(^)(k)^5^eJiJ(fc)_,r,^(k) , 



(5.2) 



em 



and eliminate the angular dependence in Eq. ( |5.l| by performing k-integral: 
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(z) 

X,tm 



47r(-i)' 







(2vr) 



5^[sgn(A)]^nJ^)(A:)ri5(fc) . 



(5.3) 



Here, we use the orthogonality relation of the spin-A spherical harmonics described in Appendix 
lAlas 



(5.4) 



Then the CMB bispectrum generated from the primordial non-Gaussianity of the scalar, vector 
and tensor perturbations is written down as 



\n=l 



(Zn) 



\n=l 

^Here, we set = 1. 
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Figure 5.1: The CMB III (left top panel), HE (right top one), I EE (left bottom one), and EEE 
(right bottom one) bispectra induced from the local-type (red solid line), equilateral-type (green 
dashed one), and orthogonal-type (blue dotted one) non-Gaussianities of curvature perturbations. 
The three multipoles are fixed as ii = £2 = £3 = i. Here we consider a power-law flat ACDM 
model and fix the nonlinear parameters as f^l^^ = /nl"'' = /nl*^"^ 
parameters are fixed to the mean values reported in Ref. (Ol. 



100. The other cosmological 



This formalism will be applicable to diverse sources of the scalar, vector and tensor non-Gaussianities, 
such as, a scalar-graviton coupling (T] (Sec. [6]), cosmic strings (2|[3], primordial magnetic fields [4]-[6] 
(Sec. [9]), and statistically-anisotropic and parity- violating interactions (7| [s] (Sees. [T] and [s]). 

To compute the CMB bispectrum composed of arbitrary perturbation modes, we have to reduce 

the expanded primordial bispectrum, {^nn=i d^m„(^n))) involving the contractions of the wave 
number vector and polarization vector and tensor, and the integrals over the angles of the wave 
number vectors. As shown later, this is elegantly completed by utilizing the Wigner symbols and 
spin-weighted spherical harmonics. 

If the initial bispectrum satisfies the rotational invariance, the CMB bispectrum is divided 
into the Wigner-3j symbol depending on the azimuthal quantum numbers and the angle-averaged 
function as 

/rTfl^^"^ \ - ( ^2 r(^i^2^^) (^6) 

\ ll"^nm„ y ~ \ mi m2 ms J -^1^2X3 a^2<?3 ■ y^-^) 



Let us focus on the CMB bispectrum from the curvature perturbations. Then, from Eqs. (4.2), 
(4.4), (4.7), (4.10), (4.11), (5.5), (5.6) and the knowledge of Appendix [C| we derive the reduced 
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Figure 5.2: The CMB III (left top panel), I IE (right top one), I EE (left bottom one), and EEE 
(right bottom one) bispectra induced from the local-type (red solid line), equilateral-type (green 
dashed one), and orthogonal-type (blue dotted one) non-Gaussianities of curvature perturbations 
with /nl*^' = /nT'' = /nl^°^ = 100- Here, we fix the two multipoles as ii = £2 = 200, and plot 
each curve as the function in terms of £3. The cosmo logical parameters are identical to those in 
Fig. 10 



bispectra as 



i,(.SSS) 



n 

n=l 



where we have introduced the / symbol as 



(5.7) 



JS1S2S3 



(2Zi + l) (2/2 + 1) (2^3 + 1) 



47r 
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S3 



(5.8) 



In Figs. 5.1 and 5.2| we plot these CMB bispectra for each multipole configuration. From these, 
one can see that depending on the shape of the primordial non-Gaussianity, the overall magnitudes 
of the bispectra do not change, but the shapes in i space change. Comparing these with the CMB 



data, the bounds on the nonlinearity parameters (4.12) are obtained. 



In this simple CMB bispectrum, there exists no dependence of the initial bispectrum on the 
polarization vector and tensor, hence we can derive the above formulae easily. Considering the 
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vector- and tensor-mode contributions in the CMB bispectrum, however, the computation becomes 
so cumbersome due to the comphcated angular dependence of the polarization vector and tensor. 
This difficulty is also true for the bispectrum where the rotational or parity invariance is broken. 
From the next section, we deal with these complicated bispectra depending on the several scenarios 
by applying the wonderful mathematical tools such as the Wigner symbols and the spin-weighted 
spherical harmonics. 
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6 CMB bispectrum induced by the two scalars and a gravi- 
ton correlator 



In this section, adapting Eq. (5.5) to the primordial non-Gaussianity in two scalars and a graviton 



correlation JlJ , we compute the CMB scalar-scalar-tensor bispectrum. This discussion is based on 
our paper [2]? 



6.1 Two scalars and a graviton interaction during inflation 

We consider a general single-field infiation model with Einstein-Hilbelt action |3] : 



S 



d 



-9 



Ml, 



(6.1) 



where g is the determinant of the metric, R is the Ricci scalar, Mpj = l/(87rG), is a scalar field, 
and X = — g^'^ d ^(pdjjcf) / 2 . Using the background equations, the slow-roll parameter and the sound 
speed for perturbations are given by 



H 
IP 



Xp, 



X 



2Xp^xx + P,x 



(6.2) 



where H is the Hubble parameter, the dot means a derivative with respect to the physical time t 
and p^x denotes partial derivative of p with respect to X. We write a metric by ADM formalism 



ds^ = -N^dt^ + a^e^^^'^dx" + N''dt){dx^ + N^dt) 



(6.3) 



where and N'^ are respectively the lapse function and shift vector, '-jab is a transverse and 

traceless tensor as = c?a7ab = 0, and e'''"*' = 5ab + lab + lad a 1 2 H . On the fiat hypersurface, 

the gauge-invariant curvature perturbation is related to the first-o rder fiuct n ation of the scalar 
field <y9 as C = ~-f^V^/0[3 Following the conversion equations (D.15) and (D.22), we decompose C, 
and into the helicity states as 



e(°)(k) = C(k), e^^^^(k) 



(k)7a.(k) . 



(6.4) 



Here, is a transverse and traceless polarization tensor explained in Appendix D The inter- 
action parts of this action have been derived by Maldacena [l] up to the third-order terms. In 
particular, we will focus on an interaction between two scalars and a graviton. This is because the 
correlation between a small wave number of the tensor mode and large wave numbers of the scalar 
modes will remain despite the tensor mode decays after the mode reenters the cosmic horizon. We 
find a leading term of the two scalars and a graviton interaction in the action coming from the 
matter part of the Lagrangian through X as 



A|3rd-ordcr D CL —'JabOa'fOb^ 



therefore, the interaction part is given by 



'S'int 3 / d^x agtsslabdaCdbC 



(6.5) 



(6.6) 



^Here, ( and jab are equivalent to TZ and hat in Eq. (2.401, respectively. 



6 CMB BISPECTRUM INDUCED BY THE TWO SCALARS AND A GRAVITON 
80 CORRELATOR 



Here, we introduce a coupling constant gtss- From the definition of C)7afe and tlie slow-roll pa- 
rameter, Qtss = e. For a general consideration, let us deal with gtss as a free parameter. In this 
sense, constraining on this parameter may offer a probe of the nature of inflation and gravity in 
the early Universe. The primordial bispectrum is then computed using in-in formalism in the next 
subsection. 



6.2 Calculation of the initial bispectrum 

In the same manner as discussed in Ref. fl], we calculate the primordial bispectrum generated 
from two scalars and a graviton in the lowest order of the slow-roll parameter: 

3 



(^(±2)(ki)e(°)(k,)e(°)(k3)> = (27r)3F±2°°(ki,k„k3)5 



n=l 
4 



rn J J \ ; , kik2 + k2h + hki kik2h 

/(fcl, fc2, fcs) = -kt H 7 \ ' 

where kf = ki + k2 + k^, and * means that it is evaluated at the time of horizon crossing, i.e., 
a^^H^, = k. Here, we keep the angular and polarization dependences, e^J^''(ki)/c2a^3fe) which have 
sometimes been omitted in the literature for simplicity (4]-[6|. We show, however, that expanding 
this term with spin-weighted spherical harmonics enables us to formulate the rotational-invariant 
bispectrum in an explicit way. The statistically isotropic power spectra of and are 
respectively given by 



(e(°)(k)e(°)*(k')> = {2n)'Psik)6ik-k') 



27r2 87r2c,*e*Mp2 



A, 



i^'\k)e'^*iy)) ^ (27r)3^5(k - k')5A,A' (for A, A' = ±2) , ^^'^^ 

k'PAk) Hi r 

- - = 2^^ ' 

where r is the tensor-to-scalar ratio and As is the amplitude of primordial curvature perturbations 
Note that the power spectra satisfy the scale invariance because we consider them in the lowest 
order of the slow-roll parameter. Using these equations, we parametrize the initial bispectrum in 
this case from Eq. (6.7) as 

F±2»»(ki,k2.k3) = /'^*^'(ii,t2,*:3)eif' (k,)4„4^, (6.9) 
/<-^-)(*„*„*3) - (6.10) 

12 3 ^ 1 



Note that fC^^^) seems not to depend on the tensor-to-scalar ratio . In F ig. 6.1 , we show the shape 
of I /ki- From this, we confirm that the initial bispectrum fC^^^) (6.10) dominates in the squeezed 
limit as ki k2 — k^ like the local-type bispectrum of scalar modes given by Eq. (4.7). 



^For Cs* — 1, these results are identical to Eqs. (2.59) and (2.681. 
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Figure 6.1: Shape of I /ki. For the symmetric property and the triangle condition, we limit the 
plot range as ki < k2 < k^ and \ki — k2\ < k^ < ki + k2- 

In the squeezed limit, the ratio of fC^^^^ to the scalar-scalar-scalar counterpart f^^^^") = 
lf^LPs{ki)Ps{k2), which has been considered frequently, reads 

f^^lOgts^^kth 20gts^l r« Hi 

/(^^^) 3/nl h ki ^ 3/nl h ■ ^ ' 

In the standard slow-roll inflation model, this ratio becomes 0{1) and does not depend on the 
tensor-to-scalar ratio because gtss and /nl are proportional to the slow-roll parameter e, and I/kt 
has a nearly flat shape. The average of amplitude is evaluated as I/kt ~ —0.6537. Therefore, it 
manifests the comparable importance of the higher order correlations of tensor modes to the scalar 
ones in the standard inflation scenario. 



6.3 Formulation of the CMB bispectrum 



Here, using Eqs. (5.5), (6.7) and (6.9), we explicitly calculate the CMB tensor-scalar-scalar bis- 
pectrum as 



(T) (S) (S) 



—2 



1 



n=l 



(2jr)» 



Ai=±2 



n=l 



n / d'^n 



:{27rff'^^''\k,,k2,ks)e'-f\k,)k2ahbSll[kA . (6.12) 



\n=l 



At first, we express all parts containing the angular dependence with the spin spherical har- 
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(6.13) 



y^dy 



^7-0 

^^iii2i3 \ Ml M2 M: 



.i=l LiMi 

Li L2 L3 

3 



where we used the relations hsted in Appendices O and D and 



JS1S2S3 



(2/i + l)(2/2 + l)(2/3 + l) / h k h 
47r V -51 S2 S3 



(6.14) 



(6.15) 



Secondly, using Eq. (C.8), we replace all the integrals of spin spherical harmonics with the Wigner 
symbols: 



rf'ki ^2r,:^^(k;)r;^,,^(ki)±2r2V(ki) = /^^f ( 
d^k; F;^^^(k2)F;^^^(k2)r4jk2) = 



2 

M 



rf^ks r;3„3(k3)F^V3(k3)r4^(k3) = 



i2 L2 1 

m2 M2 ma 

h L3 1 

m3 M3 mfe 



(6.16) 



Thirdly, using the summation formula of five Wigner-3j symbols as Eq. (C.20), we sum up the 
Wigner-3j symbols with respect to azimuthal quantum numbers in the above equations and express 
with the Wigner-9j symbol as 



E 

M1M2M3 

MmalTli, 



X 



Li L2 

Ml M2 M3 

4 Li 2 
mi Ml M 



2 1 1 

M rua rrih 

h L2 1 

1712 M2 ITLa 



h U 1 
m3 M3 mfe 



«3 

mi m2 m3 




(6.17) 



After these treatments and the summation over Ai = ±2 as 



Ai=±2 



Ai 



rAiO-Ai 



1lf^ll2 (for xi + Li + £1 = even) 
(for xi + Li + £1 = odd) 



(6.18) 



we can obtain the CMB angle-averaged bispectrum induced from the nonlinear coupling between 
two scalars and a graviton as 



B 



(TSS) 



in 



,3/2 



E(-i) 



^0 7-20-2 7-0 7-0 
-'LiL2L3-'fiLi2-'£2i2l-'4i3l 



LIL2L3 




^Equations (3.14) and (3.21) in Ref. j2j include typos. 
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n 



n=l 



/(^^^)(fci,^2,A;3)(6.19) 



Note that the absence of the summation over 1711,1712 and ins in this equation means that the 
tensor-scalar-scalar bispectrum maintains the rotational invariance. As described above, this con- 
sequence is derived from the angular dependence in the polarization tensor. Also in vector modes, 
if their power spectra obey the statistical isotropy like Eq. (6.8), one can obtain the rotational in- 
variant bispectrum by considering the angular dependence in the polarization vector as Eq. (D.ll ). 
Considering Eq. (6.18) and the selection rules of the Wigner symbols explained in Appendix |C[ 
we can see that the values of Li, L2 and L3 are limited as 




(for Xi 
(for Xi 



I,E) 
B) 



1^2 ±1| 



1^3 ±1| 



(6.20) 



|Li -L2I < < Li + L2 , 



even 



i=l 



and the bispectrum (6.19) has nonzero value under the conditions: 



i=l 



even (for Xi = J, E) 
odd (for Xi = B) . 



(6.21) 



In Figs. 6^ and 6^ we describe the reduced CMB bispectra of intensity mode sourced from 
two scalars and a graviton coupling: 



dTSS) 



dSST) 



o\-i 



/rO 0\ 



and primordial curvature perturbations (5.7): 



ASSS) 



(jO oy 



q{TSS) ^ 



-1 T^{SSS) 



(STS) 
III 1^1 



+ B 



(SST) 



(6.22) 



(6.23) 



For the numerical computation, we modify the Boltzmann Code for Anisotropies in the Microwave 
Background (CAMB) pf, 8||^ In the calculation of the Wigner-3j and 9j symbols, we use the Com- 
mon M athematical Library SLATEC i9j and the summation formula of three Wigner-6j symbols 
(C.21 ). As the radiation transfer functions of scalar and tensor modes, namely, 7J, ^. and 7J. we 
use Eq. (3.101 ). From the behavior of each line shown in Fig. 6.3 at small £3 that the reduced CMB 
bispectrum is roughly proportional to we can confirm that the tensor-scalar-scalar bispectrum 
has a nearly squeezed-type configuration corresponding to the shape of the initial bispectrum as 
discussed above. From Fig. 02, by comparing the green dashed line with the red solid line roughly 
estimated as 



lATSS) USTS) ,(55T) , 

riiiMi ' "III Me ' "III Mel 



X 10-i«|(7t 



(6.24) 



we find that \gtss\ ~ 5 is comparable to /] 
the PLANCK experiment. 



local 
NL 



5 corresponding to the upper bound expected from 



The CMB bispectra generated from the two scalars and a graviton correlator in Figs. 6.2 and 6.3 become slightly 



smaller than those in Ref. |2j due to the accuracy enhancement of the numerical calculation. 
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Figure 6.2: Absolute values of the CMB reduced bispectra of temperature fluctuation for ii = 
£2 = ^3 = L The lines correspond to the spectra generated from tensor-scalar-scalar correlation 

5 (red solid line) and the primordial non-Gaussianity in the scalar 
5 (green dashed line). The other cosmo logical parameters 
are fixed to the mean values limited from WMAP-7yr data reported in Ref. llOj. 



given by Eq. (6.22) with gt 



curvature perturbations with /^l^' 



6.4 Estimation of the signal-to-noise ratio 



Here, we compute the signal-to-noise ratio by comparing the intensity bispectrum of Eq. (6.19 ) with 



the zero-noise (ideal) data and examine the bound on the absolute value of gtss- The formulation 
of (the square of) the signal-to-noise ratio (S/N) is reported in Refs. [llj and ^12]. In our case, it 
can be expressed as 



E 

2<ii<e2<i3<i 



B 



(TSS) 



MSTS) 



MS ST) 



a 



(6.25) 



where cr^i^a^s denotes the variance of the bispectrum. Assuming the weakly non-Gaussianity, the 



variance can be estimated as 13, 14 



(6.26) 



where takes 1,6 or 2 for ii £2 7^ ^3,^1 = h = ^3, or the case that two £'s are the same, 

respectively. Q denotes that the CMB angular power spectrum included the noise spectrum, which 
is neglected in our case. 



In Fig. 6.4, the numerical result of Eq. (6.25) is presented. We find that (S/N) is a monoton- 
ically increasing function roughly proportional to i for i < 2000. It is compared with the order 



estimation of Eq. (6.25) as Ref. 12 








fSldTSS) ^(STS) 



dSST) I 



]3/2 
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Figure 6.3: Absolute values of the CMB reduced bispectra of temperature fluctuation generated 



from tensor-scalar-scalar correlation given by Eq. (6.22) {TSS + STS + SST) and the primordial 
non-Gaussianity in the scalar curvature perturbations (SSS) as a function of £3 with ii and £2 



fixed to some values as indicated. The parameters are fixed to the same values defined in Fig. 6.2 



^ X 5.4 X 10^ 



\gtss 



Here, we use Eq. (6.24) and the approximations as ^ ~ £^/2A, 
eC, ~ 6 X 10-1° 







(6.27) 



0.36 X £, and 



We confirm that this is consistent with Fig. 6.4, which justifies our numerical 
calculation in some sense. This figure shows that from the WMAP and PLANCK experimental 
data 10, 15 , which are roughly noise-free at £ < 500 and 1000, respectively, expected {S/N)/gtss 
values are 0.05 and 0.11. Hence, to obtain (S/N) > 1, we need \gtss\ > 20 and 9. The latter value 
is close to a naive estimate \gtss\ ^ 5, which was discussed at the end of the previous subsection. 



6.5 Summary and discussion 

In this section, we present a full-sky formalism of the CMB bispectrum sourced from the pri- 
mordial non-Gaussianity not only in the scalar but also in the vector and tensor perturbations. 
As an extension of the previous formalism discussed in Ref. j5], the new formalism contains the 
contribution of the polarization vector and tensor in the initial bispectrum. In Ref. [5], we have 
shown that in the all-sky analysis, the CMB bispectrum of vector or tensor mode cannot be formed 
as a simple angle-averaged bispectrum in the same way as that of scalar mode. This is because 
the angular integrals over the wave number vectors have complexities for the non-orthogonality 
of spin spherical harmonics whose spin values differ from each other if one neglects the angular 
dependence of the polarization vector or tensor. In this study, however, we find that this difficulty 
vanishes if we maintain the angular dependence in the initial bispectrum. 

To present how to use our formalism, we compute the CMB bispectrum induced by the nonlinear 
mode-coupling between the two scalars and a graviton [^. The typical value of the reduced 
bispectrum in temperature fluctuations is calculated as a function of the coupling constant between 
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Figure 6.4: Signal-to-noise ratio normalized by gtss as a function of the maximum value between 
£1,^2 and £3, namely, C.. Each parameter is fixed to the same values defined in Fig. 6.2 



scalars and gravitons (^js^: |&/// £«+&//J«£+&//7 ^ ~ ^~*^x8x 10^^^|(7fss|. Through the computation 
of the signal-to-noise ratio, we find that the two scalars and a graviton coupling can be detected 
by the WMAP and PLANCK experiment if \gtss\ ~ 0(10). Although we do not include the effect 
of the polarization modes in the estimation of gtss in this study, they will provide more beneficial 
information of the nonlinear nature of the early Universe. 
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7 VIOLATION OF THE ROTATIONAL INVARIANCE IN THE CMS BISPECTRUM 



7 Violation of the rotational invariance in the CMB bis- 
pectrum 

The current cosmological observations, particularly Cosmic Microwave Background (CMB), tell 
us that the Universe is almost isotropic, and primordial density fluctuations are almost Gaussian 
random fields. However, in keeping with the progress of the experiments, there have been many 
works that verify the possibility of the small deviation of the statistical isotropy, e.g., the so-called 
"Axis of Evil" . The analyses of the power spectrum by employing the current CMB data suggest 
that the deviation of the statistical isotropy is about 10% at most (e.g. [l]-[6]). Toward more 
precise measurements in future experiments, there are a lot of theoretical discussions about the 
effects of the statistical anisotropy on the CMB power spectrum, [7-11 , e.g., the presence of the 



off-diagonal configuration of the multipoles in the CMB power spectrum, which vanishes in the 
isotropic spectrum. 

As is well known, it might be difficult to explain such statistical anisotropy in the standard 
infiationary scenario. However, recently, there have been several works about the possibility of 
generating the statistically anisotropic primordial density fiuctuations in order to introduce non- 



trivial dynamics of the vector field. 12 24 . In Ref. 14 , the authors considered a modified hybrid 
infiation model where a waterfall field couples not only with an infiaton field but also with a mass- 
less vector field. They have shown that, owing to the effect of fluctuations of the vector fleld, the 
primordial density fluctuations may have a small deviation from the statistical isotropy and also 
the deviation from the Gaussian statistics. If the primordial density fluctuations deviate from the 
Gaussian statistics, they produces the non-zero higher order spectra (corresponding to higher order 
correlation functions), e.g., the bispectrum (3-point function), the trispectrum (4-point function) 



and so on. Hence, in the model presented in Ref. 14 , we can expect that there are characteristic 
signals not only in the CMB power spectrum but also in the CMB bispectrum. 

With these motivations, in this work, we calculate the CMB statistically anisotropic bispec- 
trum sourced from the curvature perturbations generated in the modifled hybrid inflation scenario 



proposed in Ref. 14 , on the basis of the useful formula presented in Ref. 25 . Then, we flnd 



the peculiar conflgurations of the multipoles which never appear in the isotropic bispectrum, like 



off-diagonal components in the CMB power spectrum. These discussions are based on Ref. 26 . 

This section is organized as follows. In the next subsection, we briefly review the inflation 
model where the scalar waterfall fleld couples with the vector fleld and calculate the bispectrum 



of curvature perturbations based on Ref. 14 . In Sec. 7.2 , we give an exact form of the CMB 
statistically anisotropic bispectrum and analyze its behavior by numerical computation. Finally, 
we devote the flnal subsection to the summary and discussion. 

Throughout this section, we obey the deflnition of the Fourier transformation as 



and a normalization as = (SttG) = 1. 



7.1 Statistically-anisotropic non-Gaussianity in curvature perturba- 
tions 

In this subsection, we briefly review the mechanism of generating the statistically anisotropic 
bispectrum induced by primordial curvature perturbations proposed in Ref. 114] , where the authors 



7.1 Statistically-anisotropic non-Gaussianity in curvature perturbations 
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set the system like the hybrid inflation wherein there are two scalar fields: infiaton (j) and waterfall 
field X; and a vector field coupled with a waterfall field. The action is given by 



S 



dx 



-9 



(7.2) 



Here, F^i, = d^A^, — duA^ is the field strength of the vector field A^, V^(0, x, A^) is the potential of 
fields and /(0) denotes a gauge coupling. To guarantee the isotropy of the background Universe, 
we need the condition that the energy density of the vector field is negligible in the total energy 
of the Universe and we assume a small expectation value of the vector field. Therefore, we neglect 
the effect of the vector field on the background dynamics and also the evolution of the fiuctuations 
of the inflaton. In the standard hybrid inflation (only with the inflaton and the waterfall field), 
infiation suddenly ends owing to the tachyonic instability of the waterfall fleld, which is triggered 
when the inflaton reaches a critical value 0e- In the system described using Eq. (7.2), however. 



0e may fluctuate owing to the fluctuation of the vector field and it generates additional curvature 
perturbations. 



Using the 5N formalism 27 -32 , the total curvature perturbation on the uniform-energy-density 
hypersurface at the end of infiation t = tf, can be estimated in terms of the perturbation of the 
e-folding number as 



ate 



6N{t„t 
dN 



5(f)^ + 



2 902 

dN #e(^) 



dN £(l)eiA) d^N d(f)eiA) rf0e(A) 

d^edMA^ ^ ~d^~dM dA^ 



(7.3) 



Here, is the time when the scale of interest crosses the horizon during the slow-roll infiation. As- 
suming the sudden decay of all fields into radiations just after infiation, the curvature perturbations 
on the uniform-energy-density hypersurface become constant after infiation ends Hence, at the 
leading order, the power spectrum and the bispectrum of curvature perturbations are respectively 
derived as 

{2'nfNlP4,{k{)5 (y^K 

\n=l 

,d(t)e{A) d(t)e{A) 



nc(kn) 



\n=l 



+Nt 



(M^:(ki)M:(k2)) 



nc(kn) 



dA^' dA'^ 

{2TxfNlN,,[P^{k{)P^{k2) + 2 perms.]5 ( J^k^ 



(7.4) 



\n=l 



vn=l 



+N: 



.d(t)e{A) d(t)o{A) d(t)e{A) 



dAi' dA'^ 



dAp 



{6A^M6A:{k,)SAak,)) 



4 C;0e(^)#e(^) / 1 ^Vc(^) A^ee d(P^{A) c/0e(^) 



dA'' 



NdAPdA'' m dAp dA'^ 



X [((5A^(ki)M^(k2)(M''^M'^)e(k3)) + 2 perms. 



(7.5) 



where P^{k) = Hl/{2k^) is the power spectrum of the fiuctuations of the infiaton. A"* = dN/d(j)^, 
A'^,* = d'^N/d(j)l, Ae = dN/d(j)e, Aee = d'^N/d(j)l, and -k denotes the convolution. Here, we 
assume that 50* is a Gaussian random field and {6(f)AP') = 0. 



^This C is consistent with TZ in Eq. (2.40) 
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For simplicity, we estimate the fluctuation of the vector fields in the Coulomb gauge: SAq = 
and fcjA* = 0. Then, the evolution equation of the fluctuations of the vector field is given by 



Ai - {a - a^djd^A = 
/ 



(7.6) 



where Ai = fSAi, ' denotes the derivative with respect to the conformal time, and we neglect 
the contribution from the potential term. When / oc a, with appropriate quantization of the 
fluctuations of the vector field, we have the scale-invariant power spectrum of SA^ on superhorizon 



scale as 14, 18 33 



{SAl{k,)6Ai{k,)) = {2iTfP^{k)f-'P^^{k,)6 (^E^nj , (7.7) 
where a is the scale factor, P^^{k) = S^^ — k^k\ " denotes the unit vector, and /e = f{tc)- Therefore, 



substituting this expression into Eq. (7.4), we can rewrite the power spectrum of the primordial 
curvature perturbations, (, as 



nC(k„)) = (27r)3Pc(ki)Wl]kn , 



\n=l 



vn=l 



Pc(k) = 



(7. 



where qt = dcl)e/dA\qij = dPcpd {dA^dA^). We can divide this expression into the isotropic part 
and the anisotropic part as |7| 



Pc(k) ^ Pf(A:) 



l + gp (q-k)' 



with 



Pf\k) = NlP^{k){l + P) , gp 



/3 



1 + /3 



(7.9) 



(7.10) 



where (3 = (iVe/A^*//e)^ |qp. The bispectrum of the primordial curvature perturbation given by 



Eq. (7.5) can be written as 

3 



nC(k„)) = (27r)3p^(ki,k2,k3)5 Ek„ 



(7.11) 



\n=l 



,n=l 



Pc(ki,k2,k3) 



9l 

(3 



Pf (A;i)Pf (A;2) 



^ + (jrf + w^^^' ) ^-(kon.(k.) 



+2 perms. . 



(7.12) 



Here q'^'^ = q'^'^/\ci\'^ and we have assumed that the fluctuation of the vector field 6A^ almost obeys 
Gaussian statistics; hence {5A!^{ki)5A'^{k2)5AP{ks)) = 0. 



7.2 CMB statistically- anisotropic bispectrum 
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Hereinafter, for calculating the CMB bispectrum explicitly, we adopt a simple model whose 



potential looks like an Abelian Higgs model in the unitary gauge as 14 

A 



(7.13) 



where X,g, and h are the coupling constants, m is the inflaton mass, and v is the vacuum expec- 
tation value of X- Since the effective mass squared of the waterfall field is given by 



-\v' + g'^o + h'A'Ai = 



and the critical value of the inflaton 0e can be obtained as 

g'<Pl = \v' - h'A^A, , 



(7.14) 



(7.15) 



we can express (3, g* and g'-' in Eq. (|7.12|) in terms of the model parameters as 

2 



= -A' , q 



— __ 



hA 



6'^ + A'A> 



1 (h'^A^^ 



(7.16) 



where we have used A^* ~ —N^ ^ l/\/2e with e = [dV / dcfy/VY /2 being a slow-roll parameter and 



|A| = A. Substituting these quantities into Eq. (7.12), the bispectrum of primordial curvature 
perturbations is obtained as 



F^(ki,k2,k3) = CPf°(A;0Pr(^2)A"A''5^'^Pac(ki)nd(k2) + 2 perms. 



C 



C^-gt 



2 



f-V 



(7.17) 



Note that in the above expression, we have neglected the effect of the longitudinal polarization in 
the vector field for simplicity |^ and the terms that are suppressed by a slow-roll parameter rj = 
d'^V/d(j? /V because —N^^/N^ ~ Nee/N^ — — (A^c^c)"^ — V- Since the current CMB observations 
suggest gp < C(O.l) (e.g., Refs. [l][2]) and A^"^ ~ — \/2e, the overall amplitude of the bispectrum 
in this model, C, does not seem to be sufficiently large to be detected. However, even if <^ 1 
and e ^ 1, C can become greater than unity for {g/h/AYcj)^ ^ 1. Thus, we expect meaningful 
signals also in the CMB bispectrum. Then, in the next subsection, we closely investigate the 



CMB bispectrum generated from the primordial bispectrum given by Eq. (7.17) and discuss a 



new characteristic feature of the CMB bispectrum induced by the statistical anisotropy of the 
primordial bispectrum. 



7.2 CMB statistically-anisotropic bispectrum 

In this subsection, we give a formula of the CMB bispectrum generated from the primordial 
bispectrum, which has statistical anisotropy owing to the fluctuations of the vector field, given 
by Eq. (7.17). We also discuss the special signals of this CMB bispectrum, which vanish in the 
statistically isotropic bispectrum. 



^Owing to this treatment, we can use the quantities estimated in the Coulomb gauge as Eq. (7.12). In a more 
precise discussion, we should take into account the contribution of the longitudinal mode in the unitary gauge. 
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7.2.1 Formulation 

The CMB fluctuation can be expanded in terms of the spherical harmonic function as 

AX(h) 



X 



(7.18) 



£m 



where n is a unit vector pointing toward a line-of-sight direction, and X denotes the intensity 



(= /) and polarizations (= E,B). According to Eq. (5.3), the coefficient, a^m, generated from 
primordial curvature perturbations, is expressed as 



ax,em = ATi{~iy I ^—QUk)rxAk) (for X = /, 



(27r)3 



(7.19) 
(7.20) 



where Tx/ is the time- integrated transfer function of scalar modes as described in Eq. (3.101). 



Using these equations, the CMB bispectrum generated from the bispectrum of the primordial 
curvature perturbations is given by 



with 



\n=l 



n=l 



(7.21) 



\n=l 



n=l 



(27r)35 K]k„ Fc(ki,k2,k3) . (7.22) 



,n=l 



We expand the angular dependencies which appear in the Dirac delta function, (5(ki + k2 + ks) 



and the function, F^(ki, k2, ks), given by Eq. (7.17) with respect to the spin spherical harmonics 
as 



vn=l 



8 / y'dy 

'0 



n=l LnMn 

Li L2 L3 



rO 

^1^2L3 \^ Ml M2 M3 



(7.23) 



^ / r r ' r , — n 9 ^ 



L,L',La=0,2 
MM' Ma 



L L' La 
M M' Ma 



,(7.24) 



where the 2x3 matrices of a bracket and a curly bracket denote the Wigner-3j and 6j symbols, 
respectively, and 



JS1S2S3 



(2/1 + 1) (2/2 + 1) (2/3 + 1) / /l h h 

4:71 V Si S2 S3 



(7.25) 
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Here, we have used the expressions of an arbitrary unit vector and a projection tenso r as Appendices 
[c]and[D] Note that for Yqq{A) = I/a/Stt, the contribution of = in Eq. ( 7.24 ) is independent 
of the direction of the vector field. Therefore, the statistical anisotropy is generated from the 
signals of La = 2. By integrating these spherical harmonics over each unit vector, the angular 
dependences on ki, k2, ks can be reduced to the Wigner-3j symbols as 



j2-C V* V* V* 
" l.^ra2'' L'M' 



rO 



rO 
^e2L2L' 



mi Ml M 

I2 L2 L' 
ma M2 M' 



(7.26) 



mg-'LaMs — ( — 1)'"'''^L3,^3^M'3 -ma 



From these equations, we obtain an alternative explicit form of the bispectrum of (im as 



' 3 \ „oo 

\n=l I "'0 L1L2 



"-^+■^2+% 
^LiL2h 



xPf\ki)3LAkiy)Pr{h)3L2{k2y)C3,,{hy) 
3 



x4 



47r\ 



-1) 



ma 



y rO rO rOOO 



E 7-01-1 rOl-l 
-'ill -'l'11 

L,L',La=0,2 

L L' La 
1 1 1 



M-i_M2MM'MA 

f ii Li L 
\mi Ml M 

+2 perms. . 



Li L2 

Ml M2 



*3 

-ma 



£2 L2 L' 
m2 M2 M' 



L L' La 

M M' Ma 



(7.27) 



This equation implies that, owing to the vector field A, the CMB bispectrum has a direction 
dependence, and hence, the dependence on mi,m2,m3 cannot be confined only to a Wigner-3j' 
symbol, namely. 



n^^"'»"(^'^) ) ^ {'2T^fJ^t^i2h{ki,k2,k^ 



\n=l 



^1 ^2 
mi 1712 ITT-S 



(7.28) 



This fact truly indicates the violation of the rotational invariance in the bispectrum of the primor- 
dial curvature perturbations and leads to the statistical anisotropy on the CMB bispectrum. 

Let us consider the explicit form of th e CMB bis pectr um. Here, we set the coordinate as 
A = z. Then, by substituting Eq. (7.27) into Eq. (7.21) and using the relation Yj 

•y/ {2La + 1)/ (47r)(5MA,0) the CMB bispectrum is expressed as 



LaMa^ 



n 

Kn=l 



ax, 



n j-^n' ' '"n 



y^dy 



n 



TT 



L1L2 



_n=l 



kldknTx„,i„{kn) 



+£3 rO 
^L^L2h 



94 



7 VIOLATION OF THE ROTATIONAL INVARIANCE IN THE CMS BISPECTRUM 



47r~ 



x4 1 — 1 



ErOl-l rOl-1 
-'ill -'L'11 

L,L',La=0,2 

L V La 



^ rO rO rOOO 
^^hL^L^hLiL'^llLA 111 1 



2La + 1 



All 



E 

M=~2 



-nil - M 



L2 h 

-m2 + M — ms 



nil —rrii — M M 

L L' La 
M -M 



i2 L2 L' 

ni2 -ni2 + M -M 



+ 2 perms. 



(7.29) 



By the selection rules of the Wigner symbols described in AppendixjC} the ranges of £i, £2, ^3, t^ii ^^2 
and ms, and the summation ranges in terms of Li and L2 are limited as 



E = even , ^ m„ = , 

n=l n=l 

Li = \ii-2\,ii,ii + 2 , L2 = |4-2U2,^2 + 2 

1^2-41 < ^1 < i^2 + 4 • 



(7.30) 



7.2.2 Behavior of the CMB statistically-anisotropic bispectrum 



On the basis of Eq. (7.29), we compute the CMB bispectra for the several £'s and m's. Then, we 



modify the Boltzmann Code for Anisotropics in the Microwave Background (CAMB) |34, 35 and 
use the Common Mathematical Library SLATEC 



36 



In Fig. |7.1[ the red solid lines are the CMB statistically anisotropic bispectra of the intensity 

L, and the green dashed hues are the statistically isotropic one 



mode given by Eq. (7.29) with C 



sourced from the local-type non-Gaussianity of curvature perturbations given by Eq. (4.7) 



n 

\n=l 



jO / "^1 
hh^s \ jYi^ ^2 ni3 



6 



y'^dy 



n 



n=l 



TT 



kldknTx„/r^{kn)ji„{kny) 



X Pr{ki)Pl'°{k2)-hL + 2 perms. 



(7.31) 



with /nl = 5 for £1 = £2 = ^3 and two sets of nii,ni2,ni3. From this figure, we can see that the 
red solid lines are in good agreement with the green dashed line in the dependence on i for both 
configurations of mi,ni2,ni3. This seems to be because the bispectrum of primordial curvature 



perturbations affected by the fluctuations of vector field given by Eq. (7.17) has not only the 



anisotropic part but also the isotropic part and both parts have the same amplitude. In this sense, 
it is expected that the angular dependence on the vector field A does not contribute much to a 
change in the shape of the CMB bispectrum. We also find that the anisotropic bispectrum for 
C ~ 0.3 is comparable in magnitude to the case with /nl = 5 for the standard local type, which 
corresponds to the upper bound on the local-type non-Gaussianity expected from the PLANCK 



experiment 37 
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Figure 7.1: Absolute values of the CMB statistically anisotropic bispectrum of the intensity mode 



given by Eq. ( |7.29[ ) with C 
with /nl 



1 (red solid line) and the statistically isotropic one given by Eq. (7.31 ) 
5 (green dashed line) for h = £2 = The left and right figures are plotted in the 
configurations (mi, 1712, m^) = (0, 0, 0), (10, 20, —30), respectively. The parameters are fixed to the 



mean values limited from the WMAP-7yr data as reported in Ref. 38 . 



In the discussion of the CMB power spectrum, if the rotational invariance is violated in the 
primordial power spectrum given by Eq. (7.9), the signals in the off-diagonal configurations of i also 
have nonzero values [TIIS, 10 . Likewise, there are special configurations in the CMB bispectrum 



induced from the statistical anisotropy on the primordial bispectrum as Eq. (7.17). The selection 



rule (7.30) suggests that the statistically anisotropic bispectrum (7.29) could be nonzero in the 



multipole configurations given by 



ii = \i2-h\- 4, \i2 - 41 - 2, £2 + 4 + 2, £2 + 4 + 4 



(7.32) 



and two permutations of £1,^2, 4- In contrast, in these configurations, the isotropic bispectrum 



e.g., Eq. (7.31 )) vanishes owing to the triangle condition of the Wigner-3j symbol 



nil ^2 ^3 



and the nonzero components arise only from 



\i2-i3 



+ 



(7.33) 



Therefore, the signals of the configurations (7.32) have the pure information of the statistical 
anisotropy on the CMB bispectrum. 



Figure 7.2 shows the CMB anisotropic bispectra of the intensity mode given by Eq. (7.29) 
with C = 1 for the several configurations of €s and m's as a function of £3. The red solid line 
and green dashed line satisfy the special relation (7.32), namely, £1 = i'2 + 4 + 2, [£2 — 4| — 2, 
and the blue dotted line obeys a configuration of Eq. (7.33), namely, £1 = 



+ 



From this 



figure, we confirm that the signals in the special configuration (7.32) are comparable in magnitude 
to those for £1 = ^ + 4- Therefore, if the rotational invariance is violated on the primordial 
bispectrum of curvature perturbations, the signals for £1 = £2 + 4 + 2, 1^2 — £3 1 — 2 can also become 
beneficial observables. Here, note that the anisotropic bispectra in the other special configurations: 
C-i = 4 + 4 + 4, 1^2 — £3! — 4 are zero. It is because these signals arise from only the contribution of 
L = L' = = 2, Li = £1 ± 2, L2 = -^2 ± 2 in Eq. (7.29) owing to the selection rules of the Wigner 
symbols, and the summation of the four Wigner-3j symbols over M vanishes for all £'s and m's. 
Hence, in this anisotropic bispectrum, the additional signals arise from only two configurations 
^1 = 4 + 4 + 2, 1^2 — 41 — 2 and these two permutations. 
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Figure 7.2: Absolute values of the CMB statistically anisotropic bispectra of the intensity mode 
given by Eq. (7.29) for (mi,m2,m3) = (0,0,0) (left panel) and (10,20,-30) (right one) as the 
function with respect to £3. The lines correspond to the spectra for (-^1,^2) = (102 + £3, 100) (red 
solid line), (1100-41-2, 100) (green dashed line) and (100 + £3, 100) (blue dotted line). The 
parameters are identical to the values defined in Fig. |7.1[ 



7.3 Summary and discussion 



In this section, we investigated the statistical anisotropy in the CMB bispectrum by considering 



the modified hybrid inflation model where the waterfall field also couples with the vector field 14 
We calculated the CMB bispectrum sourced from the non-Gaussianity of curvature perturbations 
affected by the vector field. In this inflation model, owing to the dependence on the direction of the 
vector field, the correlations of the curvature perturbations violate the rotational invariance. Then, 
interestingly, even if the magnitude of the parameter (7^ characterizing the statistical anisotropy 
of the CMB power spectrum is too small, the amplitude of the non-Gaussianity can become large 
depending on several coupling constants of the fields. 



Following the procedure of Sec. |5j j25], we formulated the statistically anisotropic CMB bispec- 
trum and confirm that three azimuthal quantum numbers mi, m2 and 7/13 are not confined only to 
/ £2 4 \ 

the Wigner symbol . This is evidence that the rotational invariance is violated 

Y mi m2 m3 J 

in the CMB bispectrum and implies the existence of the signals not obeying the triangle condition 
of the above Wigner symbol as {£2 — 4| < ^1 < ^2 + h- We demonstrated that the signals of the 
CMB bispectrum for £1 = ^2 + 4 + 2, |4 — 4| —2 and these two permutations do not vanish. In fact, 
the statistically isotropic bispectra are exactly zero for these configurations; hence, these signals 
have the pure information of the statistical anisotropy. Because the amplitudes of these intensity 
bispectra are comparable to those for £1 = £2 + ^3, it might be possible to detect these contributions 
of the statistical anisotropy in future experiments, which would give us novel information about 
the physics of the early Universe. Of course, also for the ii^-mode polarization, we can give the 
same discussions and results. 

Although we assume a specific potential of inflation to show the statistical anisotropy on the 
CMB bispectrum explicitly, the above calculation and discussion will be applicable to other infla- 
tion models where the rotational invariance violates. 
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8 PARITY VIOLATION OF GRAVITONS IN THE CMB BISPECTRUM 



8 Parity violation of gravitons in the CMB bispectrum 

Non-Gaussian features in the cosmological perturbations include detailed information on the nature 
of the early Universe, and there have been many works that attempt to extract them from the 
bispectrum (three-point function) of the cosmic microwave background (CMB) anisotropics (e.g., 
Refs. (l}]4]). However, most of these discussions are limited in the cases that the scalar-mode 
contribution dominates in the non-Gaussianity and also are based on the assumption of rotational 
invariance and parity conservation. 

In contrast, there are several studies on the non-Gaussianities of not only the scalar-mode 
perturbations but also the vector- and tensor-mode perturbations (5]-[7]. These sources produce the 
additional signals on the CMB bispectrum [s] and can give a dominant contribution by considering 
such highly non-Gaussian sources as the stochastic magnetic fields j9j. Furthermore, even in the 
CMB bispectrum induced from the scalar-mode non-Gaussianity, if the rotational invariance is 



violated in the non-Gaussianity, the characteristic signals appear 10 . Thus, it is very important 



to clarify these less-noted signals to understand the precise picture of the early Universe. 



Recently, the parity violation in the graviton non-Gaussianities has been discussed in Refs. 11 



12 . Maldacena and Pimentel first calculated the primordial bispectrum of the gravitons sourced 

cubic terms, namely. 
Soda et al. proved 

that the parity-violating non-Gaussianity of the primordial gravitational waves induced from WW^ 
emerges not in the exact de-Sitter space-time but in the quasi de-Sitter space-time, and hence, its 



from parity-even (parity-conserving) and parity-odd (parity-violating) Wey 
and WW"^, respectively, by making use of the spinor helicity formalism 



11 



amplitude is proportional to a slow-roll parameter. 12 In these studies, the authors assume that 
the coupling constant of the Weyl cubic terms is independent of time. 

In this section, we estimate the primordial non-Gaussianities of gravitons generated from 
and WW"^ with the time-dependent coupling parameter [isj. We consider the case where the 
coupling is given by a power of the conformal time. We show that in such a model, the parity 
violation in the non-Gaussianity of the primordial gravitational waves would not vanish even in the 
exact de-Sitter space-time. The effects of the parity violation on the CMB power spectrum have 
been well-studied, where an attractive result is that the cross-correlation between the intensity 



and i?-mode polarization is generated 14 17 . On the other hand, in the CMB bispectrum, owing 



to the mathematical property of the spherical harmonic function, the parity-even and parity-odd 



signals should arise from just the opposite configurations of multipoles 18 , 19 . Then, we formulate 
and numerically calculate the CMB bispectra induced by these non-Gaussianities that contain all 
the correlations between the intensity (/) and polarizations {E, B) and show that the signals from 
(parity-conserving) appear in the configuration of the multipoles where those from WW"^ 



(parity-violating) vanish and vice versa. These discussions are based on Ref. 20 . 

This section is organized as follows. In the next subsection, we derive the primordial bispectrum 
of gravitons induced by and WW^ with the coupling constant proportional to the power 



of the conformal time. In Sec. 8.2, we calculate the CMB bispectra sourced from these non- 
Gaussianities, analyze their behavior and find some peculiar signatures of the parity violation. 
The final subsection is devoted to summary and discussion. In Appendices |E] and [F| we describe 
the detailed calculations of the contractions of the polarization tensors and unit vectors, and of 
the intial bispectra by the in-in formalism. 

Throughout this section, we use Mpi = 1/ y/SrcG, where G is the Newton constant and the rule 
that all the Greek characters and alphabets run from to 3 and from 1 to 3, respectively. 



8.1 Parity-even and -odd non-Gaussianity of gravitons 
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8.1 Parity-even and -odd non-Gaussianity of gravitons 

In this subsection, we formulate the primordial non-Gaussianity of gravitons generated from the 
Weyl cubic terms with the running coupling constant as a function of a conformal time, /(r), 
whose action is given by 



S 



drd X 



3 fir) 



A2 



-gW^ + WW' 



•1) 



with 



W' = W^sW^'^pWP^p , 
WW^ = e'^^^^'W^u^sW^'^pW^p , 



•2) 



where W^^^s denotes the Weyl tensor, e"'^'^'^ is a 4D Levi-Civita tensor normalized as e'^^^'^ = 1, 
and A is a scale that sets the value of the higher derivative corrections 11 . Note that W^ and 
WW^ have the even and odd parities, respectively. In the following discussion, we assume that 
the coupling constant is given by 



fir) 



(8.3) 



where r is a conformal time. Here, we have set /(r*) = 1. Such a coupling can be readily realized 



by considering a dilaton-like coupling in the slow-roll inflation as discussed in Sec. 8.1.2 



8.1.1 Calculation of the primordial bispectrum 

Here, let us focus on the calculation of the primordial bispectrum induced by W^ and WW^ 



Refs. 11, 12 



of Eq. (8.1) on the exact de-Sitter space-time in a more straightforward manner than those of 



At first, we consider the tensor perturbations on the Friedmann-Lemaitre- Robertson- Walker 
metric as 



•4) 



where a denotes the scale factor and jij obeys the transverse traceless conditions; ju = d'jij/dx^ = 
Up to the second order, even if the action includes the Weyl cubic terms given by Eq. (8.1) 



11 



12 



the gravitational wave obeys the action as 

S = I dTdx^a^{'yij%j - ^ij,klij,k) 



.5) 



where ' = d/dr and = d/dx'^. We expand the gravitational wave with a transverse and traceless 
polarization tensor e^-^^ and the creation and annihilation operators a 



(A)t a(A) 



as 



7ij(x,r) 



^,J2^^sik,r)a^^'el^ik)e^^- + h.c. 



X=±2 



^jij is identical to hij in Sec. 2.6 
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(2^ 



ik-x 



A=±2 



with 



7W(k,r) 



(8.6) 



(8.7) 



Here, A = ±2 denotes the hel icity o f the gravitational wave and we use the polarization tensor 
satisfying the relations as Eq. (D.13). The creation and annihilation operators obey the 

relations as 



4'^ |0) = 



"k ' "k' 



(27r)35(k - k')<5A,A' , 



where |0) denotes a vacuum eigenstate. Then, the mode function of gravitons on the de Sitter 
space-time 7^5 satisfies the field equation as 

2 



7ds ids + k 'jds = 

T 



and a solution is given by 



Ids 



H e 



-ikr 



[l + ikr) 



.9) 



(8.10) 



where H = — (ar)~^ is the Hubble parameter and has a constant value in the exact de Sitter 
space-time. 

On the basis of the in- in formalism (see, e.g., Refs. [5} [2l]) and the above results, we calculate 
the tree-level bispectrum of gravitons on the late-time limit. According to this formalism, the 
expectation value of an operator depending on time in the interaction picture, 0{t), is written as 







J^^iJ H,ntit')dt' Q^^^rjn^-iJ Hi„t{t')dt' 







(8.11) 



where T and T are respectively time-ordering and anti-time-ordering operators and Hint{t) is the 
interaction Hamiltonian. Applying this equation, the primordial bispectrum of gravitons at the 
tree level can be expressed as 



\n=l 



dr' ( 



n=l 







(8.12) 



where : : denotes normal product. 

Up to the first order with respect to 7jj, the nonzero components of the Weyl tensor are written 

as 

1 



{HTf{iki,j - ikj,i) 

jk = -{Hrfiiikj - i^j,k) 



(8.13) 



'{.Ht^ ( ^ik'~1 jl,aa ~l~ jk,aa ~l~ ^ jk'~1il,aa ^ jl'~1ik^aoi) 



8.1 Parity-even and -odd non-Gaussianity of gravitons 
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where 'jij^aa = lij + ^"^lij- Then and WW"^ respectively reduce to 



(8.14) 



where 77*-''^ = ^^^^ . Using the above expressions and / drHint = —Smt, up to the third order, the 
interaction Hamiltonians of and WW'^ are respectively given by 



W3 



X ■^1ij,aa [ljk,l3l3lki,aa- + Qikl,iikl,j + Qiik,lijl,k " 127ife,«7M,i] 



(8.15) 



X 77*-^^ [7fcg,aa(-37j.m,/3/37jg,m + 7mj,;9/37mg,j) + '^'ipj,kipm,l{'iil,m " 7im,/)] 



Substituting the above expressions into Eq. (8.12), using the solution given by Eq. (8.10), and 
considering the late-time limit as r — > 0, we can obtain an explicit form of the primordial bispectra: 



wi 



tlfl 



n^^'"Hk„) ) = (2vr)35 K^kJ /i^;(fc,,A;2, Wi1(ki,k;,k3) 



\n=l 



int 



yn=l 



(8.16) 



(r) 
iy3 



(r) 
(a) 



H 



pi 



Re 



Aa) _ (-Ai) 
J ij 



iJ-A2) (-A3) 
2 jfc '^ki 



+ ^(^kl 



^kl «'2j^3j 



3 



_i_"„(-A2) (-As);, 7, -1 



3 (- 



(-A2)^(-A3) 



-kl 1^211^3 j 



+ 5 perms. 



H 



pi 



Im 



z?7 



ijk 



,(-Ai) 

-'kq 



(-A2)„(-A3)/: , .(-A2) (-A3)l 

"-cim "^jni "-'rraij ""oj 



I p(-Ai)„(-A2)£ _ , _ p(-^3)yf 



^31 



+ 5 perms . 



(8.17) 



Here, h = En=i = and VTW^^ "5 perms" denotes the five symmetric terms under 

the permutations of (ki, Ai), (k2, A2), and (ka, A3). From the above expressions, we find that the 
bispectra of the primordial gravitational wave induced from and WW"^ are proportional to 
the real and imaginary parts of t~"^ c/rV'^^'^e^*'^*'^', respectively. This difference comes from 
the number of ■jij^aa and ■jij,k- consists of the products of an odd number of the former 



'Here, we set that r* < 0. 
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terms and an even number of the latter terms. On the other hand, in H^yy2, the situation is 
the opposite. Since the former and latter terms contain 7^5 — k'^jds = {2Ht' /Mpi)k^^'^e~^^'^ and 
ids = i{HT' / Mp\)k^f'^e~^^'^\ respectively, the total numbers of i are different in each time integral. 
Hence, the contributions of the real and imaginary parts roll upside down in /^g and /~^2 • Since 
the time integral in the bispectra can be analytically evaluated as 

r(6 + A)fc-6(-V*)"^ , (8.18) 

/^s and f^^^^ reduce to 

where T{x) is the Gamma function. For more detailed derivation of the graviton bispectrum, see 
Appendix |Fj 

From this equation, we can see that in the casj3j3f the time-independent coupling, which cor- 
responds to the A = case, the bispectrum from WW"^ vanishes. This is consistent with a claim 
in Ref. [l2| On the other hand, interestingly, if A deviates from 0, it is possible to realize the 
nonzero bispectrum induced from WW^ even in the exact de Sitter limit. Thus, we expect the 
signals from WW'^ without the slow-roll suppression, which can be comparable to those from 
and become sufficiently large to observe in the CMB. 



cIt't 



COS I + isin {^A 



8.1.2 Running coupling constant 

Here, we discuss how to realize f oc r"^ within the framework of the standard slow-roll inflation. 
During the standard slow-roll inflation, the equation of motion of the scalar fleld 0, which has a 
potential V, is expressed as 

~ ±v^Mpir-i , (8.21) 

where = [dV/d(j)/ (3Mpiif^)]^/2 is a slow-roll parameter for 0, -|- and — signs are taken to be for 
dV/d(j) > and dV/dcj) < 0, respectively, and we have assumed that aH = — 1/r. The solution of 
the above equation is given by 

= 0* ± v/2e^M"piln (^-^^ . (8.22) 
Hence, if we assume a dilaton-like coupling as / = e^'^~'^*^^^'^ , we have 



T 



A 



/(r)=(^-J , A = ±v/2^^, (8.23) 

where M is an arbitrary energy scale. Let us take r^, to be a time when the scale of the present 
horizon of the Universe exits the horizon during inflation, namely, |r*| = ~ 14Gpc. Then, the 



In Ref. 



12 



, the authors have shown that for ^ = 0, the bispectrum from WW^ has a nonzero value upward 



in the first order of the slow-roll parameter. 



8.2 CMB parity-even and -odd bispectrum 
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coupling /, which determines the amphtude of the bispectrum of the primordial gravitational wave 
induced from the Weyl cubic terms, is on the order of unity for the current cosmological scales. 



From Eq. ( |8.23[ ), we have A = ±1/2 with M = ./Se^Mpi. As seen in Eqs. ( |8.19[ ) and ( |8.2Q[ ), 
this leads to an interesting situation that the bispectra from and WW^ have a comparable 



±/~|^2. Hence, we can expect that in the CMB bispectrum, the signals from 



magnitude as /^g 
these terms are almost the same. 

In the next subsection, we demonstrate these through the explicit calculation of the CMB 
bispectra. 



8.2 CMB parity-even and -odd bispectrum 

In this subsection, following the calculation approach discussed in Sec. [6| we formulai^ the CMB 
bispectrum induced from the non-Gaussianities of gravitons sourced by and WW'^ terms 
discussed in the previous subsection. 



8.2.1 Formulation 

Conventionally, the CMB fluctuation is expanded with the spherical harmonics as 

AX(n) 



X 



(8.24) 



where n is a unit vector pointing toward a line-of-sight direction, and X means the intensity 
(= J) and the electric and magnetic polarization modes (= B). By performing the line-of-sight 
integration, the coefficient, a^m, generated from the primordial fluctuation of gravitons, 7*-^^'', is 



given by [corresponding to Eq. (5.3) 

47r(- 



O-Xlm 



Txm 



(2vr 



\=±2 



(8.25) 
(8.26) 



where x discriminates the parity of three modes: x = for X = I,E and x = 1 for X = B, and 



Tx e is the time-integrated transfer function of tensor modes (3.101). Like Eq. (5.5), we can obtain 



the CMB bispectrum generated from the primordial bispectrum of gravitons as 




n jt-n 'I ''71 



n=l 



kldk 



h 



A„=±2 



(8.27) 



\n=l 



In order to derive an explicit form of this CMB bispectrum, at first, we need to express all the 
functions containing the angular dependence on the wave number vectors with the spin spherical 



harmonics. Using the results of Appendix lEl and f—'^ "^^^ calculated as 



L', L" =2,3 M,M',M" 



2 L' L" 
M M' M" 
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X Ai (kl ) A2 ^L'M' (^2 ) As ^L"M" (^3 ) 

1 /7 



X 



20 V 3 



X 



{a) 



(Stt) 



3/2 




+27r 



+ 5 perms , 



(8.28) 



L',L"=2,3M,M',M" 

X Ai (ki ) A2 Yl'M' (^2 ) A3 Yl"M" (ka ) ( - 1) /^?°2 




27r /7 f 2 L' L" 
12 2^ 

(8.29) 

where the 2x3 matrix of a bracket, and the 2x3 and 3x3 matrices of a curly bracket denote 
the Wigner-3j, 6j and 9j symbols, respectively, and 



JS1S2S3 



(2/1 + 1) (2/2 + 1) (2/3 



h h h 

Sl S2 S3 



(8.30) 



The delta function is also expanded as 

3 



y^dy 
x/°°° 



n=l LnMn 

Li L2 L3 



'ili2L3 \^ Ml M2 Ms 

Next, we integrate all the spin spherical harmonics over ki, k2, ks as 



(8.31) 
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(8.32) 



Through the summation over the azimuthal quantum numbers, the product of the above five 
Wigner-3_7' symbols is expressed with the Wigncr-9j symbols as 
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M1M2M3 

MM'M" 



Li L2 Ls 
Ml M2 M3 
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X 



V 



nil Ml M I \ m2 Ms M' I \ M3 M" 



L" 



^1 ^2 ^3 

mi 1712 ITT'S 




(8.33) 



Finally, performing the summation over the helicities, namely Ai, A2 and A3, as 



A=±2 



E 

A=±2 

A^" 
2 



7-AO-A 
HL2 



E 

A=±2 



A 



rAO-A rAO-A 
^ILL' -'l'12 
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2J^^-^ {£ + L + x = even) 



(8.34) 



and considering the selection rules of the Wigner symbols as described in Appendix [C} we derive 
the CMB bispectrum generated from the non-Gaussianity of gravitons induced by as 



\n=l I vi/3 
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(8.35) 



and WW'^ as 
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X 



2^ 



3,' 



X 



I 7-20-2 
"r-'L'12 



X 




+ 5 perms 



(8.36) 



Here, "5 perms" denotes the five symmetric terms under the permutations of (£1, mi,xi), (£2, ^^2, X2] 
and (£3, ms, X3), and we introduce the filter functions as 



(8.37) 



where the superscripts (e) and (o) denote = even and = odd, respectively. From Eqs. (8.35) 



in a Wigner-3j symbol as 



. This guarantees the rotational invariance of the 



and (|8.36|), we can see that the azimuthal quantum numbers mi, 1712, and are confined only 

il i2 h 

nil ^2 ^3 

CMB bispectrum. Therefore, this bispectrum survives if the triangle inequality is satisfied as 

\il-i2\<i3<il+i2- 



Considering the products between the V functions in Eq. (8.35) and the selection rules as 



Y2n=i ~ 6ven, we can notice that the CMB bispectrum from W does not vanish only for 



J2i^n + Xn) = eyen . (8.38) 



n=l 

Therefore, contributes the ///, HE, lEE, IBB, EEE, and EBB spectra for ^^^^ = even 
and the I IB, lEB, EEB, and BBB spectra for X]n=i ~ odd. This property can arise from any 
sources keeping the parity invariance such as_W^. On the other hand, in the same manner, we 
understand that the CMB bispectrum from WW^ survives only for 

3 

J2i^n + Xn) = odd . (8.39) 

n=l 

By these constraints, we find that in reverse, WW'^ generates the IIB,IEB,EEB, and BBB 
spectra for Yln=i = even and the ///, HE, I EE, IBB, EEE, and EBB spectra for Yln=i = 



odd. This is a characteristic signature of the parity violation as mentioned in Refs. 18, 19]. Hence, 
if we analyze the information of the CMB bispectrum not only for Yln=i = even but also for 
Yln=i = odd, it may be possible to check the parity violation at the level of the three-point 
correlation. 

The above discussion about the multipole configurations of the CMB bispectra can be easily 
understood only if one consider the parity transformation of the CMB intensity and polarization 
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fields in the real space ( |8.24[ ). The ///, HE, LEE, IBB, EEE and EBB spectra from W^, and 
the J/S, lEB, EEB, and BBB spectra from WW'^ have even parity, namely, 



n 



AX,(rii) 



n 



AX,(-rii) 



140) 



Then, from the multipole expansion (8.24) and its parity flip version as 
AX(-h) 



X 



ax,imyim{-ii-) = y^^i-lYax/myemjii) 



(8.41) 



em 



one can notice that Yln=i = even must be satisfied. On the other hand, since the I IB, lEB, EEB, 
and BBB spectra from W^, and the III, HE, lEE, IBB, EEE, and EBB spectra from WW^ 



have odd parity, namely. 



n 

one can obtain J2n=i = odd. 



AXAui) 



X, 



n 



14 = 1 



AX,(-rii) 



(8.42) 



In Sec. 8.2.3, we compute the CMB bispectra (8.35) and (8.36) when A = ±1/2,0,1, that is, 
the signals from become as large as those from WW"^ and either signals vanish. 



8.2.2 Evaluation of f^l and /^''^ 



Here, to compute the CMB bispectra (8.35) and (8.36) in finite time, we express the radial func- 



(r) (r) 

tions, and /^^^a' with some terms of the power of ki,k2, and ^3. Let us focus on the dependence 



on ki,k2, and /cs in Eqs. (8.19) and (8.20) as 



f^'''^ CY f^''^ nr k-^(-kT ]-^- SA{ki,k2,k^ 



where we define Sa to satisfy Sa oc k ^ as 

SAiki,k2,k2) 



{k^k2k^YI\-T,Y 

{hk2k,yi^ 



(8.43) 



(8.44) 



In Fig. 8.1, we plot Sa for A = —1/2,0, 1/2, and 1. From this, we notice that the shapes of Sa 



are similar to the equilateral-type configuration as Eq. (4.10) 

1 1 



'S'cquil(fcl, ^2, ^3) 



6 



J,3i,3 1,3 i,3 i,3j,3 t2L2iL2 
12 3 3 1 12 3 



1 



1 2 
1 



+ 



^1 ^2 ^3 ^1 ^3 ^2 ^2 ^3 ^1 ^2 ^1 ^3 ^3^1 ^2 ^3 ^2 ^1 



. (8.45) 



To evaluate how a function S* is similar in shape to a function S", we introduce a correlation 
function as [3, 22 



cos(5' • S') = 



S-S' 



(5- 5)1/2(5'. 5')l/2 



(8.46) 
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(8.47) 



where the summation is performed over all ki, which form a triangle and P{k) oc k~^ denotes the 
power spectrum. This correlation function gets to 1 when S = S'. In our case, this is calculated 
as 

ro.968 , (A = -1/2) 
0.970 , {A = 0) 
0.971 , {A = 1/2) 
,0.972, (A=l) 



cos(S'a ■ Si 



cquil ) 



< 



(8.48) 



that is, an approximation that 5*^ is proportional to Sgquii seems to be vahd. Here, we also calculate 
the correlation functions with the local- and orthogonal-type non-Gaussianities |4] and conclude 
that these contributions are negligible. Thus, we determine the proportionality coefficient as 
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Sa ■ S, 



equil 



'S'equil ■ 'S'equil 
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equil 
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(8.49) 
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Substituting this into Eqs. (8.19) and (8.20), we obtain reasonable formulae of the radial functions 
for A = 1 /2 as 



fir) ^ fir) 

and for A = — 1 /2 as 



TT 



^^M„,\^ 10395 1.42 X lO^^^equii 
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2 {-T^Y/^hhhyi^ 
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Here, we also use 



rr \ 2 2 
H \ TT^ 



pi 



(8.52) 



where the amplitude of primordial curvature perturbations and r is the tensor-to-scalar 



ratio ji, 8 
function o 



. For A = 0, the signals from WW'^ disappear as f~ ^ = and the finite radial 
is given by 



(r) 
W3 



-rAc 



A 



960 X 2.50 X 10-^S'equii • 



(8.53) 



In contrast, for A = 1, since = 0, we have only the parity-violating contribution from WW'^ 



as 



TT 



-rAc 



A 



pM .V«n 8.09 X lO-^^equil 

' 5760 X 



(8.54) 



8.2.3 Results 



On the basis of the analytical formulae (8.35), (8.36), (8.50), (8.51), (8.53) and (8.54), we compute 
the CMB bispectra from and WW^ for A = —1/2,0,1/2, and 1. Then, we modify the 



Boltzmann Code for Anisotropics in the Microwave Background (CAMB) 231 24 . In calculating 



the Wigner symbols, we use the Common Mathematical Library SLATEC [25] and some analytic 
formulae described in Appendices C and D From the dependence of the radial functions and 



/~:|^2 on the wave numbers, we can see that the shapes of the CMB bispectra from and WW^ 
are similar to the equilateral-type configuration. Then, the significant signals arise from multipoles 
satisfying £i ~ £2 — ^3- We confirm this by calculating the CMB bispectrum for several €s. Hence, 
in the following discussion, we give the discussion with the spectra for £1 ~ £2 — ^3- However, we 
do not focus on the spectra from Yln=i = odd for ii = £2 = h because these vanish due to the 
asymmetric nature. 



In Fig. 8.2, we present the reduced CMB III,IIB,IBB, and EBB spectra given by 



5x1X2X3/1^2^3 



mini2m3 




OX„,£„m„ 



.55) 



n=l 
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Figure 8.2: Absolute values of the CMB III,IIB,IBB, and BBB spectra induced by and 
WW^ for A = —1/2, 0, 1/2, and 1. We set that three multipoles have identical values as £i — 2 = 
^2 — 1 = ^3- The left figures show the spectra not vanishing for Yln=i = even (parity-even 
mode) and the right ones present the spectra for Yln=i^n = odd (parity-odd mode). Here, we fix 
the parameters as A = 3 x lO^GeV, r = 0.1, and = —k^^ = — 14Gpc, and other cosmological 
parameters are fixed as the mean values limited from the WMAP 7-yr data [4]. 
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for 



- 1 



G 



Here, the G symbol is defined by [To] 



2v/4(4 + 1)^2(4 



^i(^i + l)- ^2(^2 + 1) 




-1 1 



(8.57) 



At first, from this figure, we can confirm that there are similar features of the CMB power spectrum 
of tensor modes 



26 27 



In the III spectra, the dominant signals are located in £3 < 100 due to 
the enhancement of the integrated Sachs- Wolfe effect. On the other hand, since the fluctuation of 
polarizations is mainly produced through Thomson scattering at around the recombination and 
reionization epoch, the BBB spectra have two peaks for £3 < 10 and £3 ~ 100, respectively. 
The cross-correlated bispectra between / and B modes seem to contain both these effects. These 
features back up the consistency of our calculation. 

= —1/2,0,1/2, and 1, respectively. We notice 



The curves in Fig. 8.2 denote the spectra for A 



that the spectra for large A become red compared with those for small A. The difference in tilt 

of i between these spectra is just one corresponding to the difference in A. The curves of the left 

^3 



8.2.1 



we 



and right figures obey ^^=1 = even and = odd, respectively. As mentioned in Sec 
stress ^again that in the £ configuration where the bispectrum from vanishes, the bispectrum 
from WW'^ survives, and vice versa for each correlation. This is because the parities of these 
terms are opposite each other. For example, this predicts a nonzero III spectrum not only for 
En=i = even due to but also for XlLi = odd due to WW^. 

We can also see that each bispectrum induced by has a different shape from that induced 
by WW^ corresponding to the difference in the primordial bispectra. Regardless of this, the overall 
amplitudes of the spectra for A = ±1/2 are almost identical. However, if we consider A deviating 
from these values, the balance between the contributions of and WW^ breaks. For example, 
if —1/2 < A < 1/2, the contribution of dominates. Assuming the time-independent coupling, 
namely, A = 0, since /^|y2 ~ 0' CMB bispectra are generated only from W^. Thus, we 
will never observe the parity violation of gravitons in the CMB bispectrum. On the other hand, 
when —3/2 < A < —1/2 or 1/2 < A < 3/2, the contribution of WW^ dominates. In an extreme 
case, if A 

invariance. Then, the information of the signals under Yln=i 
in the analysis of the III spectrum. 

we focus on the /// spectra from for 



odd, since f^j^l = 0, the CMB bispect ra arise only from WW^ and violate the parity 

odd will become more important 



In Fig. 8.3 



'1 = ^2 = ^3 = ^ to compare these with 
the III spectrum generated from the equilateral-type non-Gaussianity of curvature perturbations 
given by 



dSSS) 



y^dy 



n 

n=l 



kldknTj-jl{kn)jl„{kny) 



x-f^^\27r^AsfS,^UkiMM) 



■^The conventional expression of the CMB-reduced bispectrum as 



h h h 

nil 1^2 ^3 



\-i 



(8.58) 



(8.56) 



Here, replacing the / symbol with 




breaks down for — odd due to the divergence behavior of (/°^ 

the G symbol, this problem is avoided. Of course, for J2n=i = even, Gi^i^i^ is identical to li^i^i^ 
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Figure 8.3: Absolute value of the CMB III spectra generated from for A = —1/2 (red solid 
line), (green dashed one) an d 1/2 (blue dotted one), and generated from the equilateral-type 
non-Gaussianity given by Eq. ( [S^SSt with = 300 (magenta dot-dashed one). We set that 

three multipoles have identical values as £i = £2 = ^3 = ^- Here, we fix the parameters as the same 



values mentioned in Fig. 8.2 



where /^l"^^ is the nonlinearity parameter of the equilateral non-Gaussianity and 



is the transfer 

3 « 



function of scalar mode described in Eq. (3.101 ). Note that these three spectra vanish for X]n=i — 
odd. From this figure, we can estimate the typical amplitude of the III spectra from at large 
scale as 



m\ 



X 3.2 X 10" 



/GeVy /^y 



^.59) 



This equation also seems to be applicable to the III spectra from WW"^. On the other hand, the 
CMB bispectrum generated from the equilateral-type non-Gaussianity on a large scale is evaluated 
with /nl"'' as 



X 4 X 10 



-15 



f equil 

/nl 



300 



(8.60) 



From these estimations and ideal upper bounds on /^l"'' estimated only from the cosmic variance 

0.1, we find a rough limit: 



100 



28 



30 



namely /X'^ ^ 300 and 



7" 



y the signals for 



= 1 



A > 3 X lO^GeV. 

even due to the comparison with the parity-conserving 



for ^ < 

Here, we use on 

bispectrum from scalar-mode non-Gaussianity. Of course, to estimate more precisely, we will have 
to calculate the signal-to-noise ratio with the information of Yln=i = odd 19 



8.3 Summary and discussion 

In this section, we have studied the CMB bispectrum generated from the graviton non-Gaussianity 
induced by the parity-even and parity-odd Weyl cubic terms, namely, and WW'^, which have a 
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dilaton-like coupling depending on the conformal time as / oc r^. Through the calculation based 
on the in-in formalism, we have found that the primordial non-Gaussianities from WW^ can have 
a magnitude comparable to that from even in the exact de Sitter space-time. 

Using the explicit formulae of the primordial bispectrum, we have derived the CMB bispectra 
of the intensity (J) and polarization {E, B) modes. Then, we have confirmed that, owing to the 
difference inj^he transformation under parity, the spectra from vanish in the i space where 
those from WW"^ survive and vice versa. For example, owing to the parity-violating WW"^ term, 
the III spectrum can be produced not only for Yl^=i = even but also for Yl^=i = odd, and 
the IIB spectrum can also be produced for Y^n=i = even. These signals are powerful lines of 
evidence the parity violation in the non-Gaussian level; hence, to reanalyze the observational data 
for Y2n=i = odd is meaningful work. 

When A = —1/2,0,1/2, and 1, we have obtained reasonable numerical results of the CMB 
bispectra from the parity-conserving and the parity-violating WW^. For A = ±1/2, we have 
found that the spectra from and WW"^ have almost the same magnitudes even though these 
have asmall difference in the shapes. In contrast, if A = and 1, we have confirmed that the signals 
from WW^ and vanish, respectively. In the latter case, we will observe only the parity- violating 
signals in the CMB bispectra generated from the Weyl cubic terms. We have also found that the 
shape of the non-Gaussianity from such Weyl cubic terms is quite similar to the equilateral-type 
non-Gaussianity of curvature perturbations. In comparison with the III spectrum generated from 
the equilateral-type non-Gaussianity, we have found that if r = 0.1, A > 3 x lO^GeV corresponds 
approximately to f^^^ < 300. 

Strictly speaking, to obtain the bound on the scale A, we need to calculate the signal-to-noise 
ratio with the information of not only X]n=i ~ even but also X]n=i ^" ~ ^"^"^ f*-*^ each A by the 
application of Ref. |19||. This will be discussed in the future. 
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9 CMB bispectrum generated from primordial magnetic 
fields 

Recent observational consequences have shown the existence of O(10~^)G magnetic fields in galax- 
ies and clusters of galaxies at redshift z ~ 0.7 — 2.0 (l}]3). One of the scenarios to realize this 
is an amplification of the magnetic fields by the galactic dynamo mechanism (e.g. [i]), which re- 
quires C(10"^°)G seed fields prior to the galaxy formation. A variety of studies have suggested 
the possibihty of generating the seed fields at the infiationary epoch [5||6], the cosmic phase tran- 
sitions [7||8], and cosmological recombination [9 11 and also there have been many studies about 
constraints on the strength of primordial magnetic fields (PMFs) through the impact on the cosmic 
microwave background (CMB) anisotropies, in particular, the CMB power spectrum sourced from 
the PMFs 



12 17 



The PMFs excite not only the scalar fiuctuation but also the vector and tensor 
fiuctuations in the CMB fields. For example, the gravitational waves and curvature perturbations, 
which come from the tensor and scalar components of the PMF anisotropic stresses, produce addi- 
tional CMB fiuctuations at large and intermediate scales 15, 17 . In addition, it has been known 



that the magnetic vector mode may dominate the CMB small-scale fiuctuations by the Doppler 
effect (e.g. fll|l5)). 



The PMF anisotropic stresses depend quadratically on the magnetic seed fields. Thus, assuming 
the gaussianity of the PMF, the anisotropic stress and CMB fiuctuation obey the highly non- 
Gaussian statistics 18 , 19 . Owing to the Wick's theorem, the CMB bispectrum contains the pure 
non-Gaussian information. Hence, to extract the information of the PMF from the CMB fields, 
the analysis of the CMB bispectrum is of great utility. Recently, in Refs. 20 23 , the authors 
investigated the contribution of the scalar-mode anisotropic stresses of PMFs to the bispectrum of 
the CMB temperature fiuctuations. From current CMB experimental data, some authors obtained 
rough limits on the PMF strength smoothed on IMpc scale as -BiMpc < 0(l)nG. However, since in 
all these studies, the complicated angular dependence on the wave number vectors are neglected, 
there may exist any uncertainties. In addition, the authors have never considered the dependence 
on the vector- and tensor-mode contributions and hence more precise discussion including these 
concerns should be realized. 

With these motivations, we have studied the CMB scalar, vector and tensor bispectra induced 
from PMFs and firstly succeeded in the exact computation of them with the full-angular depen- 
26 by applying the all-sky formulae for the CMB bispectrum 27 ^ In our studies, we 



24 



dence 

also updated constraints on the PMF strength. 

In this section, after reviewing the impact of PMFs on the CMB anisotropies, we present the 
derivation of the CMB bispectra induced from PMFs and discuss their behaviors. In addition, we 



put limits on the PMFs by considering the WMAP data and the expected PLANCK data 29 , 30 



Finally, we mention our future works. These discussions are based on our studies 24 26 



9.1 CMB fluctuation induced from PMFs 

The PMFs drive the Lorentz force and the anisotropic stress, and change the motion of baryons 
(protons and electrons) and the growth of the gravitational potential via the Euler and Einstein 



equations. Consequently, the photon's anisotropy is also affected. We illustrate this in Fig. 9.1 



In the following discussion, we summarize the impacts of PMFs on the CMB fluctuations in detail 



^In Ref. 28 , after us, the authors presented an analytic formula for the CMB temperature bispectrum generated 
from vector anisotropic stresses of the PMF. 
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Figure 9.1: Interaction between several components in tlie Universe if tlie PMF exists. 



and current constraints on the PMFs obtained from the CMB power spectrum. 



9.1.1 Setting for the PMFs 

Let us consider the stochastic PMFs -B''(x, r) on the homogeneous background Universe which is 
characterized by the Friedmann-Robertson-Walker metric, 



(9.1) 



where r is a conformal time and a(r) is a scale factor. The expansion of the Universe makes the 
amplitude of the magnetic fields decay as 1/a? and hence we can draw off the time dependence as 
-B*(x, r) = i?^(x)/a^. Each component of the energy momentum tensor (EMT) of PMFs is given 

by 



1 



-pB 



Ana* 



= , 



52 ( 



-p,(r)A^(x'^) 



(9.2) 



p,(r) [AB{xnS', + U' 



^ ( 
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The Fourier components of the spatial parts are described as 

^-W-S^/lW^™"-'''- (9.3) 

nL(k) = -J / ^s^(k')i?.(k - k') , 

where we have introduced the photon energy density in order to include the time dependence 
of and p.yfi denotes the present energy density of photons. In the following discussion, for 
simplicity of calculation, we ignore the trivial time-dependence. Hence, the index is lowered by 6bc 
and the summation is implied for repeated indices. 

Assuming that i?''(x) is a Gaussian field, the statistically isotropic power spectrum of PMFs 
Psik) is defined by[^ 

{B,ik)B,ip)) = (27r)3^P.,(k)5(k + p) , (9.4) 

with a projection tensor 

Pafe(k) = J2 e^Hk)e[-")(k) = Sab - kah , (9.5) 

(T=±l 



which comes from the divergence free nature of PMFs. Here k denotes a unit vector and e'^^^ is a 
normalized divergenceless polarization vector which satisfies the orthogonal condition; k"'e^^^^ = 0. 
The details of the relations and conventions of the polarization vector are described in Appendix 



D Although the form of the power spectrum Psik) is strongly dependent on the production 



mechanism, we assume a simple power law shape given by 

PB{k)=ABk^^, (9.6) 

where Ab and hb denote the amplitude and the spectral index of the power spectrum of magnetic 
fields, respectively. In order to parametrize the strength of PMFs, we smooth the magnetic fields 
with a conventional Gaussian filter on a comoving scale r: 

Bl^ r^-^e-'^^'^PBik), (9.7) 



27r2 



then, Ab is calculated as 



where T{x) is the Gamma function and kj. = 27r/r. 

We focus on the scalar, vector and tensor contributions induced from the PMFs, which come 
from the anisotropic stress of the EMT, i.e., IlBab- Following the definition of the projection 



^Here we neglect the helical component. This effect will be considered in Ref. 



31 
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operators in Appendix [Dl tlie PMF anisotropic stress fluctuation is decomposed into 



nS(k) = 






2 


4f Hk) = 


2 



(9.9) 



Tliese act as sources of tlie CMB scalar-, vector- and tensor-mode fluctuations as follow. 
9.1.2 Scalar and tensor modes 

If the seed magnetic fields exist in the early Universe, the scalar and tensor components of the 



PMF anisotropic stress are not compensated prior to neutrino decoupling 15, 17 , and the scalar 
and tensor metric perturbations generated from them survive passively. These residual metric 
perturbations generate the CMB anisotropics of the scalar and tensor modes. These kind of CMB 
anisotropics are so called "the passive mode" and may dominate at intermediate and large scales 



depending on the PMF strength 17 



To estimate the curvature perturbations generated from PMFs, let us consider the metric 
perturbation 



ds'^ = a^[-{l + 2A)dT'^ - 2CidTdx' + {6ij + 2Hij)dx'dx^] 
Then, the Gauge-invariant Bardeen potentials [32] are 



(9.10) 



3 ^ 



A;2 

' k fc2 



(9.11) 



where we obey the vector and tensor decomposition of Appendix |D] as 
a(k) = C7(°)(k)Or+5^CW(k)OW, 

A=±l 

^a.(k) = -\h-\\.)6^, + ff(°)(k)Oi;) + /f(^)(k)Oi^) + ^ 

A=±l A=±2 

The energy momentum tensor is given by 

T% = -p(l + A) , 



(A) 
ab 



(9.12) 



P 



P 



w 



r + 5'i + n 



(9.13) 



where w = p/p is the coefficient of the equation of state and = p/p is the sound speed. The 
anisotropic stress is also decomposed into 



n„,(k) = ni°)(k)oi? + Y: ni^Hk)oit) + Y ni'^(k)oiJ) . 



A=±l 



A=±2 



W 



Before neutrino decoupling, the Universe is dominated by the radiative fluid with 
The fluid is tightly coupled to the trace amount of baryons and can not create any anisotropic 



(9.14) 



1/3. 
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stress. Hence, in this period, the total anisotropic stress comes from only PMFs, namely, the 



constant H^]. Until neutrino decoupling, there is no mechanism to compensate Il^j, and it will 
be a source of the gravitational potential via the Einstein equation. Using % = and taking 
into account the radiation-dominated limit, the Einstein equation is written as 



(0) 



3{kT) 



(9.15) 



where ~ 0.6 is the ratio by the energy density of photons to all relativistic particles and we 
assume F = 0. This can be solved exactly and the superhorizon limit, when /cr <^ 1, is 



Cl 



— -c.-^i?,nSlnr 



{kry (kr)^ Qkr 
We have the interest in the contribution on the comoving curvature perturbation 



(9.16) 



C = $ + 



3(l + w) 



n 



(9.17) 



From the above estimation, this superhorizon limit becomes 



C(r) = ars) - i?,ng] 



yrBj 2r 2 



(9.18) 



where tb is the conformal times at the generation of the PMF. At = IMeV^ , neutrinos 
decouple from the radiative fluid. By considering the Boltzmann equations of neutrinos, we find 
the solution for r ^ Tj, as fli"'' = — (i?^/i?j,)n^] with Ri, ~ 0.4 being the ratio by the energy 
density of neutrinos to all relativistic particles. This implies that the neutrino anisotropic stress 
compensates the magnetic anisotropic stress. Hence, when this compensation is effective, the 
source term of the Einstein equation become zero and the gravitational potentials stop growing. 
Including the further bare growth, the total curvature perturbation in the superhorizon limit for 
r ^ r,^ is evaluated as 



C(k)^C(k,rB) 



R H^°^ 



(k) In 



Tb 



(9.19) 



Note that this logarithmic growth and saturation of the curvature perturbation even on superhori- 
zon scales is caused by only the property of the sources as p oc a~^, hence the above discussion is 
valid for the general radiation fluid other than PMFs 33 . This curvature perturbation acts as a 
source of the CMB fluctuation of the scalar mode. 

Similar discussion can be given in the tensor perturbation. Here, we shall consider the spatial 
metric and decomposition into the helicity states as |^ 



9ij 



eif)(k)/i„,(k,r) 



(9.20) 



As is well known, the gravitational potential of tensor modes can be generated from anisotropic 
stresses via the Einstein equation. Like the scalar case, after neutrino decoupling, the anisotropic 



^On superhorizon scales, this C is equal to 
4/i(±2) is equal to -VSHt of Refs. [l5 



17 



-7^ in Eq. (2.40). 
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10 (r ; 



stresses of PMFs vanish via the compensation of those of neutrinos. However, prior to this epoch, 
there is no compensation process due to the absence of the neutrino anisotropic stresses. Hence, 
from the Einstein equation, we find the evolution equation of the tensor-mode metric perturbations 
as 

(±2) < < 

,~ '''' , (9.21) 

From this equation, we can find a superhorizon solution of the tensor metric perturbation for 

r > Tj, as 

/i(±2)(k) ^ /i(±2)(k^ + 6i?^Hg'^(k) In . (9.22) 

This gravitational wave also generates additional CMB fluctuation. 

Note that the effects of PMFs on the transfer functions 7J/ and 7J/ are inconsiderable 



at larger scales 17, hence it is safe to use the non-magnetic transfer functions (3.101) in the 
computation of the CMB spectra induced from scalar and tensor modes of the PMF anisotropic 
stress. 

9.1.3 Vector mode 

The vector mode has no equivalent passive mode as the gravitational potential of the vector mode 
decays away via the Einstein equation posterior to neutrino decoupling. Thus, in the vector mode, 
we need to consider the impact of the PMF on the transfer function. In Refs. [T3| [M], [25| |28] , it is 
discussed that the temperature fluctuations are generated via Doppler and integrated Sachs- Wolfe 
effects on the CMB vector modes. On the basis of them, we derive the transfer function of the 
vector magnetic mode as follows. 

When we decompose the metric perturbations into vector components as 

Sgoc = Sqco = a^Ac , 

59., = {d^^P + d,hp) , ^^-^^^ 

we can construct two gauge-invariant variables, namely a vector perturbation of the extrinsic 
curvature and a vorticity, as 

V = A-h , , , 

ri = V - A , 

where v is the spatial part of the four-velocity perturbation of a stationary fluid element and a 
dash denotes a partial derivative of the conformal time r. Here, choosing a gauge as h = 0, we 
can express the Einstein equation 

V + 2«V = -'°-^^-°'"^'^' + "-'^' + "i''^'>. (9.25) 

and the Euler equations for photons and baryons 

fl^ + k{^r^ - Vb) = , (9.26) 
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Ob + nQh 



K 



L(v) 



Vb 



Here L*^"^) = kp^^oU^^ is the Lorentz force of vector mode and 11^ = —ikbPaJlbc, P is the isotropic 
pressure, the indices 7, z/ and b denote the photon, neutrino and baryon, n is the optical depth, and 
R = (pf, + pb)/{p'y + p-y) — 3pb/ (4p^). In the tight-coupling limit as ~ Vb, the photon vorticity 
is comparable to the baryon one: fl^ ~ fib = O. Then, the Euler equations (9.26) and (9.27) are 
combined into 

1 + R)n + Run = , ^ , , (9.28) 



a^{Pb + Pb) 



(9.27) 



and this solution is given by 



0(k,r) 



«^(P7 

rL(^)(k) 

^1 + R){p-yfl+P-yfl) 



(9.29) 



Note that Eq. (9.28) and the above solution are valid for perturbation wavelengths larger than 
the comoving Silk damping scale Ls = 2ii/ks, namely, k < ks, where photon viscosity can be 
neglected compared to the Lorentz force. For k > ks, due to the effect of the photon vorticity, the 



Euler equation (9.28) is changed as 12 



{1 + R)n + R{'H + 



apb 



L(V) 



a^(P7 +P7) 



(9.30) 



where x = (4/15)p^L^a is the photon shear viscosity coefficient and = k ^ is the photon 
comoving mean-free path. We can obtain this analytical solution: 



0(k,r) 



L(v)(k) 



(PL^/5)(p^,o +P7,o) 
Hence, we can summarize the vorticity of the baryon and photon fluids as 

0(k,r)^/3(fc,r)ng')(k) , 



(9.31) 



/3(fc,r) 



P7,0 



P7,0 + ^7,0 



X 



kT/{l + R) for k<ks 
hkjk for k > ks 



(9.32) 



As mentioned above, the CMB temperature anisotropics of vector modes are produced through 
the Doppler and integrated Sachs- Wolfe effects as 



A/(n) 



n 



\T0 



drV-n 



(9.33) 



where tq is today and t^, is the recombination epoch in conformal time, pk,n = k ■ n, x = A;(ro — r), 
and n is a unit vector along the line-of-sight direction. Because of compensation of the anisotropic 
stresses, a solution of the Einstein equation (9.25) expresses the decaying signature as V oc 
after neutrino decoupling. Therefore, in an integrated Sachs- Wolfe effect term, the contribution 
around the recombination epoch is dominant. Furthermore, neglecting dipole contribution due to 
V today, we can form the coefficient of anisotropics as 



''Ilm 



n 
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(2^ 



d^n[ng^)(k)-n]F4(n)/3(fc,r,)e- 



(9.34) 



In tlie transformation n — )■ {fik,n, (pk,n), tlie functions are rewritten as 



A=±l 



F;„(n) ^ (^(k)) Y,un,,, 

m! 

(fh — )■ dflk n , 



(9.35) 



where we use the relation: ni^^ = Ea=±i -^nj^^jei^^ and the Wi mer D matrix under the rotational 
transformation of a unit vector parallel to z axis into k corresponding to Eq. (3.96). Therefore, 



performing the integration over flk,n in the same manner as Sec. [3| we can obtain the explicit form 
of cij £m Slid express the radiation transfer function introduced in Eq. (|3.100|) as 



4Z = M-') 



(27r)3 



J2 A-Ar4(k)ng](k)r,;P(A:) 



A=±l 



{e-iy. 



.1/2 



(9.36) 



This is consistent with the results presented in Refs. 15 , 34 . In the same manner, the vector-mode 



transfer functions of polarizations are derived 14 , 19 . Then, we can also express as 

d^k 



^ ' A=±l 



(9.37) 



9.1.4 Expression of a^m's 



As the above results and Eq. (5.3), the CMB intensity and polarization fluctuations induced from 
PMFs are expressed as 



k^dk 
(27r)3 



dW.Uk) 



R-y In 



d^k^ir4(k)ng;'Hk) 



nS(k) 



d%,Y;^{k) 



6Rj In 



4?Hk) 



(9.38) 



where we take C(k, r^) = /i(=^^)(k, r^) = 0, and regard and ^^^'^^ as — C(= Tl) and h^^'^\ 
respectively 



In Refs. mlp35^, we equate to C- 
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9.1.5 CMB power spectrum from PMFs 



From the formulae (9.38), the CMB power spectra from PMFs are written as 

2 



n 

\n=l 



n=l 



Kdkn^iz„) 



(2 



71 



'7'ifA(^n)5^[sgn(A„)]^"+- 



. (9.39) 



\n=l 



To compute the initial power spectrum {^nn=i v„m„(^n)/' "^^ need to deal with the power spectrum 
of the anisotropic stresses as 



\-2 



'7,0 J 



n 

n=l 



(2^ 



x^[P,,(k\)n,(kT^'i) + P„c(k'i)n,(kT^'i)] . 



(9.40) 



Note that this equation includes the convolution integral and the complicated angular dependence. 
In Refs. 



16, 36 38 



the authors performed the numerical and analytical computation of this 
convolution integral over k'j^ and provided the fitting formulae with respect to the magnitude of 
the wave numbers ki for each value of the magnetic spectral index hb- 



In Fig. 9.2, we plot the power spectra of the intensity anisotropics (9.39) for the scalar, vector 
and tensor modes when magnetic spectrum is nearly scale invariant as = —2.9. Here 



we 



assume that the PMFs generate from the epoch of the grand unification = 10 GeV to that 
of the electroweak phase transition r^^ = lO^GeV. Firstly, we will see that the shapes of the 
tensor and scalar power spectra are similar to those of the non-magnetic case coming from the 
scale- invariant primordial spectra shown in Fig. |3.5 This is because PMFs impact on only the 
primordial gravitational waves and primordial curvature perturbations, and do not change the 
transfer functions of the tensor and scalar modes. For i < 100, the tensor mode dominates over the 
intensity signal. The scalar mode seems to dominate in the intermediate scale as 100 ^ i ^ 2000. 
The vector-mode spectrum monotonically increases for larger than Silk damping scale, namely, 
£ < ksTQ ~ 2000, and decreases for i > 2000. These features seem to trace the scaling relation 



of the transfer function in terms of the wave number (9.32). Hence, we can understand that the 



latter damping effect arises from the viscosity of photons. The vector mode seems to show up for 
very small scale, namely, i > 2000. 

In this figure, for comparison, we also plot the CMB intensity power spectrum from the primor- 
dial curvature perturbations not depending on the PMF. In principle, comparing this spectrum 
with that sourced from PMFs leads to bounds on the PMF parameters. Actually, the researchers 
perform the parameter estimation by the Malkov Chain Monte Carlo approach 36, 37 39 43 . So 



far, the most stringent limit on the PMF strength from the CMB two-point function data of the 
intensity and polarizations are -BiMpc < 5nG and ub < —0.12 |41j. In Refs. 36, 43 , combining 
the CMB data with the information of the matter power spectrum, tighter bounds are gained. 

As discussed above, conventionally, the CMB power spectra from PMFs are computed by using 
the fitting formulae for the power spectra of the magnetic anisotropic stresses. However, without 
these formulae, we can obtain the CMB power spectra by applying the mathematical tools such 
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Figure 9.2: Power spectra of the CMB intensity fluctuations. Tlie red solid, green dashed and blue 
dotted lines correspond to the spectra generated from the tensor, vector and scalar components of 
the PMF anisotropic stress for ub = —2.9, respectively. The upper (lower) line of the red solid and 
blue dotted ones are calculated when t^/tb = 10^^(10^). The magenta dot-dashed line expresses 
the spectrum sourced from the primordial curvature perturbations. The strength of PMFs is fixed 
to -BiMpc = 4.7nG and the other cosmological parameters are fixed to the mean values limited from 



WMAP-7yr data reported in Ref. 29 



as the Wigner symbols 25 . In the remaining part, focusing on the vector mode, we present this 



new approach and show the consistency with the conventional result. Then, we use some colors in 
the equations for readers to follow the complex equations more easily. 



From Eq. (9.39), the CMB power spectrum of the intensity mode induced from the magnetic- 



vector-mode anisotropic stress is formulated as 

2 



(.V) (V)* 
^I,eimi^I,e2m2 



n=l 



tv^atvn ,-i-{V), 



I I In V « 



X 



(2vr: 



Ai,A2=±l 

'n,ii_^h/2^mi,m2 ; (9.41) 

where the initial power spectrum, which is expanded by the spin spherical harmonics, is 



n 



(Ai) 

Bv,£imi 



{A2)* (I. ^ 



-4vrp^,o) ^ / ki / k2_Ai>"£tmi(kl)-A2>"£2m2(k2) 



X 



k'^dk'PB{k[) / k'^dk'^PB{k'^) 
Jo 

X / d^k; / rf2kV(ki - k; - k'2)5(k2 - k'2 - k'l) 
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(9.42) 



where we use Eq. (9.9) and the definition of the vector projection operator, O^^^^(k), in Ap- 
pendix D Note that we rewrote the power spectrum (9.40) as more symmetric form in terms of 
ki,k2 and ka. Then, we should simplify this initial power spectrum. 



For the first part in two permutations, we calculate (5-functions and the summations with 
respect to a, b, c and d: 



M1M2M3 

l' +3L' +3L' 

B'dB {-^)^^^ll'^U3L'ShB)3vSkWL'SKB) 

Tl T' T' 



X y?.M;(kJni„aki)F4_„;(ki)(-l)"; ( \ ^, ) ,(9.43) 

a=±l ma,m^=±lfi ^ 



'47r\ 



k,J-'^\k,)PUK) = E E ( 

o-'=±l mc,mi,=±l,0 



X Yi,m^ (kl ) A2 Ylma ) -a'^am, (^l ) a^lm^ (k'l ) 

1 (-Ai) 



47r\^ 



3 / 



X Fi^^(k2)_,,n^,(k;)_<,,r;^^(k^).,y* (k^) , 
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perform the angular integrals of the spin spherical harmonics: 



= E(-i)" 



a+nia tO—ct—S tO—u—S 

LMS 



1 L 



L2 

M2 



1 L 

-rua M 



J d'^^'^-a'yimJ^L'^M^a'y'lnnXLs-Ms 

E/i \a'+mc jO-a'-S' jO-a'-S' 

L'M'S' 



L3 ^ L' \f L'2 
-M3 rub M' ) \ M'^ 



1 



-rUc M' 



^Li^iLfc ^llLfc 1 Mfc i V 



1 1 Lk 
ma nib Mk 



1712 Mp 



E/_-| \m2+A2 rOA2-Sp / 

Lp Mp Sp 

sum up the Wigner-3j symbols over the azimuthal quantum numbers: 



-vfic -rud Mp 



M1M2M3 



X 



■a rrib Ml 



L3 1 L' 
-Ms nib M' 



^ (_-]\M+ii+L3+L+l ( L' L 

^ ' \M' -M ni] 



Li L2 \ 

Ml -M2 -M3 j 

1 L \ / Li £1 Lfe \ 
, -ma M J \ Mi mi Mk J 

V L L \ { ^' ^ 1 



L2 

M2 - 



-nic -md Mp 
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-2 1 L' 
M^ -mc M' 



1 L 



-M's md M 
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L'l £2 Lp 
-M[ m2 Mp 

^ (_-l\M'+e2+L'2+L+l+Lp f L' L 

^ ' \M' -M J 1 ^ 

sum up the Wigner-3j symbols over M and M': 

\-(_^^M+M' ( L' L ii \f L' L £2 \_ (-ir% . 

^ V M' -M mi [ M' -M ni2 / ~ 2i, + 1^''^' ' 

MM' ^ ^ ^ /II 




(9.44) 



(9.45) 



Following the same procedures in the other permutation and calculating the summation 



(9.46) 



over Lr, 
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E7-OA2-A2 T-OAa-Aai+ljO^ j r, r, r] 1 

-'LifaLp ^lli-P 2 I ^ ^ ^ I 

Lp [ 1 1 Lp J 

we can obtain the exact solution of Eq. (9.42) as 



OA2-A2 



-' T.I 




(9.47) 
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^L'^IL ^L2lL -'L3IL' ^L'^IL' 



\^ \ \ 7-OA1-A1 7-OAi-Ai 7-OA2-A2 



(9.4J 



Note that in this equation, the dependence on the azimuthal quantum number is included only in 
^mi,m2- In the similar discussion of the CMB bispectrum, this implies that the CMB vector-mode 
power spectrum generated from the magnetized anisotropic stresses is rotationally-invariant if the 
PMFs satisfy the statistical isotropy as Eq. (9.4). 
Furthermore, using such evaluations as 



E( _'\ \L2+L'+Ls+L' 7-05-5 7-05-5 7-05'-5' 7-05'-5' 
I ^) -'l^jIL ^L2lL -'L3IL' ^L'^IL' 
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(9.49) 
(9.50) 
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(9.51) 



the CMB angle-averaged power spectrum is formulated as 



c 



{V) 

III 



2ti /32(27r) 



/(2£+l) 



47r 



k'^dk 
(27r)3 
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L'Lg ' 
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3 ^, + ^; 
1 2 



-L' +L3 7-0 7-000 7-01-1 7-01-1 7-01-1 7-01-1 




(9.52) 



Here, we use the thin LSS approximation described in Sec. |9.3.1[ This has nonzero value in the 
configurations: 



),(1,^ 



L[ = \i±2\ 



(Lfc,L0 = (2,|^±2|),(2, 
\i- L\ < L' < i + L , 

{L2, L'3) = (|L - 1|, \L ± 1|), (L, L), (L + 1, |L ± 1|) , 
(4, L3) = (|L' - 1|, |L' ± 1|), {L', L'), {L' + 1, \L' ± 1|) , 
Li + L2 + -L3 = even , L[ + L'2 + L'^ = even , 
\Li — L2I <^ Li -\- L2 , — < L3 < L\ + L2 . 



(9.53) 



This shape is described in Fig. 9.3 From this figure, we confirm that the amplitude and the 



overall behavior of the red solid line are in broad agreement with the green dashed line of Fig. 9.2 



and the previous studies (e.g. 15-17, 44]). For ^ < 2000, using the scaling relations of the Wigner 
symbols at the dominant configuration L ~ L' ~ 1 as discussed in Sec. 9.5, we analytically find 
that Cf2 oc r^+l This traces our numerical results as shown by the green dashed line. 

This computation approach is of great utility in the higher-order correlations. In the next 
subsection, in accordance with this approach, we compute the CMB bispectra sourced from PMFs. 



9.2 Formulation for the CMB bispectrum induced from PMFs 

In this subsection, we derive the explicit form of the CMB bispectra induced from PMFs by 
calculating the full-angular dependence which has never been considered in the previous studies 



20 23, 28 . The following procedures are based on the calculation rules discussed in Ref. 25 . Note 



that we use some colors in the following equations for readers to follow the complex equations more 
easily. 



9.2.1 Bispectrum of the anisotropic stress fluctuations 



According to Eq. (9.2), EMT of PMF at an arbitrary point, T^„{'k), depends quadratically on the 



Gaussian magnetic fields at that point. This non-Gaussian structure is identical to the local- type 
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Figure 9.3: CMB power spectra of the temperature fluctuations. The lines correspond to the spec- 



tra generated from vector anisotropic stress of PMFs as Eq. (9.52) (red solid line) and primordial 
curvature perturbations (blue dotted line). The green dashed line expresses the asymptotic power 
of the red solid one. The PMF parameters are fixed to ns = —2.9 and -BiMpc = 4.7nG, and the 
other cosmological parameters are fixed to the mean values limited from WMAP-7yr data reported 
in Ref . 



29 



non-Gaussianity of the curvature perturbations as mentioned in Sec 4^, hence it is expected that 
the statistical properties of the magnetic fields obey those of the local-type non-Gaussianity. This 



will be automatically shown in Sec. 9.4 



Using Eq. (9.4) and the Wick's theorem, the bispectrum of the anisotropic stresses is calculated 



as 



{UBab{kl)UBcd{k2)'nBef{k3)) 



n 

n=l 



(27r)= 



X (i?„(k;)i?,(ki - k\)B,{k'^)B,{k2 - k'^)B,{k'^)Bf{k3 - k'3)) 



n=l 



X5(ki - k'l + k'3)5(k2 - k'2 + k'i)<5(k3 - k'3 + k'2) 

xl[Pad{ki)Pbe{K)Pcf{k2) + {a ^ b 01 C ^ d OT 6 ^ f}] , 
o 

(9.54) 

as k~p^ ~ C(0.1)Mpc and the curly 



45 



46 



where ka is the Alfven-wave damping length scale 
bracket denotes the symmetric 7 terms under the permutations of indices: a -H- 6, c -H- ci, or e -H- /. 
Note that we express in a more symmetric form than that of Ref. [Ts] to perform the angular 



integrals which is described in Sec. 9.2 To avoid the divergence of {IlBab{ki)IlBcdik2)IlBefiks)) 
in the IR limit, the value range of the spectral index is limited as ub > —3. We note that this 
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bispectrum depends on the Gaussian PMFs to six, hence this is highly non-Gaussian compared 
with the bispectrum of primordial curvature perturbations proportional to the Gaussian variable 
to four as shown in Sec. 14.21 



9.2.2 CMB all-mode bispectrum 



Following the general formula (5.5) and using Eq. (9.38), the CMB bispectra induced from PMF 
are written as 



Ha 

\,n=l 



fj47r( 



n=l 



(27r) 



An 



\n=l 



(9.55) 



Remember that the index Z denotes the mode of perturbations: Z = S (scalar), = V (vector) or 
= T (tensor) and its helicity is expressed by A; A = for {Z = S), = ±1 for {Z = V) or = ±2 for 
{Z = T), X discriminates between intensity and two polarization (electric and magnetic) modes, 
respectively, as X = I,E,B and x is determined by it: x = for X = I,E or = 1 for X = B. In 
the following discussion, we calculate (nn=i ^i^m (^n) / ^^"^ fi^'i explicit formulae of the CMB 



bispectra corresponding to an arbitrary Z. 



Using Eqs. (9.38) and (9.54), we can write 



\n=l 



>-3 



n=l 



n / ci'k„_,„r,:„„(k„) / k'^dk'MK 



ku 



x5(ki - k'l + k'3)5(k2 - k; + k;)5(k3 - k'3 + k'2) 

xc^-A,oi,^^)(k;)c:,^oL-'^^(k;)ci,3o(/^)(k3) 

xPad{K)Phe{K)PcA^2) , 



(9.56) 



with 

[1^7 In (^) (A = 0) 

C'x^ll (A = ±1) . (9.57) 

[3i?,ln(^) (A = ±2) 

Let us consider this exact expression by expanding all the angular dependencies with the spin- 
weighted spherical harmonics and rewriting the angular integrals with the summations in terms of 
the multipoles and azimuthal quantum numbers. 

In the first step, in order to perform all angular integrals, we expand each function of the wave 
number vectors with the spin- weighted spherical harmonics. By this concept, three delta functions 
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are rewritten as 



POO 

<5(ki - k; + k'3) = 8 / A'dA i-^)'"''^lllLnAhA)j,MA)jL^^^^^ 

■^0 L1L2L3 
M1M2M3 
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^0 T I T I T ! 
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T n T n T II 



L-^ 7^2 -^3 



// I 5 



where 



(2('i + l)(24 + l)(2('3 + l) / <, <2 <3 



The other functions in Eq. (9.56), which depend on the angles of the wave number vectors, can be 
also expanded in terms of the spin-weighted spherical harmonics as 

{-' 



Oi;'^^(kl)OL-'^^(k;)0(7^^)(k3)Pa.(k\)ne(k;,)Pc/(k\) 



3 „ Un^ 2 
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2 



5 



(9.60) 

where we have used Eqs. (9.5), some relations in Appendices [C| and [Pj and the definition (D.23) 



as 



Almairit 

-2 (A = 0) (9-61) 
Cx = { 2VSX (A = ±1) . 
2V3 (A = ±2) 
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In the second step, let us consider performing all angular integrals and replacing them with the 
Wigner-3j symbols. Three angular integrals with respect to k'-L,k2 and kg are given as 



LM S=±l ^ 
(T\rncT0-a2-S'T0-a2-S' 
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rrid M 
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L" 



-me M" 



(9.62) 



where we have used a property of spin- weighted spherical harmonics given by Eq. (C.7). We can 
also perform the angular integrals with respect to ki, k2 and kg as 



j2u V* V* V* rOAi-Ai 
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(9.63) 



At this point, all the angular integrals in Eq. (9.56) have been reduced into the Wigner-3j symbols. 



As the third step, we consider summing up the Wigner-3j symbols in terms of the azimuthal 
quantum numbers and replacing them with the Wigner-6j and 9j symbols, which denote Clebsch- 
Gordan coefficients between two other eigenstates coupled to three and four individual momenta 



27, 47-49 . Using these properties, we can express the summation of five Wigner-3j symbols with 
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a Wigner-9j symbol: 
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(9.64) 
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Furthermore, we can also sum up the renewed Wigner-3j symbols arising in the above equations 
over M, M' and M" with the Wigner-6j symbol as [50] 
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(9.65) 



With this prescription, one can find that the three azimuthal numbers are confined only in the 

il i2 i3 

mi m2 m^ 

ensures the rotational invariance of the CMB bispectrum as pointed out above. 

Consequently, we can obtain an exact form of the primordial angular bispectrum given by 



Wigner-3j symbol as 



^ . This 3j symbol arises from (^nn=i ^etmA^ri)^ and exactly 
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(9.66) 



where 
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Here, the coefficients have been calculated as 
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(9.67) 
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(9.68) 



Substituting Eq. (9.66) into Eq. (9.55), we can formulate the CMB bispectra generated from 
arbitrary three modes such as the scalar-scalar-vector and tensor-tensor-tensor correlations with 
the / function as 
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(9.69) 



From this form, it can be easily seen that due to the sextuplicate dependence on the Gaussian 
PMFs, the Wigner-6j symbol connects the true multipoles (£1,^2 and £3) and the dummy ones 
(L, L' and L"), and the 1-loop calculation with respect to these mulipoles is realized as illustrated 
in the left panel of Fig. |9.4[ Due to the extra summations over L, L' and L", it takes a lot of time 
to compute this compared with the tree-level calculation presented in the previous sections. 



9.3 Treatment for numerical computation 

In order to perform the numerical computation of the CMB bispectra, we give the explicit angle- 



averaged form of Eq. (|9.69|) as 

,(^1^223) 



B 



Cz^Cz2Czz (-47rp^,o) ^ W^/ ^ } 

LL'L" ^ ^ 
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Figure 9.4: Diagrams with respect to multipoles 35 . The left panel corresponds to Eq. (9.69). 
Due to the Wigner-6j symbol originated with the sextuplicate dependence on the Gaussian PMFs, 
the true multipoles £1,^2 and £3 are linked with the dummy ones L,L', and L" and the 1-loop 
structure is realized. The right panel represents the tree-structure diagram, which arises from the 
CMB bispectrum induced by the four-point function of the Gaussian fields as mentioned in the 
previous sections. 
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(9.71) 
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We consider performing the summations with respect to the hehcities. By considering the selection 
rules of the Wigner-3j symbol, the summations over S, S' and S" are performed as 

E/ 1 \L2+L'+L''+La+L'+L'' tOS~S tOS-S tOS'-S' tOS'-S' tOS"-S" tOS"-S" 
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(9.72) 



By the same token, the summations over Ai, A2 and A3 are given by 
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(9.73) 



where Ns is the number of the scalar modes constituting the CMB bispectrum]^ Thus, we rewrite 
the bispectrum as 
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where we introduce the filter functions as 

Ql'^,L2,L = {h'^,L+l + h'^,\L-l\){h2,L+l + h2,\L-l\) + h'^,Lh2,L 
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(9.74) 



(9.75) 



^Caution about a fact that |A| is determined by Z, namely, |A| = 0, 1, 2 for Z = S, V,T, respectively. 
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The above analytic expression seems to be quite useful to calculate the CMB bispectrum induced 
from PMFs with the full-angular dependence. However, it is still too hard to calculate numerically, 
because the full expression of the bispectrum has six integrals. In addition, when we calculate the 
spectra for large £'s, this situation becomes worse since we spend a lot of time calculating the 
Wigner symbols for large £'s. The CMB signals of the vector mode appear at ^ > 2000, hence we 
need the reasonable approximation in calculation of the CMB bispectra composed of the vector 
modes. In what follows, we introduce an approximation, the so-called thin last scattering surface 
(LSS) approximation to reduce the integrals. 



9.3.1 Thin LSS approximation 

Let us consider the part s of the integrals with respect to A, B, C, k',p' and q' in the full expression 
of the bispectrum (9.74) of Bfj^^j^^^. In the computation of the CMB bispectrum, the integral in 

terms of k, {p and q) appears in the form as J k'^dkTj\^{k)jL^{kA). We find that this integral is 
sharply-peaked at A ~ tq — r,,,, where tq is the present conformal time and is the conformal time 



of the recombination epoch. According to Sec. [9.1.3 the vorticity of subhorizon scale sourced by 
magnetic fields around the recombination epoch mostly contributes to generate the CMB vector 
perturbation. On the other hand, since the vector mode in the metric decays after neutrino de- 
coupling, the integrated Sachs- Wolfe effect after recombination is not observable. Such a behavior 
of the transfer function would be understood on the basis of the calculation in Sec. 19.11 and we 
expect TjJ^{k) oc ji^{k{TQ — r*)), and the fc-integral behaves like 6{A — {tq — t^)). By the numerical 
computation, we found that 
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r*)), (9.76) 



is a good approximation for Li = £i ± 2,ii as described in Fig. |9.5[ Note that only the cases 
Li = £i±2, ii should be considered due to the selection rules for Wigner-3j symbols as we shall see 



later. From this figure, we can find that the approximation (the right-handed term of Eq. (9.76)) 
has less than 20% uncertainty for £i ~ Li > 100, and therefore this approximation leads to only 
less than 10% uncertainty in the bound on the strength of PMFs if we place the constraint from 
the bispectrum data at £1,^2, ^3 ^ 100 Using this approximation, namely A = B = C^to — 
and / dA = f dB = f dC — t- r*/5, the integrals with respect to A, B, C, k',p' and q' are estimated 
as 
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(9.77) 



^Of course, if we calculate the bispectrum at snialler multipoles and the CMB bispectra are produced by other 
modes than the vector one, we may perform the full integration without this approximation. 



9.3 Treatment for numerical computation 



141 





Figure 9.5: The ratio of the left-hand side (exact solution) to the right-hand side (approximate 
solution) in Eq. (9.76). The lines correspond to the case for Li = ii + 2 (red sohd line), for Li = ii 
(green dashed one), and for Li = f i — 2 (blue dotted one). 



Here the function JCjJ, is defined as 

POO 

Jo 



n y^-i r(iv)r(^±^^) 

2y 2^ p^ ^-/'+i+A^ ^p^ -/+/'+i+Af 



(for y,N,l + l' + 2- N >0) 



N -rr i-i'+i+N \-r^ -i+i'+i+N \-rr i+i'+2+N -'^ 

(9.78) 

which behaves asymptotically as lCfii{y) oc for / ~ /' ^ 1. Here we have evaluated the k' 
integrals by setting ko — ^ oo. This is also a good approximation because the integrands are 
suppressed enough for k\p\q' < kn ^ C(10)Mpc~^. 



9.3.2 Selection rules of the Wigner-3j symbol 

From the selection rules of the Wigner symbols as described in Appendix [C} we can further limit 
the summation range of the multipoles as 

\L-i2\<L' <L + i2 , Max[|L -ii\,\L' - 41] < L" < Mm[L + 4, L' + 4] , 
Li + L2 + L3 = even , L\ + -|- L3 = even , -|- + L3 = even , (9.79) 
l-Li — L2I < < i^i + ) l-^'i ~ -^2! — -^3 — -^'1 ~^ -^2 5 1-^1 ~ -^2! — -^3 — -^1 ~^ -^2 5 
and from the above restrictions the multipoles in the bispectrum, £1,^2 and £3, are also limited as 

14-41 <4 < 4 + 4 • (9.80) 

Therefore, these selection rules significantly reduce the number of calculation. In these ranges, 
while L' and L" are limited by L, only L has no upper bound. However, we can show that the 
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summation of L is suppressed at ~ ^2 ~ ^3 ^ -Z^ as follows. When the summations with respect 
to L,L' and L" are evaluated at large L,L' and L", namely ii,i2,^3 <^ L V L",L2 ~ L'^ 
L, L2 ~ Ln ~ L' and L2 L", we get 



3 ~ 



LL'L" ^ ^ LiL'L''"'^ 



L'^L'^Ls 

k'idk'2PBiK)jL'^ ik',B)jL'.l ik',C) I k'idk'^PB{k'^)jL'i {KC)JL, {k'^A) 

Jo 

' 1 \Ylt=i 9 ' rO rO rO rOl-l rOl-l rOl-l rOl-1 rOl-l rOl-l 
^ ^1 ^ 111213-^ L'^L'2L'^-^L'IL'^L'^'^L';^lL'^L2lL-^LpL''^L'2lL''^L3lL'''^L'^lL'' 

K h 
I/3 1/2 I/i 

1 1 2 

oc {LL' L"Y^+^/'' . (9.81) 

LL'L" 

Therefore, we may obtain a stable result with the summations over a limited number of L when 
we consider the magnetic power spectrum is as red as ~ —2.9, because the summations of L' 
and L" are limited by L. Here, we use the analytic formulas of the I symbols which are given by 




5 } ^ (LL'^'T'^" ' ^^Sr'^ « L--^' , I 1 t t I cx (L"L)-V2 , (9.82) 

as described in detail in Appendix [C| 

Using the thin LSS approximation and the summation rules described above, we can perform 
the computation of the CMB bispectrum containing full-angular dependence in a reasonable time. 



9.4 Shape of the non-Gaussianity 

In this subsection, in order to understand the shape of the non-Gaussianities arising from PMFs, 
we reduce the bispectra of the PMF anisotropic stress by the pole approximation [35]. 

Let us focus on the structure of the bispectrum of the PMF anisotropic stresses (9.54). If the 
magnetic spectrum is enough red as ra^ ~ —3, the integral over the wave number vectors is almost 
determined by the behavior of the integrand around at three poles as k[,k2,k'^ ~ 0. Considering 
the effects around at these poles, we can express the bispectrum of the PMF anisotropic stresses 
approximately as 



(nBab(ki)nBcd(k2)nije/(k3)) ~ {-Arcpj^o 

1 



,-3 



aA 



,n=l 



X ■ 



ub + 3 3 

PB{kl)PB{k2)SadPl>e{^l)Pcf{\^2) + Pb {k2) Pb (k^) Pad{\^2) Pbei^s) 6 cf 
+PBiki)PB{k3)Pad(kl)SbePcfO<-3) + {a <-> b OT C ^ d OT 6 ^ f} , 



where we evaluate as 



d-^k'PB{k')Pab{k') ~ a / k"dk'PB{k') / d'k'Pabik' 



aA 



+ 3 



<87r 

y 



(9.83) 



(9.84) 
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Note that a is an un know n parameter and should be determined by the comparison with the exact 
bispectra [(9.69) or (9.74)], and we take k^: = lOMpc"^. 

Under this approximation, the angular bispectrum of the primordial tensor perturbations 
(Ai, A2, As = ±2) is given by 



\n=l 



n=l '' 



R^lYljTy / Tb) 

4vrp^,o 



vn=l 



+PB{k,)PBik3)Pad{kl)5bePcfik3 

(- 
ef 



x(-27)el;^^)(k;)ei-^^)(k;)ei7^=')(k3). 



(9.85) 



Using Eq. (D.23), we reduce the contraction of the subscripts in the TTT spectrum to 

0[-,''\\^l)0i-/^\i^2)0i/^\ks) [PB{k,)PB{k2)6a,PU^l)Pcf{k,) 

+PB{k2)PB{k3)Pad{k2)Pbe{^s)Scf + PB{kl)PB{k^)Pad{^l)5bePcf{^z) 

= ei;'^)(ki)ei/^^)(k3)e(-/^)(k2)[Pi.(A:i)Pi,(A:2) + 2 perms.] 



6TT 



Y^llli'llll' {111} \PB{h)PB{k2) + 2 perms.] 



X 



Ai (kl ) A2 ^2A/' (^2 ) A3 ^2A/" (^s) 



MM'M" 



2 2 2 
M M' M" 



(9.86) 



The delta function is also expanded with the spin spherical harmonics as Eq. (6.14) 

3 



.1=1 



8 / y'^dy 



rO 

^i^^La I Ml M2 M3 



Then, the angular integrals are performed as 

T^* -i^* -i^* 7-AiO-Ai 
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(9.87) 
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(9.88) 



and all the Wigner-3j symbols are summed up as 



M1M2M3 

MM'M" 



Li L2 
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^1 ^2 ^^3 

mi 7712 fns 




(9.89) 



Thus the initial bispectrum (9.85) is rewritten as 
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X 



%S'Kf'^[PB{h)PB{k2) + 2perms.] 



(9.90) 



Comparing the exact initial bispectrum of the tensor modes (9.66) with this equation, we can 
see that the number of the time-integrals and summations in terms of the multipoles decreases. 
This means that corresponding to the pole approximation, the 1-loop calculation (the left panel of 
Fig. 9.4) reaches the tree-level one (the right one of that figure). This approximation seems to be 
applicable to the non-Gaussianity generated from the chi-squared fields without the complicated 
angular dependence 51 . Note that the scaling behaviors of these initial bispectra with respect to 
/ci, ^2 and ^3 are in agreement with that of the local-type non-Gaussianity (4.7). Thus, if the pole 
approximation is valid, we can expect that the PMFs generate the CMB bispectra coming from the 
local-type non-Gaussianity. Via the summation over Ai, A2 and A3 as Eq. (9.73), the approximate 
CMB bispectra of the tensor modes are quickly formulated: 
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where the multipoles are limited as 



(9.91) 



even 



n=l 



and the triangle inequality imposes 



\Li — < < Li + L2 



< 



<il+i2 . 



(9.92) 



(9.93) 



In the next subsection, we compare these approximate spectra with the exact spectra given by 
Eq. (9.74) and evaluate the validity of the pole approximation. 



9.5 Analysis 
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In this subsection, we show the result of the CMB intensity-intensity-intensity spectra induced 
from the auto-correlations of the each-mode anisotropic stress. In order to compute numerically, we 



insert Eq. (9.74) into the Boltzmann code for anisotropics in the microwave background (CAMB) 
151152 . We use the transfer functions shown in Sec. 9.1 In the calculation of the Wigner-3j, 6j and 



9j symbols, we use a common mathematical library called SLATEC 53 and analytical expressions 
in Appendix [Cj 



In Fig. 9.6 we plot the CMB reduced bispectra of these modes defined as 54 



i{ZlZ2Z3) 



(tOOO\ 



'1 75(^1.^2^3) 

^III,hl2l3 ' 



(9.94) 



for ^1 = ^2 = ^3- Here, for comparison, we also write the bispectrum generated from the local-type 
primordial non-Gaussianities of curvature perturbations given by Eq. (4.7). 

From the red solid lines, we can find the enhancement at £ < 100 in tensor-tensor-tensor 
bispectra. It is because the ISW effect gives the dominant signal like in the CMB anisotropics of 
tensor modes 



17, 55 



From the green dashed line, one can see that the peak of the vector- vector- 
vector bispectrum is located at £ ~ 2000 and the position is simil ar to that of the angular power 

At small scales, the vector 



9.1 



spectrum cf^^ induced from the vector mode as calculated in Sec. 
mode contributes to the CMB power spectrum through the Doppler effect. Thus, we can easily 
find that the Doppler effect can also enhance the CMB bispectrum on small scale. From the blue 
dotted lines, we can see that the scalar-scalar-scalar bispectra are boosted around at £ ~ 200 due 
to the acoustic oscillation of the fluid of photons and baryons. On the other hand, as ^ enlarges, 
the spectra are suppressed by the Silk damping effect. These features are also observed in the 
non- magnetic case (the magenta dot-dashed line), however, owing to the difference of the angular 
dependence on the wave number vectors in the source bispectra, the location of the nodes slightly 
differs. Comparing the behaviors between the three spectra arising from PMFs, we confirm that 
the tensor, scalar and vector modes become effective for £ < 100, 100 < £ < 2000 and i > 2000, 
respectively, like the behaviors seen in the power spectra. Thus, for £ < 1000, namely the current 

we expect that 



30 



instrumental limit of the angular resolution such as the PLANK experiment 
the auto- and cross-correlations between the scalar and tensor modes will be primary signals of 
PMFs in the CMB bispectrum. 



The overall amplitudes of 2^ and bf^J^^^ seem to be comparable to 



C 



III 



3/2 



and 



C 



(T) 



-1 3/2 



However, we find that the amplitude of bf^]^^J, is smaller than the above expectation. This is 
because the configuration of multipoles, corresponding to the angles of wave number vectors, is 
limited to the conditions placed by the Wigner symbols. We can understand this by considering 
the scaling relation with respect to i at high i. If the magnetic power spectrum given by Eq. 

L" 



(9.4) is close to the scale-invariant shape, the configuration that satisfies L 



and L' ~ 1 



contributes dominantly in the summations. Furthermore, the other multipoles are evaluated as 



) -^2 



L'r, ~ L'i ~ 1 



(9.95) 



from the triangle conditions described in Appendix 



C 



Then we can find b^^^^e 
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Figure 9.6: Absolute values of the normalized reduced bispectra of temperature fluctuations for a 
configuration ii = £2 = h = The red solid, green dashed, and blue dotted lines correspond to 
the spectra generated from the auto-correlations of the PMF tensor, vector, and scalar anisotropic 
stresses for ub = —2.9, respectively. The upper (lower) lines of the red solid and blue dotted 
lines are calculated when Tu/tb = 10^^(10^). The magenta dot-dashed line expresses the spectrum 
sourced from the primordial non-Gaussianity with f^^^^ = 5. The strength of PMFs is fixed to 
-BiMpc = 4.7nG and the other cosmological parameters are fixed to the mean values limited from 



WMAP-7yr data reported in Ref. 29 



< 1000, where we have also used the following relations 



edkr}'^\k)jeMro-n))oci 



V ("B+1) 




(9.96) 



which, except for the first relation, are also coming from the triangle conditions of the Wigner 3-j 



symbols^ 
in Sec. 



9.1 



herefore, combining with the scaling relation of the CMB power spectrum mentioned 



we find that h 



{VVV) . 



In Fig. 



9.7 



/// m is suppressed by a factor 

'{VVV) 



(ns-l)/2 fro^ C^V)^/\ 



we show for 



?i = ^2 = ^3 for the different spectral index hb- Red solid and 
green dashed lines correspond to the bispectrum with the spectral index of the power spectrum 
of PMFs fixed as ub = —2.9 and —2.8, respectively. From this figure, we find that the CMB 
bispectrum becomes steeper if becomes larger, which is similar to the case of the power spec- 
trum. These spectra trace the scaling relation in the above discussion. These will lead to another 
constraint o n th e strength of PMFs. 



In Figs, 
with setting 



9.8 



and 



9.9 



we show the reduced bispectrum bf^^^J^ 



respectively. From Fig. 



9.8 



we can see that b 



and bfjJJ^g^^^ with respect to 

(VVV) 



for 



> 100 is 
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Figure 9.7: Absolute values of the normalized reduced temperature-temperature-temperature spec- 
tra arising from the auto-correlation between the PMF vector anisotropic stresses for a configu- 
ration ii = £2 = ^3 = ^- The lines correspond to the spectra generated from vector anisotropic 
stress for = —2.9 (red solid line) and —2.8 (green dashed line), and primordial non-Gaussianity 
with /nl*^' = 5 (blue dotted line). The strength of PMFs is fixed to -BiMpc = 4.7nG and the other 



cosmological parameters are identical to the values used in Fig. 9.6 



nearly flat and given as 

hih + mih + lWillUl - 2x10-19(^1 . (9.97) 
We can understand this by the analytical evaluation as follows. As mentioned above, in the sum- 



-BlMpc A 

4.7nG / 



mations of Eq. (9.74), the configuration that L ~ £i,L' ~ 1 and L" ~ £3 contributes dominantly. 



By using this and the approximations that 



'1 ) 



(9.98) 



which again come from the triangle conditions from the Wigner symbols, the scaling relation 
of £3 at large scale is evaluated as b'^n'ilxhh °^ ^1^^^ ■ From this estimation we can find that 
£i(£i + 1)4(4 + ^)bfnJXh °^ = -2.9, and for ub = -2.8, respectively, which 



match the behaviors of the bispectra in Fig. |9.8 
b 



From Fig. 9.9, we can also see that if the PMF spectrum obeys the nearly scale-invariant shape, 
nJJLia fo'^ h,^2,h< 100 is given by 



4(4 + 1)4(4 + l)ldU 



(130 - 6) X 10-- ( ^ , 



(9.99) 



where the factor 130 corresponds to the t^/tb = 10^'' case and 6 corresponds to 10^. In order to 
obtain a rough constraint on the magnitude of the PMF, we compare the bispectrum induced from 
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Figure 9.8: Absolute values of the normalized reduced temperature-temperature-temperature bis- 
pectra induced by the auto-correlation between the PMF vector anisotropic stresses and generated 
by primordial non-Gaussianity given by Eq. (9.100) as a function of £3 with li and £2 fixed to some 



value as indicated. Each parameter is fixed to the same value defined in Fig. 9.6 



the PMF with that from the local-type primordial non-Gaussianity in the curvature perturbations, 



which is typically estimated as 56 



4 X lO""*^*^ -flocal 



NL 



(9.100) 



By comparing this with Eq. (|9.99|), the relation between the magnitudes of the PMF with the 

(9.101) 



nearly scale-invariant power spectrum and /nl*^' is derived as 

-BlMpc 



InG 



(1.22 -2.04) l/^lj'^'^^l^/'' 



Using the above equation, we can obtain the rough bound on the PMF strength. As shown in 



Fig. 9.6, because the tensor bispectrum is highly damped for £ > 100, we should use an upper 



bound on obtained by the current observational data for ^ < 100, namely /nl^^ < 100 57 



This value is consistent with a simple prediction from the cosmic variance 54 . From this value, 
we derive -BiMpc < 2.6 — 4.4nG. 

From here, let us discuss the validity and possibility of the CMB bispectra under the pole 



approximation (9.91). Figure 9.10 shows the shapes of the CMB tensor-tensor-tensor spectra 



based on the exact form (9.74) and approximate one (9.91). Both spectra seem to have a good 



agreement in the shape of the i space. To discuss more precisely, using the correlation 



b-b' bmb'i 



m 5 



we calculate a correlation coefficient as 

^ex . ^app 



^(^^ex . ^ox^j^^app . ^app) 



0.99373 , 



(9.102) 



(9.103) 
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Figure 9.9: Absolute values of the normalized reduced temperature-temperature-temperature bis- 
pectra induced by the auto-correlation between the PMF tensor anisotropic stresses and generated 



from primordial non-Gaussianity in curvature perturbations given by Eq. (9.100) as a function of 
£3 with ^l = £2- Each parameter is identical to the values defined in Fig. 9.6 



where b'^^ and b^^^ are the exact and approximate reduced bispectra, respectively. This fact, which 
this quantity approaches unity, implies that the pole approximation can produce an almost exact 
copy. An unknown parameter, a, is derived from the relation as 

a = ""-^ ^ ■ ^"^(" = ^) = 2991 (9 104) 

6-PP(a = l) b-PP{a = l)-b-PP{a = l) ' 



The cases other than the tensor-tensor-tensor spectrum will be presented in Ref. 35 . 

As shown in the previous subsections, the CMB bispectra from PMFs arise from the six-point 
correlation of the Gaussian magnetic fields and have one-loop structure due to the summation over 
the additional multipoles. Hence, it takes so long hours to estimate all £'s contribution and it is 
actually impossible to compute the signal-to-noise ratio. However, using the pole approximation, 
since the summation reaches the tree-level calculation, we will obtain more precise bound through 
the estimation of the signal-to-noise ratio including the contribution of the cross-correlations be- 



tween scalar and tensor modes 35 . 



9.6 Summary and discussion 



In this section, on the basis of our recent works 24-26 , we presented the all-sky formulae for 
the CMB bispectra induced by the scalar, vector, and tensor non-Gaussianities coming from the 
PMFs by dealing with the full-angular dependence of the bispectrum of the PMF anisotropic 
stresses. Then, expressing the angular dependence with the spin-weighted spherical harmonics and 
converting the angular integrals into the Wigner symbols were key points of the formulation. From 
the practical calculation, it is found that the CMB bispectra from the magnetic tensor, scalar, and 
vector modes dominate at large {i < 100), intermediate (100 < < 2000), and small {i > 2000) 
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Figure 9.10: Absolute values of the normalized reduced bispectra of temperature fluctuation for a 
configuration £i = £2 = ^3 = ^- The red solid and green dashed lines represent the exact and ap- 
proximate spectra arising from the tensor-tensor-tensor correlation of the PMF anisotropic stresses 
for -BiMpc = 4.7nG,nB = —2.9 and t^/tb = 10^^, respectively. The cosmological parameters are 



identical to the values defined in Fig. 9.6 



scales. For the discussion about the shape of the non-Gaussianity in the PMF anisotropic stresses, 
we performed the pole approximation, which is the evaluation of the convolutions at around the 
divergence points of the integrands, and found that the bispectra of the PMF anisotropic stresses 
are classified as the local-type configuration. Owing to this, we had some significant signals of the 
CMB bispectra on the squeezed limit also in the multipole space. Compared with the exact formula, 
the approximate one reduces the computing time, hence we expect the use for the calculation of 



the signal-to-noise ratio 35 . We also investigated the dependence of the CMB bispectrum on the 
spectral index of the PMF power spectrum and confirmed that the CMB bispectrum induced from 
the PMFs is sensitive to it. Since the characteristic scale varies with the value of the spectral 
index, it is important to consider not only the contribution from the scalar mode, but also those 
from the vector and tensor modes. 

By translating the current bound on the local-type non-Gaussianity from the CMB bispectrum 
into the bound on the amplitude of the magnetic fields, we obtain a new limit: -BiMpc < 2.6— 4.4nG. 
This is a rough estimate coming from the large scale information of the tensor mode and a precise 
constraint is expected if one considers the full i contribution by using an appropriate estimator of 
the CMB bispectrum induced from the primordial magnetic fields. 

Because of the complicated discussions and mathematical manipulations, here we restrict our 
numerical results to the intensity bispectra of auto-correlations between scalar, vector and tensor 
modes despite the fact that our formula for the CMB bispectra (9.69 ) contains the polarizations and 
the cross-correlations between scalar, vector and tensor modes. However, like the non-magnetic 
case 



58 



the modes other than our numerical results, such as -Sj^^^^^^^^^ , will bring in more 



reasonable bounds on the PMFs 35 . Furthermore, the effect on the CMB four-point correlation 



(trispectrum) is just beginning to be roughly discussed 59 . Applying our studies, this should be 
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taken into account more precisely. 

References 

[1] M. L. Bernet, F. Miniati, S. J. Lilly, P. P. Kronberg, and M. Dessauges-Zavadsky, Strong 
magnetic fields in normal galaxies at high redshifts, Nature 454 (2008) 302-304, 
[arXiv: 0807 .334 7] . 

[2] A. M. Wolfe, R. A. Jorgenson, T. Robishaw, C. Heiles, and J. X. Prochaska, An 84 
microGauss Magnetic Field in a Galaxy at Redshift z=0.692, Nature 455 (2008) 638, 
[arXiv: 0811 .240 81. 

[3] P. P. Kronberg et ai, A Global Probe of Cosmic Magnetic Fields to High Redshifts, 
Astrophys. J. 676 (2008) 70-79, [ arXiv : 0712 . 0435] . 

[4] L. M. Widrow, Origin of Galactic and Extragalactic Magnetic Fields, Rev. Mod. Phys. 74 
(2002) 775-823, |a stro-ph/0207240] . 

[5] J. Martin and J. Yokoyama, Generation of Large- Scale Magnetic Fields in Single-Field 
Inflation, JC^P 0801 (2008) 025, [arXi v: 071 1.4307] . 

[6] K. Bamba and M. Sasaki, Large-scale magnetic fields in the inflationary universe, JGAP 
0702 (2007) 030, |astro- ph/061170l|. 

[7] T. Stevens and M. B. Johnson, Theory of Magnetic Seed-Field Theory of Magnetic Seed- 
Field Generation during the Gosmological First-Order Electroweak Phase Transition, 
[arXiv: 1001 .36941 

[8] T. Kahniashvili, A. G. Tevzadze, and B. Ratra, Phase Transition Generated Gosmological 
Magnetic Field at Large Scales, Astrophys. J. 726 (2011) 78, [arXiv: 090T7 bl97|. 

[9] K. Ichiki, K. Takahashi, H. Ohno, H. Hanayama, and N. Sugiyama, Gosmological Magnetic 
Field: a fossil of density perturbations in the early universe. Science. 311 (2006) 827-829, 
|astro-ph/0603631 . 

[10] S. Maeda, S. Kitagawa, T. Kobayashi, and T. Shiromizu, Primordial magnetic fields from 
second-order cosmological perturbations:Tight coupling approximation. Glass. Quant. Grav. 
26 (2009) 135014, |arXiv: 0805.0169] . 

[11] E. Fenu, C. Pitrou, and R. Maartens, The seed magnetic field generated during 

recombination, Mon.Not.Roy.Astron.Soc. 414 (2011) 2354-2366, [arXiv: 1012.2958 . 

[12] K. Subramanian and J. D. Barrow, Microwave Background Signals from Tangled Magnetic 
Fields, Phys. Rev. Lett. 81 (1998) 3575-3578, |astro-ph/980326lj . 

[13] R. Durrer, T. Kahniashvili, and A. Yates, Microwave Background Anisotropics from Alfven 
waves, Phys. Rev. D58 (1998) 123004, |astro-ph/9807089]. 

[14] A. Mack, T. Kahniashvili, and A. Kosowsky, Vector and Tensor Microwave Background 
Signatures of a Primordial Stochastic Magnetic Field, Phys. Rev. D65 (2002) 123004, 
|astro-ph/0105504 . 



152 



REFERENCES 



A. Lewis, CMB anisotropics from primordial inhomogeneous magnetic fields, Phys. Rev. 
D70 (2004) 043011, |astro-ph/0406096|. 

D. Paoletti, F. Finelli, and F. Paci, The full contribution of a stochastic background of 
magnetic fields to CMB anisotropics, Mon. Not. Roy. Astron. Soc. 396 (2009) 523-534, 
|arXiv:0811.0230| . 



J. R. Shaw and A. Lewis, Massive Neutrinos and Magnetic Fields in the Early Universe, 



Phys. Rev. D81 (2010) 043517, [arXiv: 0911 .2714 



I. Brown and R. Crittenden, Non-Gaussianity from Cosmic Magnetic Fields, Phys. Rev. 
D72 (2005) 063002, |astro-ph/0506570| . 



I. A. Brown, Primordial Magnetic Fields in Cosmo/og'y, arXiv: 0812. 1781, Ph.D. Thesis. 



T. R. Seshadri and K. Subramanian, CMB bispectrum from primordial magnetic fields on 
large angular scales, Phys. Rev. Lett. 103 (2009) 081303, [ arXiv: 0902 .4066] . 

C. Caprini, F. Finelh, D. Paoletti, and A. Riotto, The cosmic microwave background 
temperature bispectrum from scalar perturbations induced by primordial magnetic fields, 
JC^P 0906 (2009) 021, [ arXiv: 0903. 1420] . 

R.-G. Cai, B. Hu, and H.-B. Zhang, Acoustic signatures in the Cosmic Microwave 
Background bispectrum from primordial magnetic fields, JCAP 1008 (2010) 025, 
|ar Xiv: 1006 .2985] . 

P. Trivedi, K. Subramanian, and T. R. Seshadri, Primordial Magnetic Field Limits from 
Cosmic Microwave Background Bispectrum of Magnetic Passive Scalar Modes, Phys. Rev. 
D82 (2010) 123006, [a rXiv: 1009 .27 241. 

M. Shiraishi, D. Nitta, S. Yokoyama, K. Ichiki, and K. Takahashi, Cosmic microwave 
background bispectrum of vector modes induced from primordial magnetic fields, Phys. Rev. 



D82 (2010) 121302, [arXiv: 1009.3632 



M. Shiraishi, D. Nitta, S. Yokoyama, K. Ichiki, and K. Takahashi, Computation approach for 
CMB bispectrum from primordial magnetic fields, Phys. Rev. D83 (2011) 123523, 
|arXiv :110"r5287|. 

M. Shiraishi, D. Nitta, S. Yokoyama, K. Ichiki, and K. Takahashi, Cosmic microwave 
background bispectrum of tensor passive modes induced from primordial magnetic fields, 
Phys. Rev. D83 (2011) 123003, | arXiv: 1103.4103 |. 

M. Shiraishi, D. Nitta, S. Yokoyama, K. Ichiki, and K. Takahashi, CMB Bispectrum from 
Primordial Scalar, Vector and Tensor non-Caussianities, Prog. Theor. Phys. 125 (2011) 
795-813, larXiv: 1012. 1079] . 

T. Kahniashvili and G. Lavrelashvili, CMB two- and three-point correlation functions from 
Alfven waves, arXiv: 1010.45431 

WMAP Collaboration, E. Komatsu et al.. Seven- Year Wilkinson Microwave Anisotropy 
Probe (WMAP) Observations: Cosmological Interpretation, Astrophys. J. Suppl. 192 (2011) 
18, [arXiv: 1001 .45381 . 



REFERENCES 



153 



PLANCK Collaboration, The Scientific programme of planck, astro-ph/0604069 



M. Shiraishi, Parity violation of primordial magnetic fields in the CMB bispectrum, J CAP 
1206 (2012) 015, I arXiv: 1202.2847]. 

J. M. Bardeen, Gauge Invariant Cosmological Perturbations, Phys.Rev. D22 (1980) 
1882-1905. 

K. Kojima, T. Kajino, and G. J. Mathews, Generation of Gurvature Perturbations with 
Extra Anisotropic Stress, JGAP 1002 (2010) 018, | arXiv:0910. 1976] . 

A. Lewis, Observable primordial vector modes, Phys. Rev. D 70 (Aug, 2004) 043518. 

M. Shiraishi, D. Nitta, S. Yokoyama, and K. Ichiki, Optimal limits on primordial magnetic 
fields from GMB temperature bispectrum of passive modes, JGAP 1203 (2012) 041, 
|arXiv: 1201 . 0376]. 

D. G. Yamazaki, K. Ichiki, T. Kajino, and G. J. Mathews, Effects of a Primordial Magnetic 
Field on Low and High Multipoles of the GMB, Phys.Rev. D77 (2008) 043005, 
I arXiv : 080 1.2 572]. 

F. Finelli, F. Paci, and D. Paoletti, The Impact of Stochastic Primordial Magnetic Fields on 
the Scalar Contribution to Gosmic Microwave Background Anisotropics, Phys.Rev. D78 
(2008) 023510, [arXiv: 0803. 1246]. 



I. A. Brown, Goncerning the statistics of cosmic magnetism, arXiv: 1005.2982 



M. Giovannini and K. E. Kunze, Generalized GMB initial conditions with pre-equality 
magnetic fields, Phys.Rev. D77 (2008) 123001, | arXivTo'8027l 053 1 . 

K. E. Kunze, GMB anisotropics in the presence of a stochastic magnetic field, Phys.Rev. 
D83 (2011) 023006, |arXiv: 1007.3163]. 

D. Paoletti and F. Finelli, GMB Constraints on a Stochastic Background of Primordial 
Magnetic Fields, Phys.Rev. D83 (2011) 123533, [ arXiv: 1005.014 8|. 

D. G. Yamazaki, K. Ichiki, T. Kajino, and G. J. Mathews, New Constraints on the 
Primordial Magnetic Field, Phys. Rev. D81 (2010) 023008, [arXiv: 1001 . 2012]. 



J. R. Shaw and A. Lewis, Constraining Primordial Magnetism, arXiv: 1006.4242 



D. Yamazaki, K. Ichiki, T. Kajino, and G. J. Mathews, Constraints on the Evolution of the 
Primordial Magnetic Field from the Small-Scale Cosmic Microwave Background Angular 
Anisotropy, Astrophys. J. 646 (2006) 719-729, tastro-ph/0602224| . 

K. Jedamzik, V. Katalinic, and A. V. Olinto, Damping of Cosmic Magnetic Fields, Phys. 
Rev. D57 (1998) 3264-3284, |astro- ph/9606080]. 



K. Subramanian and J. D. Barrow, Magnetohydrodynamics in the early universe and the 
damping of noninear Alfven waves, Phys. Rev. D58 (1998) 083502, |astro-ph/9712083| . 



W. Hu, Angular trispectrum of the cosmic microwave background, Phys. Rev. D64 (2001) 



083005, |astro-ph/0105117 



154 



REFERENCES 



R. Gurau, The Ponzano-Regge asymptotic of the 6j symbol: an elementary proof, Annates 
Henri Poincare 9 (2008) 1413-1424, |arXiv: 080873533] . 

H. A. Jahn and J. Hope, Symmetry properties of the wigner 9j symbol, Phys. Rev. 93 (Jan, 
1954) 318-321. 

"The wolfram function site." http://functions.wolfram.com/. 

D. H. Lyth, Non-gaussianity and cosmic uncertainty in curvaton-type models, JCAP 0606 
(2006) 015, [astro-ph/0602285] . 

A. Lewis, A. Challinor, and A. Lasenby, Efficient Computation of CMB anisotropics in 



closed FRW models, Astrophys. J. 538 (2000) 473-476, |astro-ph/9911177i . 



"Slatec common mathematical library." http://www.netlib.org/slatec/. 

E. Komatsu and D. N. Spergel, Acoustic signatures in the primary microwave background 
bispectrum, Phys. Rev. D63 (2001) 063002, |astro-ph/0005036l. 

J. R. Pritchard and M. Kamionkowski, Cosmic microwave background fluctuations from 
gravitational waves: An analytic approach. Annals Phys. 318 (2005) 2-36, 
|astro-ph/041258l|. 

A. Riotto, The Quest for Non-gaussianity, in Inflationary Cosmology (M. Lemoine, 

J. Martin, & P. Peter, ed.), vol. 738 of Lecture Notes in Physics, Berlin Springer Verlag, 

pp. 305-+, 2008. 

K. M. Smith, L. Senatore, and M. Zaldarriaga, Optimal limits on f^^ff^ from WMAP 5-year 
data, JC^P 0909 (2009) 006, |arXiv: 0901 . 2572]. 

D. Babich and M. Zaldarriaga, Primordial bispectrum information from CMB polarization, 
Phys. Rev. D70 (2004) 083005, | astro-ph/0408455 1 . 

P. Trivedi, T. Seshadri, and K. Subramanian, Cosmic Microwave Background Trispectrum 
and Primordial Magnetic Field Limits, Phys. Rev. Lett. 108 (2012) 231301, 
I arXiv: 1111.0744]. 



155 



10 Conclusion 

The main purpose of this thesis was to present the formahsm for the CMB bispectrum induced by 
the non-Gaussianities not only in the standard scalar-mode perturbations but also in the vector- 
and tensor-mode ones where the violation of the rotational or parity invariance is also involved, 
and some attempts to prove the nature of the early Universe by applying our formalism. To do 
this, we have discussed the following things. 

In Sec. [T| we gave the introduction of this thesis. Then, we quickly summarized the history 
of the Universe, the paradigm in the early Universe, and the concept of this thesis. In Sec. [2| we 
summarized how to generate the curvature perturbations and gravitational waves and the consis- 
tency relations in the slow-roll inflation. In Sec. [3} we showed how to construct the a^m's generated 
from the primordial scalar, vector and tensor sources in order to formulate the CMB bispectrum 
easily. We also summarized the constraints on several key parameters, which characterize the 
nature of inflation and the dynamics of the Universe, obtained from the current CMB data. In 
Sec. |4| we focused on the topic of the primordial non-Gaussianities. In Sec. [5| we gave the general 
formulae for the CMB bispectrum coming from not only scalar-mode but also vector- and tensor- 
mode perturbations, which includes both the auto- and cross- correlations between the intensity 
and polarizations. Next, applying this formalism, we computed the CMB bispectra from several 
kinds of the non-Gaussianities. In Sec. |6| we treated the two scalars and a graviton correlator and 
obtained the CMB bispectrum including the tensor-mode perturbation. Here, we had a bound 
on the nonlinear scalar-scala-tensor coupling by the computation of the signal-to-noise ratio. In 
Sec. [7| we considered the non-Gaussianity which has the preferred direction. Through the analysis, 
we found that the finite signals arise from the multipoles except for the triangle inequality. We 
furthermore confirmed that these special signals are comparable in magnitude with the signals 
keeping the triangle inequality. In Sec. |8j we dealt with the graviton non-Gaussianity arising from 
the parity-conserving and parity-violating Weyl cubic terms. Calculating the CMB intensity and 
polarization bispectra, we clarified that the intensity-intensity-intensity spectrum from the parity- 
violating non-Gaussianity obeys the condition as Yln=i = odd. These configurations will be very 
beneficial to check the parity violation of the Universe in the non-Gaussian level observationally. 
In Sec. [9} we took into account the effect of the non-Gaussianities due to the primordial magnetic 
fields. Depending quadratically on the magnetic fields, the magnetic anisotropic stresses obey the 
chi-square distributions. Since these non-Gaussian anisotropic stresses become the sources of the 
CMB fluctuations, their bispectra have the finite values. Computing the CMB intensity-intensity- 
intensity spectra, we clarified that the tensor (vector) mode dominates at large (small) scales and 
the scalar mode shows up at intermediate scales. By the pole approximation, we also found that 
the bispectrum of the magnetic anisotropic stresses is similar to the local-type bispectrum. Com- 
paring the theoretical results with the observational limit on the local-type non-Gaussianity, we 
obtained a bound on the strength of the magnetic fields, -BiMpc < 2.6 — 4.4nG. We expect that 
this bound will be updated by considering the impacts of the cross-correlations between scalar, 
vector and tensor modes, and the additional information from polarizations. 

Our formalism for the CMB bispectrum is general enough to be applicable to the non-Gaussian 
sources other than the above ones. Moreover, this will be easily extended to the higher-order 
correlations. Therefore, the studies in this thesis will be very beneficial to quest for the true 
picture of the origin of the Universe. 
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A SPIN- WEIGHTED SPHERICAL HARMONIC FUNCTION 



A Spin-weighted spherical harmonic function 

Here, we review the theory of spin-weighted functions and their expansion in spin-s spherical 
harmonics. In the past, these functions was mainly applied to the analysis of the gravitational 
wave (see e.g. Ref. [^). This discussion is based on Refs. [2]-[4]. 

For any given direction on the sphere specified by the angles {0,(j)), one can define three or- 
thogonal vectors, one radial and two tangential to the sphere. Let us denote the radial direction 
vector with n and the tangential with ei, 62- The latter two are only defined up to a rotation 
around n. 

A function sf{G, 0) defined on the sphere is said to have spin-s if under a right-handed rotation 
of (61,62) by an angle ip it transforms as s/'(^, 0) = e~**'^s/(6', 0). For example, given an arbitrary 
vector a on the sphere the quantities, a ■ ei + ia ■ 62, n ■ a and a ■ ei — ia ■ 62 have spin 1,0 and -1, 
respectively. 

A scalar field on the sphere can be expanded in spherical harmonics, Yim{0,(j)), which form a 
complete and orthonormal basis. These functions are not appropriate to expand spin weighted 
functions with s 7^ 0. There exist analog sets of functions that can be used to expand spin-s 
functions, the so called spin-s spherical harmonics sYim{0,(f)). These sets of functions (one set for 
each particular spin) satisfy the same completeness and orthogonality relations, 

27r ^1 

dcOs9sYif^,{9,(t>)sYim{9,(t)) = 6v^l5m'.m , 

(A.l) 

sYiUe, 0).n^(^', 0') = 5(0 - 0')5(cos e - COS 6'). 

Im 

An important property of spin-s functions is that there exists a spin raising (lowering) operator 
d {d ) with the property of raising (lowering) the spin-weight of a function, [d sf)' = e"**^*"'"^^^^ sf, 
{d sf)' = e~'^^'^~^^'^ 6 sf- Their explicit expression is given by 



d sf{e, 0) = - sin^ e[de + t CSC ed^] sm-^ Osfie, 0) 
5 sf{e, 0) = - sin-^ e[dg~i CSC 68^] sin^ J(e, 0) 



(A.2) 



^ '2f{f^, 0) = [-d, + j [(1 - ^^)2/(/i, 0)] 

^ '-2/(^, 0) = (-5, - ' [(1 - M^)-2/(/i, 



Specifically, the d and d operators acting twice on the spin-±2 function ±2/(/^, 0) such as the 
CMB polarization fields can be expressed as 



(A.3) 



where /i = cos 9 and ±2f{fJ', 0) = ±2/(/^)e*'"'^. With the aid of these operators, we can express gYim 
in terms of the spin-0 spherical harmonics Yim, which are the usual spherical harmonics, 

(A.4) 
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Table A.l: Dipole (/ = 1) 


harmonics for spin-0 and 1 




m 




2^2m 


±2 




(1-f cos0)%™<^ 


±1 


— m 


sin^ cosOe"'''^ 


— 1 sin 6* (1 — m cos ^) 







1 

2 y 


/^(3cos2^-l) 


sin^^ 

4 y DVT 



Table A. 2: Quadrupole {I = 2) harmonics for spin-0 and 2. 



where these equations contain 



6 sYim = [(/ - + s + 1)]2 ,+iyi„ , 

^s>/m = -[(/ + s)(/-S + l)]^s-l>;. 

^.5,yz^ = -(/-s)(/ + s + i),ii^ . 



(A.5) 



The above relations lead to an explicit expression for the spin weighted spherical harmonic 
function: 



sYlmiQ-, ' 



(/ + m)!(/-m)!2/ + l 
_ {l + sy.{l-s)\ 



1/2 



sin2'(0/2) 



X 



l-s 
r 



I + s 
r + s — m 



(A.6) 



The following properties of spin-weighted harmonics are also useful 



sYirr^TT - 6*, + TT 



(A.7) 



Finally, we give the specific expressions for some simple cases in Tables A.l and A.2 
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B WIGNER D-MATRIX 



B Wigner D-matrix 



we introduce the properties of the Wigner D-matrix -D^"*^/, 



Here, on the basis of Refs. js 5 6 
which is the unitary irreducible matrix of rank 21 + 1 that forms a representation of the rotational 
group as SU{2) and 5*0(3). With this matrix, the change of the spin wighted spherical harmonic 
function under the rotational transformation as n — )■ Rh is expressed as 

sY;^{m) = Y,D^2AR)sYUh) . (B.i) 

m' 

This satisfies the relation as 

D^2M = i-^r-"^' D^Xm'iR) = dI^'JR-') ■ (B.2) 

When we express the rotational matrix with three Euler angles (a, /3, 7) under the z — y — z 
convention as 

(cos a cos (3 cos 7 — sin a sin 7 — cos (3 sin 7 cos a — cos 7 sin a cos a sin /3 \ 
cos a sin 7 + cos 7 cos /3 sin a cos a cos 7 — cos /3 sin a sin 7 sin /3 sin a , (B.3) 
— cos 7 sin 13 sin 7 sin j3 cos /3 J 

we can write a general relationship between the Wigner D-matrix and the spin weighted spherical 
harmonics as 



^£(«,A7) = i-ir^^-sY;^W,a)e-^ . (B.4) 

Like the spin weighted spherical harmonics, there also exists the orthogonality of the Wigner 
Z)-matrix as 

da j dcosP d'yDll}*,{a,P,'y)D^^{a,P,'y) = ^, , ,i5m\m5s' ,s ■ (B.5) 
^-1 ^0 zi + 1 
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C Wigner symbols 

Here, we briefly review the useful properties of the Wigner-3j, 6j and 9j symbols. The following 



discussions are based on Refs. 5l 7 12 . 



C.l Wigner-3j symbol 

In quantum mechanics, considering the coupling of two angular momenta as 

l3 = ll+l2, (C.l) 

the scalar product of eigenstates between the right-handed term and the left-handed one, namely, 
a Clebsch-Gordan coefficient, is related to the Wigner-3j symbol: 

k k ^3 \ (-l)^i-'2+m3 {l^m^l^m2\{kk)hm:,) 
nil rn2 -m^ J ~ ^/2k~+l 

This symbol vanishes unless the selection rules are satisfied as follows: 

\mi\ < k , 1^21 < k , l^^sl < ^3 , mi + m2 = ma , 
1^1 — ^2! ^ k ^ k + k (the triangle condition) , k + k + k ^ ■ 

Symmetries of the Wigner-3j symbol are given by 

k k \ = ( ^2 k k \^fi\j:Uhf k k k 

mi m2 ms J y m2 mi ms / v 

(odd permutation of columns) 
k k k \ _ I k k k 

1712 ^3 ITT'l J \ TTT'S '^'1 ^2 

(even permutation of columns) 

k k k 
-nil — m2 —ma 

(sign inversion of nii,m2,ni3) . (C.4) 



(C.3) 



The Wigner-3j symbols satisfy the orthogonality as 



J2{2k + l)( 

/a 7713 



^1 k k \ f k k l. 



'3 



mi m2 ms I \ ni[ m'2 mg J "^2-^2 ' 



(C.5) 



(2/3 + I) y ( k \( k k k )=S,^,6^,rr^' . 

V 3 ^ ^ \ mi m2 ms / V mi m2 m'o / '"''^ "'^''"s 

mir?i2 ^ ^ ' 

For a special case that Yl^=i ~ even and mi = m2 = = 0, there is an analytical expression as 

(k k k\_^ ^p_^=i. {EU'i)Wi-h + k + k)\V{k^k+W (^_g) 

\ J / -h+f2+/3 N| ( h-h+l3 )l ( h+h-h )] / ^ , 
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C WIGNER SYMBOLS 



This vanishes for Yli=i ~ odd. The Wigner-3j symbol is related to the spin-weighted spherical 
harmonics as 

2 

JJ^,y«,m,(n) = ^ ,^Yil^^{n)Ii^\ j , (C.7) 



i=l h'msss 

which leads to the "extended" Gaunt integral including spin dependence: 



h h h 

h h \ jji^ jYi^ ^2 



(C.8) 



Here P^'^'^ = , / (2'i+i)(2^2+i)(2Z3+i) ( h h h 
hhh V 471 I 5^ S2 S3 



C.2 Wigner-6j symbol 

Considering two other ways in the coupling of three angular momenta as 

I5 = I1 + I2 + I4 (C.9) 

= I3 + I4 (C.IO) 

= li + le, (C.ll) 

the Wigner-6j symbol is defined using a Clebsch-Gordan coefficient between each eigenstate of I5 
corresponding to Eqs. (C.IO) and ( |C.ll ) as 



^4 h k J ^^(2/3 + 1) (2/6 + 1) 

This is expressed with the summation of three Wigner-3j symbols: 



(C.12) 



(_i)Et4 ( ^5 ^1 ^6 \ 

^ I 7715 —mi — mg / 



h h h \ I h h U \ f U h h 

777-6 —"7714 / I —7773 — 7775 



h h ^3 \ J h h h 
7771 m2 ms J \ h h k 

hence, the triangle conditions are given by 

1^1 ~ ^2! ^ ^3 ^ ^1 + ^2? 1^4 ~ ^5! ^ ^3 ^ ^4 + ^5 ) 
1^1 ~ ^sl ^ ^6 ^ ^1 + ^5) 1^4 ^ /2I ^ ^6 ^ ^4 + ^2 • 

The Wigner-6j symbol obeys 24 symmetries such as 



(C.13) 



(C.14) 



(permutation of columns) 



/l ^2 ^3 1 n2 /l /3 1 f k k h 

h h h I 1 ^5 h h j \ h h U 

h h h \ _ \ h h h 
h h h } \ U h h 
(exchange of two corresponding elements between rows). (C.15) 

Geometrically, the Wigner-6j symbol is expressed using the tetrahedron composed of four triangles 



obeying Eq. (C.14). It is known that the Wigner-6j symbol is suppressed by the square root of 



the volume of the tetrahedron at high multipoles. 



C.3 Wigner-9j symbol 
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C.3 Wigner-9j symbol 

Considering two other ways in the couphng of four angular momenta as 



where I3 = li + I2, le 



Ic 



14 + 15,17 = ll 



ll + 12 + 
13 + 16 

I7 + Is , 



L + Is 



(C.16) 
(C.17) 
(C.18) 

the Wigner 9j symbol expresses a 



I4, Ig = I2 

Clebsch-Gordan coefficient between each eigenstate of I9 corresponding to Eqs. (C.17) and (C.18) 
as 



^1 h ^3 
h h h 
h h h 



v/(2/3 + 1) (2/6 + 1) (2/7 + 1) (2/8 + 1) 



This is expressed with the summation of five Wigner-3j symbols: 



E 

m4m5m6 
m^mgrng 



h 


h 


)( 


h 




me 




m-j 


h 


h 


)( 


h 


mi 


mi 






-{ 


h 


h 


^3 




mi 


m2 





la lo 



la In 



l& h h 
me mg 777.3 




and that of three Wigner-6j symbols: 




5^(-l)2-(2x + l)| 



l\ Z4 Ij 



k k 
X k 



k k k 
X k k 



hence, the triangle conditions are given by 

|/i - /2I <k<k + k , \k - k\ <k<k + k 
\k - k\ <k<k + k , \k - k\ <k<k + k 
The Wigner- 9j symbol obeys 72 symmetries: 



\k-k\ <k<k + k 
\k-k\ <k<k + k 




-1) 





k k k] ( k k 

^-^H k k k [ = (-i)^-''<^ k k 

k k k ) k k 

(odd permutation of rows or columns) 

k ] \ k k k 
k ( ~ { k k k 
k } (k k k 

(even permutation of rows or columns) 

k ] \ k k k 
k ( ~ ] k k k 
k J k k k 

(reflection of the symbols) . 



(C.19) 



(C.20) 



; (C.21) 



(C.22) 



(C.23) 
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C WIGNER SYMBOLS 



C.4 Analytic expressions of the Wigner symbols 

Here, we show some analytical formulas of the Wigner symbols. 
The / symbols, which are defined as 



^.,...3 = / (2/1 + 1) (2/2 + 1) (2/3 + 1) / /i h h \ 24^ 
'i'2'3 V 47r \si S2 ) ' ^ ' ' 



are expressed as 



X 



7-0 _ ■ / n?=i(2/i + 1) . -,NEf-i4^ 

(E-=l I)! V^Hr+ ^2 + l3)Wih -l2 + h)Wih + /2 - /3)! 

(^)! (i^)!^(EL /^ + l)! 

(for /i + /2 + /s = even) 
(for /i + /a + /a = odd) , (C.25) 

rOl-l _ l^r Uh + l /(^2-l)(/2 + l) 



15 , / ^2 + 1/2 



-1) A 7^ TTT^TT^^wT^^ (for /i = /a, /s = 2) 



167r' ' V (/2-l/2)(/2 + 3/2) 




./^(-i)'3+iy;^ (fo, /i = /3 - 1, /2 = 1) 

V oTT 

|^(-1)'^+^ V'I^+172 (for h = /3, ^2 = 1) 

|-(-l)^^+V^3 (for/i = /3 + l,/2 = l) . (C.26) 

OTT 



The Wigner-9j symbols are calculated as 
h h h 

I4 Irj 1% 

1 1 2 



^ 2(/3±l) + l f /i /4 1 1/ ^2 /si..., , , 

= \/ ^ [ ;3±2 /3±1 /5 ri /3±1 h h ^ (fo^^^ = ^3±2) 



(2/3-l)(2/3 + 2)(2/3 + 3) r /i 1 1 r /2 /5 1 



30(2/3)(2/3 + 1) ( /3 /3-I /5ji/3-l /3 h 



+ 1 



+ i 



/2 (2/3~l)(2/3 + l)(2/3 + 3) r /i k 1 1 f ^2 /5 1 
15(2/,,)(2/, + 2) \ /3 h k j \ h h h 

' (2/3-i)(2/3)(2/3 + 3) r /i /4 11 r /2 /5 1 

30(2/3 + l)(2/3 + 2) \ /3 /3 + I k }\h + l h h 

(for k = /3) , (C.27) 



C.4 Analytic expressions of the Wigner symbols 
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where these Wigner- 6j symbols are analytically given by 

|/l /2 1 1 /i;^^^ (f0rZ2 = /l-l,/5 = /4 + l) 

I4 h h ) y 2I4+3P3 2Zi+l-r3 

^_^y^+l,+l,+, l 2{h + k + k + 2){h + U-k + l) 
y 2U+3-f3 

l{-ll+U + k + l)ih-U + k) ,r J , , , , -.x 

X\ ^ {ioi l2 = ll,h = k + l) 

y 2Zi+2-r3 

^ ^_-^y,+U+le+l -h+lMlP2 h-U+l,+lP2 (foj.Z2 = i, + l,i5 = i4 + l) 

y 2/4+3-'3 2/i+3-f3 

= (_l)'l+U+;6 + l ^^l^ + 1) + _ + 1) _ + 1) _ ^^(^^ + 

^ / 2(/i+/4 + /6 + l)(/l + /4-/6) (f^ 1 -1 1 7 _ 7 ^ 

''VHl + ^4 + /6 + l)(/l-/4 + /6)2U+2P3 2,+lP3 " " 1' " 

^y,+^+;,+l /4(^4 + 1) + + l)(/4 + 1) - kik + 1) - + l)k 

\/ 2k+2-P3 2ii+2-P3 

(for l2 = li,k = U) 

= + 1) + (/i + + 2)(/4 + 1) - /6(/6 + 1) - + 1)^4] 

. / 2(-/i + /4 + /6)(^l-^4 + /6 + l) ^f. 7 _ 7 ^ 1 7 _ 7 ^ 

(C.28) 

Using these analytical formulas, one can reduce the time cost involved with calculating the CMB 
bispectrum from PMFs. 



164 



D POLARIZATION VECTOR AND TENSOR 



D Polarization vector and tensor 



We summarize the relations and properties of a divergenceless polarization vector ei'^^'' and a 



transverse and traceless polarization tensor e^^^'' j6 



11 



The polarization vector with respect to a unit vector n is expressed using two unit vectors 9 
and perpendicular to n as 



1 



V2 



[^,(n)±^4(n)] . 



(D.l) 



This satisfies the relations: 



n'^ei±i)(n) = 0, 

Va' (for A,A' = ±1) . 



(D.2) 



By defining a rotational matrix, which transforms a unit vector parallel to the z axis, namely z, 
to n, as 



cos On cos 0„ — sin (pn sin On cos (pn 
S{n) = I cos6'„sin0„ cos0„ sin6'„sin0 
— sin On cos On 



(D.3) 



we specify and as 



(n) = S'(n)x, <j){ii) = S{Yi)y 



(D.4) 



where x and y are unit vectors parallel to x- and y-axes. By using Eq. (D.l), the polarization 
tensor is constructed as 



n) = v^ei±^)(n)e 



(±1), 



n 



(D.5) 



To utilize the polarization vector and tensor in the calculation of this thesis, we need to expand 



Eqs. (D.l) and (D.5) with spin spherical harmonics. An arbitrary unit vector is expanded with 



the spin-0 spherical harmonics as 




-m((5„,i + 5m -i) 

V^Snifl 



(D.6) 



Here, note that the repeat of the index implies the summation. The scalar product of a™ is 
calculated as 



m m / 1 "l^A r\/''^n S 



(D.7) 



Through the substitution of Eq. (D.4) into Eq. (D.6), is expanded as 



Un) = $^«-Fi,„(^(n)) = 5^a-5^Di;i:,(5(n))rw(i) 
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= + J]OFi„(n) . (D.8) 

* m 

Here, we use the properties of the Wigner D-matrix as described in Appendix [B] |3| [5| [6| [Ts] 



Air 



2i+l 



1/2 



(D.9) 



:-l)^_.F;^(n) . 



In the same manner, 6 is also calculated as 



(D.IO) 



hence, the explicit form of Eq. (D.l) is calculated as 



(D.ll) 



Substituting this into Eq. (D.5) and using the relations of Appendix C and /2'i^/^^ 
polarization tensor can also be expressed as 



a 



rrib 



Mrriamf, 



2 11 

M TUa nib 



This obeys the relations: 



,(±2)* 



(±2), 



-n 



eii)(n)elf (n)=2VA' (for A, A' = ±2) 
Using the projection operators as 



^(0) ik-x 
Q(±l)g.k.x 



5aA 



(±2) ik-x _ ^(±2)/f N ik-x 



the arbitrary scalar, vector and tensor are decomposed into the helicity states as 
r/(k) = r/(0)(k), 

A=±l 



^ , the 



20f 



(D.12) 



(D.13) 



(D.14) 



(D.15) 
(D.16) 
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A=±l A=±2 



ab 



(D.17) 



Then, using Eq. (D.2) and (D.13), we can find the inverse formulae as 



a; 



(0) 



(±1) 



(k) 
(k) 



X(°)(k) 
X(^^)(k) 

x(^^)(k) 



^Oi?(k)x..(k) , 

^oir^)(k)x..(k) . 



From these, we can derive the relations of several projection operators as 



(D.18) 
(D.19) 

(D.20) 
(D.21) 
(D.22) 



oi?(k) 



Oif^(k) = er(k) 



oifHk)ne(k) 



-kah + -Sab 



.(±l)^:^ , L.(±i), 



a. 



2 1 1 

M uia nil, 



E «i^v(t)"r<- ( I, ,1 ) 



Mmaini, 



,(±2), 



Pab(k) = 6ab - kah 



M rria nih 



2 E ^Lu' E ^LM(k)«r«r ( 

L=0,2 Mmamt, ^ 



L 1 1 

M rua ruh 



Oi",^(k)n,(k) = -Pacik) , 

Oit'^(k)nc(k) = fc.6(±^)(k) 



.(±2) 



(k) 



r^''''kae^\k) = T^tf'\k) . 



ir'e^^'\k)e[-'\k) , 
(±1), 



(D.23) 
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E Calculation of f^l and 



Here, we calculate each product between th e wave number vectors and the polarization tensors of 
f^l and f~^^ mentioned in Sec. 



8.2.1 



14 



Using the relations discussed in Appendix Dl the all terms of /i^g are written as 



Vl/3 



(-Ai) {-A2) (-A3) 
«i jk ki 



{-AO (-A2) {-A3)r 



(_Ai) (-A2) (-As),- £ 
^ij ^7 «^2«fc3fc 



(-Ai) (-A2) {-A3)r £ 
^ij e^,, 



2/«'3j 



X 



[Snf A.1^2M(kl)A.V'2V(k2)A3>^2V(k3) 

M,M',M" 

^ /7 / 2 2 2 
10 V 3 V ^ ^' M" 

3/2 



X 



X 



(Stt 

15 
(Stt 

47r 

y 

-(Stt 

4:71 



X 



X 



Ai^2A/(kl)A2^L'A/'(k2)A3^L"A/"(k3) 

L',L"=2,3M,M',M" 



_-|nL' rAaO-Aa rA30-A3 
-■-J -'l'12 -'l"12 



2 L' L" 
2 1 1 



2 L' L" 
M M' M" 

'^^^ X] X] Ai'i^2A/(kl)A2^L'Af'(k2)A3^L"A/"(k3 

L',L"=2,3 M,M',M" 

' 2 L' L" 

L' ,-A20-A2jA30-A3 



-^J -'L'12 -'L"12 



2 L' L" 
M M' M" 



1 1 2 
1 2 1 



3/2 



A2 ^L' A/' ( 1^2 ) A3 '5^L" A/" ( ks ) 

L',L"=2,3 Af,A/',A/" 

2 L' L" 
M M' M" 



L' 7-A2O-A2JA3O-A3 



12 -'L"12 





1 






V L" \ 


1= 


1 






1 1 J 



(E.l) 



In the calculation of f^^2^ need to consider the dependence of the tensor contractions 

on 7]^^''. Making use of the relation: 



1 1 1 

nia rUb rric 



(E.2) 



the first two terms of f~ ^ reduce to 




^^^^■'efc'^^ej.-^^^i'^^)^™ = -(87r)3/2 ^ ^ A.r4,(k;)A2>^2V(k2)A3>^£"M"(k3) 

L"=2,3 Af,Af',Af" 

27r, ,,r».x„n_x, / 2 2 L" \ ( 2 2 L" 
^ M M' M" J \ 1 2 1 

^f^' E E A.1^2M(kl)A2V'2M'(k2) A3^L"Af"(k3) 
L"=2,3 Af,Af',Af" 

L" rAaO-Aa / 2 2 ^" ^ J 1 1 ^ 



,Vifcp(-^l)p(-^2)p{-A3)L" 
'''/ "^feij ^mi ^m<j "'3j 



X 2V27r(-l)^ /^?; 



12 



(E.3) 
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For the other terms, by using the relation 



we have 



L',L"=2,3 M,M',M" 

y ^(_\\L" 7-A20-A2 T-AgO- As 
^ 2 V -^J -'L'12 -'L"12 



X 



2 L' L" 
M M' M" 




2 

L',L"=2,3 M,M',M" 
^ Ai ^2*M (kl ) A2 ^L'M' 0^2) X3 Yl"M" (^s) 

27r 
X — 
15 



X 



1^^'"7-A20-A2 7-A3O-A3 

2 L' L" \ ( 2 L' L" 
M M' M" I \ 1 2 2 



(E.5) 
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F Graviton non-Gaussianity from the Weyl cubic terms 

Here, let us derive the bispectra of gravitons coming from the parity-even and parity-odd Weyl 



cubic terms, namely, Eqs. (8.16) and (8.17) |14|. For convenience, we decompose the interaction 
Hamiltonians of and WW'^ {S.lbn into 



TT _ rrin) 
J^mt — / ^ J^int • 

n=l 

Depending on this, we also split the graviton non-Gaussianity as 



(F.l) 



\n=l I int rn=l \n=l 

In what follows, we shall show the computation of each fraction. 



F.l 



(m) 



int 



(F.2) 



The bracket part of Eq. (8.12) in terms of is expanded as 
^<3(0:,fl7^'"Hk„,r) 

n=l 



n=l 



T \ 1 



r. / 4 



n=l 



X 



: 



. n=l 



-: 



-A" 

x6 
3 



x2ilm 



n / E 

J lm=l 
m=l J l.n=l 

(Hrr {^yii^-r^ (Ekn) eSf Vki)4^Vk;)ei^Vk3) 
Ylilds - kl'jds){kn,T')Ydsikn,r)^ - ^Yllds{kn,T){^*as - klYds)ikn,r') 

YlildS - kl-ids){kn, ^'h*dsikn, t) 



.n=l 

-2/ rr'\2 



2 ' \T, 



n=l 



(F.3) 
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Here, we use 



n 

n=l 



„(A;)„(A„)t 



) = (2x)'5(ki + kyij,,^, i(k'i + k3)5j;,j,5(k2 + ki)ij„j. + 6 perms. 



n 

n=l 

el7'^(k) = e^)(-k). 



(F.4) 



Furthermore, since 



7d5 - A;^7(i5 



2ii'r' 



^3/2g-ifer' 



pl 



(F.5) 



n=l 



the time integral at r — > is performed as 

i\ A 3 



Im 



kfklkllm (fl^*,sikn,r) 



8H^ 



r-)-0 8if^ 
)■ — TrRe 

Ml, 



\n=\ 



(F.6) 



where kt = kn- Thus, the graviton non-Gaussianity in the late time limit arising from 

is 



(1) 



\n=l 



.71=1 



pl 



A 



Re 



.-A 



(F.7) 



The bracket part in terms of H^i is given by 



n=l 



= 



n=l 
A 



0-0 



n=l 



- ^^-'iHr'f i-\ k,k,{2^f5 ( 5^kJ 4.4,eS7^^)(k;)ei7^^Hk,)ei7^^)(k3) 



vn=l 



x2ilm 



{idS - kl-ids){kuT')'^ds{k2,T')'^ds{k^,r')Y\^'^ls{kn,T) 



n=l 



+ 5 perms. . (F.8) 



F.2 WW^ 
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Using 



pi 



(F.9) 



we can reduce the time integral to 



Im 



r dr'iHr'f (-\ fc2fc3(7d5-^h<i5)(A:i,T')7d5(A;2,T')7d5(A;3,T')n7d5(^n,T) 



>■ — TTR-R-e 







J-tb+A-iktr' 



oo 



(F.IO) 



and obtain 

3 



3 \ / 3 \ 6 



A 



Re 



x^A;2,e;;'^)(ki)A;3,e^7'^^^(k2)e^7^^^(k3) + 5 perms. . 



(F.ll) 



The graviton non- Gaussianities from H'^', and i/^^i are derived in the same manner as that from 



(4) 



(2) 

4 / 3 



E n7'^"'(k: 



m=3 \n=l 



X 



(-As)/ 



2i 



^3,e(7^^)(ki)4"'^Hk2)e^7'^^^(k3)fc2, + 5 perms. 



(-As), 



(F.12) 



F.2 WW'^ 



At first, we shall focus on the contribution of H~ The bracket part is computed as 

WW-' 



n=l 



n=l 



.n=l 



- j d^x'K-^Hr'fy-] (-3) 



n / He-^"' E 

xr^^^'^elf (k;)efi) (k^eg^^ (k^,) {ik'.J 

X [{idS - k'^lds) {k[, r') (7<i5 - k'^'^ds) ik'2, r')^asik's, r') 



X : ( 



(An)t 



. m=l 



.n=l 
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X : ( 



n=l 
A 



. m=l 



: 



.n=l 



X 



{idS - kl-fds) {h, t') {^ds - kl-fds) {h, T')^ds{k3, r') < 7*ds{kn, r) 



,n=l 



- {%s - khds) {k,,T') {%s - khls) {k^yWdsik.y) Y[lds{kr.,T) 

ln=l 

+5 perms. 

{idS - kl^ds) (ki, t') {%s - kj^ds) {k2, T')^ds{k3, r') < ^dsikn, r 



x2ilm 

+5 perms. . 

Via the time integral: 



.n=l 



(F.13) 



Im 



dr'iHr'f (^^^ ks {^ds - k^ds) {ki,T') {%s - khds) {k2, r') 

Xids{k3,T') \YlYdsikn:r) 



,n=l 



(F.14) 



we have 

3 



(1) 



\n=l 



.n=l 



Im 



x(-3i)[r?^^-^ei;'^)(k;)e(;;,'^Hk2)eS;'^)(k3)4^ + 5perms.] . (F.15) 



Like this, we can gain the second counterpart: 

3 \ (2) 



n=l I WW^ \n=l / \ Pi/ \ 



A 



Im 



'^'^■'ei;'^^(k;)eL-/^)(k;)eL-/^)(k3)4, + 5 perms.l . (F.16) 



The bracket part with respect to H~ , is 

WW-' 



n=l 



F.2 WW^ 
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n=l 



T 



T 



n=l 



A' =±2 



3 



: 



: 



(Am)t 



, n=l 



, m=l 



. m=l 



,n=l 



^ * ' \n=l 



X 



3 



n^^5(^n,r')7ds(/cn,r) 



,n=l 



,n=l 



5 perms. 



vn=l 



■ 3 

x2ilm Yl'yds{kn:r')-fas{kn:r) + Sperms.. 



n=l 



(F.17) 



The time integral is 



Im 



J -00 



> — ^im 



(F.18) 



so that the bispectrum of gravitons becomes 

i \ (3) / 3 



\n=l 



vn=l 



pi 



A 



Im 



r/^^-^e|,/^)(ki)e(-^^)(k2)e;-'^\k3)fcu4/4,n + 5 perms.] . (F.19) 
Through the same procedure, the bispectrum from H~^^ is estimated as 



n^^'-Hkn) 



\n=l 



(4) 



pi 



A 



Im 
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^3)kikhik3i + 5 perms. 
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